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Water  Conservation  Area  1  is  a  57,234  ha  remnant  of  the 
Northern  Everglades,  presently  used  as  both  a  flood  control 
retention  area  and  a  wildlife  refuge.   It  is  a  unique  wetland 
ecosystem  in  which  the  dynamics  of  water,  vegetation,  fire,  and 
peat  are  tightly  linked.   These  interacting  factors  produce  a 
landscape  pattern  that  provides  excellent  habitat  for  wading 
birds  and  other  Everglades  wildlife. 

This  doctoral  research  focused  on  the  landscape  of  Water 
Conservation  Area  1  and  the  processes  determining  its  pattern. 
Few  ecological  studies  had  been  conducted  in  Water  Conservation 
Area  1,  so  previous  studies  of  the  Everglades  and  the  Okefenokee 
Swamp  were  reviewed  for  relevant  data.   The  landscape  was 
observed  and  described  in  the  field.   Aerial  photography, 
satellite  images,  and  a  computer  geographic  information  system 
(GIS)  were  used  to  quantify  the  landscape. 
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The  GIS  was  also  used  to  perforin  a  spatial  analysis  of 
tree  islands  in  Water  Conservation  Area  1.   The  analysis  showed 
more  small  islands  than  large  ones,  and  a  steep  increase  in  the 
sizes  of  islands  above  a  certain  size.   Simple  models  built  to 
generate  island  size  and  frequency  distributions  demonstrated 
mathematical  functions  that  could  create  a  size  distribution 
similar  to  that  found  in  Water  Conservation  Area  1.   A 
conceptual  model  proposing  an  explanation  for  the  tree  islands' 
size  distribution  was  developed.   The  theory  proposed  that  tree 
islands  grow  in  size  over  time,  but  fire  reverses  island  growth. 
The  ability  of  fire  to  limit  island  growth  is  reduced  once 
islands  reach  some  size  threshold.   Therefore,  the  growth  rate 
of  large  islands  is  greater  than  that  of  small  islands. 

A  spatial  model  of  vegetation  succession  and  fire  was  made 
to  illustrate  the  conceptual  model  in  two  dimensions,  using  a 
map  of  the  actual  landscape.   The  GIS  was  used  to  construct  and 
run  the  spatial  model  using  a  classified  satellite  image  as 
input.   The  model  produced  maps  of  landscape  pattern  under 
different  fire  frequency  and  severity  regimes.   The  landscape 
composition  and  heterogeneity  of  the  maps  were  compared  with  the 
current  landscape.   The  results  illustrated  the  importance  of 
fire  in  determining  landscape  pattern. 
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CHAPTER  1 
INTRODUCTION 


Twenty-five  years  ago,  as  the  Apollo  8  space  mission 
prepared  to  orbit  the  moon,  many  people  envisioned  that  the 
future  of  humanity  lay  in  outer  space.   For  three  decades  the 
public  imagination  had  entertained  visions  of  a  future  society 
in  space.   It  was  expected  that  the  achievement  of  a  lunar 
landing  would  lead  to  new  efforts  in  space  exploration  and 
colonization.   However,  while  the  television  images  of  Apollo  11 
astronauts  on  the  moon  inspired  awe  and  pride,  another  image  was 
just  as  powerful--the  planet  Earth,  suspended  in  the  void  of 
space.   Along  with  the  astronauts,  viewers  saw  the  beauty  and 
isolation  of  the  earth,  and  felt  a  new  appreciation  of  our 
dependence  upon  it.   Since  that  time,  public  interest  in  space 
colonization  has  waned,  while  the  use  of  space  technology  to 
monitor  our  home  planet  has  increased.   Photos  from  manned  space 
craft  and  images  from  satellites  have  shown  us  the  colors  and 
patterns  of  the  atmosphere,  the  oceans,  and  the  land  as  never 
before . 

It  was  such  an  image  of  Earth  that  inspired  the  work  that 
follows:  a  satellite  image  of  the  Florida  Everglades  showing  a 
remarkable  pattern  of  tree-covered  islands  emerging  from  a 
background  of  marsh.   The  islands  resemble  a  school  of  fish 
swimming  upstream.   A  particularly  striking  example  of  the 
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Figure  1-1.   Satellite  image  of  Water  Conservation  Area  1 


pattern,  contrasting  with  the  pattern  of  nearby  urban 
development,  appears  in  a  northern  remnant  of  the  Everglades 
known  as  Water  Conservation  Area  1 . 

The  islands'  pattern  inspired  many  questions.   How  were 
these  islands  formed?   Why  were  the  islands  oriented  in  the  same 
direction?   Why  were  a  few  islands  large  and  oblong,  while  most 
others  were  smaller  and  round?   Were  the  islands  changing  in 
size  over  time  and  were  they  subject  to  disturbances?   Were  the 
processes  responsible  for  the  islands'  pattern  still  operating, 
and  if  so,  what  effects  would  recent  human  alteration  of  the 
region's  ecology  have  on  them? 

A  review  of  the  ecological  literature  on  Everglades 
islands  showed  that  these  questions  have  not  been  definitively 
answered,  and  are  very  difficult  if  not  impossible  to  answer 
using  standard  ecological  methods.   Island  formation  is  too  slow 
or  infrequent  to  document  in  the  time  frame  of  a  scientific 
study.   Some  of  the  processes  that  formed  and  disturbed  this 
landscape  have  ceased  functioning.   Examination  and  dating  of 
the  peat  record  could  provide  clues  to  the  islands'  histories, 
but  collecting  and  analyzing  an  adequate  sample  size  is 
prohibitively  expensive. 

These  constraints  are  common  in  landscape  ecology.   Yet 
the  "big-picture"  viewpoint  that  remote  sensing  has  given  us, 
and  a  new  emphasis  on  conserving  whole  ecosystems  rather  than 
species,  increasingly  raise  questions  to  be  answered  at  the 
landscape  scale.   New  methods  must  be  developed  to  address  these 
questions.   In  addition,  as  we  seek  to  understand  larger 


evolving  ecosystems,  our  expectation  of  definitive  answers  may 
have  to  be  revised.   We  must  develop  skills  to  extract  useful 
knowledge  from  ecosystems  without  statistical  significance  and 
in  spite  of  uncertainty  and  unpredictability. 

In  an  Indian  folktale,  three  blind  men  attempt  to  identify 
an  elephant.   Because  each  man  examines  only  one  part  of  the 
animal,  they  in  turn  describe  it  as  a  snake,  a  brush,  and  a 
tree.   By  feeling  all  three  parts,  a  fourth  blind  man  correctly 
identifies  the  elephant,  even  though  he  cannot  feel  the  entire 
elephant  at  once.   This  doctoral  work  took  a  similar  approach, 
attempting  to  synthesize  an  understanding  of  the  dynamics  of  the 
landscape  of  Water  Conservation  Area  1  by  piecing  together  data 
from  parts  of  the  system. 

Rather  than  trim  the  subject  of  island  pattern  in  this 
piece  of  the  northern  Everglades  into  a  smaller  question  on  a 
smaller  scale,  this  work  remained  focused  on  the  landscape.   The 
tools  of  remote  sensing  and  a  geographic  information  system 
(GIS)  were  used  to  observe  and  quantitatively  analyze  the 
landscape  (see  Figure  1-2).   Data  from  existing  ecological 
studies,  many  of  which  were  carried  out  at  the  community  level, 
were  synthesized  with  ground  observations  to  form  a  conceptual 
model  of  the  interaction  between  ecosystem  processes  and 
structures.   A  portion  of  this  conceptual  model  was  translated 
into  a  landscape  model  on  the  computer  using  remotely  sensed 
data  as  input.   A  geographic  information  system  ran  the  model, 
simulating  fire  and  vegetation  succession. 
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Figure  1-2.   Diagram  of  approach  to  landscape  analysis 


This  dissertation  describes  the  progression  of  this 
undertaking.   Chapter  1  gives  an  overview  of  previous  studies  of 
landscape  pattern,  and  describes  the  ecosystem  of  Water 
Conservation  Area  1.   Chapter  2  describes  the  collection  of  data 
from  Water  Conservation  Area  1  and  the  processing  of  satellite 
images  used  to  quantify  landscape  pattern. 

Chapter  3  describes  a  spatial  analysis  of  the  tree  islands 
in  Water  Conservation  Area  1.   The  analysis  focused  on  the 
distribution  of  tree  islands  of  different  sizes,  as  recorded  in 


a  1987  classified  satellite  image.   The  size  and  frequency 
distributions  of  tree  islands  raised  three  questions.   Could  the 
ongoing  formation  and  growth  of  islands  create  a  size 
distribution  like  that  found  in  the  landscape  of  Water 
Conservation  Area  1?   Why  were  more  small  islands  than  large 
islands  found  in  the  satellite  image  of  the  landscape?   What 
might  be  causing  a  discontinuity  found  in  the  size  distribution 
of  islands? 

The  first  two  questions  were  addressed  in  chapter  3  with 
simple  models  built  to  generate  island  size  and  frequency 
distributions  for  different  growth  functions.   The  models 
demonstrated  that  the  creation  and  growth  of  islands  over  time 
could  create  a  distribution  like  that  found  in  the  landscape  of 
Water  Conservation  Area  1.   The  third  question  was  not  easily 
answered.   Theories  of  the  relationship  between  ecosystem 
processes  and  ecosystem  architecture  suggested  that  a  change  of 
some  process  affecting  the  system  could  cause  a  discontinuity  in 
the  island  size  distribution. 

Chapter  4  outlines  a  conceptual  model  of  landscape  pattern 
dynamics  that  proposes  a  mechanism  for  the  discontinuity  in  the 
island  size  distribution.   The  conceptual  model  was  also 
converted  into  a  spatial  model  of  vegetation  disturbance  and 
landscape  pattern.   Chapter  4  describes  the  methods  and  results 
of  the  spatial  model,  used  to  simulate  landscapes  under 
different  fire  frequency  and  severity  regimes. 

Chapter  5  summarizes  and  discusses  the  entire  work.   The 
work  produced  both  an  example  of  a  holistic,  landscape  approach 


to  an  ecosystem,  and  a  tool  to  spatially  visualize  a  theory  of 
landscape  dynamics. 

The  Study  of  Landscape  Pattern  Dynamics 

Landscape  Ecology 

The  study  of  the  spatial  distribution  of  vegetation  and 
its  relationship  to  environmental  processes  can  be  carried  out 
at  many  levels  of  resolution:  from  the  study  of  plants  in  a 
vegetation  community,  to  communities  within  a  landscape,  to 
landscapes  on  the  globe.   A  landscape  is  defined  as  area  of  land 
containing  multiple  ecological  communities,  "a  heterogeneous 
land  area  composed  of  a  cluster  of  interacting  ecosystems  that 
is  repeated  in  similar  form  throughout"  (Forman  and  Godron 
1986:594) .   Landscape  ecology  examines  and  analyzes  the  spatial 
patterns  and  dynamics  of  landscapes  over  time. 

One  area  of  study  in  landscape  ecology  is  the  relationship 
between  vegetation  structure  (i.e.,  the  size  of  community 
patches,  patch  connectivity,  edge  length  and  diversity) ,  and 
environmental  disturbance,  (i.e.,  fire,  wind  and  human  activity) 
(Weins  et  al.  1985,  Turner  1987).   The  majority  of  research  on 
this  topic  has  focused  on  the  effect  of  disturbance  on  landscape 
pattern  (Potter  and  Kessell  1985).   For  example,  Suffling  et  al. 
(1988)  investigated  the  control  of  landscape  diversity  by  fire 
disturbance  in  boreal  forests.   They  found  that  forests  with  an 
intermediate  fire  frequency  had  a  higher  landscape  diversity 
than  forests  with  high  or  low  fire  frequencies. 


Fewer  researchers  have  addressed  the  reciprocal  effect  of 
landscape  pattern  on  disturbance  (Franklin  and  Forman  1987). 
Turner  et  al.  (1989)  used  spatial  percolation  theory  models  on 
simple  abstract  landscapes  to  study  how  landscape  pattern 
affects  the  spread  of  disturbance.   They  found  that  the 
frequency  of  disturbance  (the  probability  of  a  disturbance  event 
occurring  in  a  given  time  period)  and  the  intensity  of 
disturbance  (the  probability  of  a  disturbance  spreading) 
influence  the  area  of  landscape  disturbed.   However,  this  effect 
differs  according  to  the  amount  and  connectivity  of  patches  of 
habitat  susceptible  to  disturbance.   Disturbance  and  landscape 
structure,  therefore,  have  a  complex  two-way  interaction. 
Simulation  modeling  is  required  to  gain  an  understanding  of 
these  interactions,  due  to  the  complexity  of  landscape  dynamics. 

Models  of  Landscape  Pattern 

The  study  of  landscape  dynamics  is  difficult  because  of 
the  large  scale  in  space  and  time  over  which  landscapes  change. 
Scientific  methods  for  the  study  of  single  species  or 
communities  within  ecosystems  are  well  developed,  while  there 
are  few  methods  designed  for  landscapes.   The  amount  of 
information  contained  in  a  large  natural  system  is  too  great  to 
be  analyzed  using  standard  ecological  sampling  and  analysis 
techniques.   Therefore  landscapes  are  frequently  studied  by 
simplifying  them  into  models,  such  as  ecosystem/energy  models 
(Odum  1983),  vegetation/forest  stand  models  (Pearlstine  et  al. 
1985,  Menaut  et  al.  1990),  patch  dynamics  (gap)  models  (Shugart 


and  Noble  1981,  Shugart  and  West  1977)  and  gradient  models 
(Kessel  1979) . 

The  term  "landscape  model"  is  sometimes  used  to  refer  to 
models  which  simulate  processes  that  determine  where  entire 
landscapes  of  different  types  are  found.   These  models  operate 
at  a  higher  hierarchical  level  than  the  models  referred  to  in 
this  work.   "Landscape  model"  is  used  here  to  indicate  models 
which  simulate  dynamics  within  a  landscape. 
Appropriate  scaling  in  models 

Ecological  computer  models  can  be  packed  with  information 
and  simulate  a  great  number  of  processes,  but  this  complexity 
does  not  necessarily  produce  realistic  behavior  or 
understandable  results.   A  model  which  limits  information  and 
processes  to  the  temporal  and  spatial  scales  pertinent  to  the 
model's  subject  is  more  comprehensible  and  efficient  than  one 
which  attempts  to  explicitly  model  as  many  processes  as 
possible.   In  general,  models  are  best  restricted  to  variables 
ranging  over  no  more  than  three  orders  of  magnitude  in  space  or 
time  (Gunderson  1992) .   This  limits  a  model  to  variables  within 
one  or  two  scale  domains. 

The  scale  of  a  variable  or  process  can  be  identified  by 
its  frequency,  by  its  area  of  influence,  or  by  its  energy. 
Areas  best  measured  in  centimeters  to  tens  of  meters,  and  time 
frequencies  of  minutes  to  decades,  fall  into  the  domain  termed 
the  "microscale . "   Areas  of  tens  of  meters  to  hundreds  of  meters 
and  time  frequencies  of  decades  to  centuries  fall  into  the 
"mesoscale"  domain.   Areas  of  hundreds  of  meters  to  thousands  of 
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kilometers  and  time  frequencies  of  centuries  to  millennia  fall 
into  the  "macroscale"  domain  (Holling  1992a). 

In  the  Everglades,  for  example,  the  species  composition  of 
a  wet  prairie  community  changes  seasonally  during  an  annual 
cycle  (Kushlan  1990)  at  the  microscale.   However,  this 
fluctuation  is  not  indicative  of  long-term  community  change  at 
the  landscape  level,  because  the  species  involved  are  all 
normally  members  of  the  wet  prairie  community.   To  add  this 
fluctuation  to  a  landscape  pattern  model  would  unnecessarily 
increase  model  complexity.   Instead,  community  types  rather  than 
species  can  be  used  in  the  model,  understanding  that  species 
composition  varies  within  a  community. 

Landscape  patterning  processes  in  the  Everglades  operate 
at  the  mesoscale  (Gunderson  1992).   These  processes  include 
vegetation  succession  and  disturbances  like  fire,  droughts  and 
storms  (DeAngelis  and  White  1994) .   A  landscape  pattern  model 
should  focus  on  changes  at  the  mesoscale.   It  should  track  a 
measure  of  pattern  or  other  landscape  attribute  which  reflects 
these  processes. 

At  the  macroscale,  slow  climatic  and  geologic  processes, 
such  as  sea  level  rise  and  geologic  uplift,  determine  where 
entire  landscape  types  are  found  within  the  Everglades  (Davis  et 
al.  1994).   Rather  than  explicitly  include  macroscale  processes 
in  a  landscape  pattern  model,  these  slow  global  changes  can  be 
assumed  static  for  the  time  period  of  the  model. 
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Spatial  models 

Most  ecological  models  simulate  population  dynamics  or 
ecosystem  functions:  the  flows  of  energy,  materials,  species,  or 
individuals  in  the  system.   The  majority  of  these  models  are 
spatially  aggregated.   This  means  that  while  they  calculate  the 
quantities  and  flows  of  variables  within  a  system,  they  ignore 
the  location  of  those  variables  (Sklar  and  Costanza  1991,  Turner 
and  Dale  1991) .   Spatial  models  differ  from  spatially  aggregated 
models  in  that  they  explicitly  incorporate  the  influence  of 
spatial  structure  into  model  functions,  usually  by  means  of  a 
grid.   While  including  spatial  information  can  mean  the 
sacrifice  of  functional  details,  this  is  justifiable  when 
spatial  factors  such  as  edge  effects  are  critical.   A  landscape 
model  must  be  spatial  if  landscape  pattern  is  to  be  analyzed. 

Individual  cell  models.   One  approach  to  spatial  modeling 
is  to  create  a  function  model  for  each  individual  cell  within  a 
spatially-referenced  grid  (Richardson  1988).   Each  individual 
cell  model  uses  information  from  neighboring  cells  to  simulate 
flows  between  it  and  other  cells  and  to  change  the  value  of  what 
is  stored  within  the  cell.   Rates  of  change  are  therefore  a 
function  of  spatial  relationships. 

The  cell  model  method  can  perform  complex  computations  for 
each  cell  and  handles  continuous  variables  well.   It  is  often 
used  in  models  in  which  an  important  component  is  the  simulation 
of  water  flow  (Costanza  et  al.  1990,  Walters  et  al.  1992, 
Fennema  et  al.  1994) .   However,  the  computing  requirements  to 
run  a  model  of  this  type  are  very  high  if  numerous  cells  are 
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used.   For  example,  a  cell  model  of  the  Atchafalaya  delta 
(Costanza  et  al.  1990)  covered  2479  km^ ,  and  used  a  1  km^ 
gridcell.   A  typical  22-year  run  of  the  model  required 
approximately  24  hours  on  a  VAX  11/780  computer,  or  15  minutes 
on  a  CRAY  X/MP  supercomputer.   If  the  model  had  instead  used  a 
10  m^  gridcell  capable  of  showing  landscape  pattern  but  was 
otherwise  unchanged,  the  same  run  would  take  approximately  27 
years  on  the  VAX  and  104  days  on  the  CRAY  (Costanza  et  al . 
1990)  . 

The  number  of  gridcells  in  a  model  can  be  reduced  by 
modeling  small  areas  (less  than  20  km^ )  or  using  a  coarse 
resolution  (greater  than  1  km^ ) .   However,  if  a  natural 
landscape  is  being  modeled,  a  small  area  may  not  adequately 
illustrate  landscape  behavior.   A  course  resolution  may  not 
capture  spatial  variation  in  the  environment  necessary  to  model 
processes  important  in  determining  landscape  pattern.   A  cell 
model  of  the  hydrology  of  "trie  Everglades  using  a  4  km^  cell 
(Walters  et  al.  1992)  simulated  hydrologic  behavior  well,  but  a 
vegetation  component  added  to  the  model  did  not  predict 
vegetation  change  satisfactorily  (Gunderson  1992). 

Cartographic  models.   Another  variety  of  spatial  model 
uses  a  geographic  information  system  (GIS) .   The  GIS  approach 
emerged  from  the  landscape  architecture  and  geography 
disciplines.   Some  GIS  models  use  vector  (line)  maps  as  input. 
However,  the  majority  of  ecological  GIS  models  use  sets  of 
georef erenced  gridcells  known  as  rasters  or  "data  layers."   This 
approach  is  also  known  as  cartographic  modeling  (Tomlin  1990)  . 
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Rather  than  using  a  single  grid,  with  a  model  and  multiple 
storages  for  each  cell,  cartographic  models  store  different 
types  of  values  in  separate  layers  and  use  a  single  model  to 
superimpose  and  manipulate  them.   These  models  can  perform 
spatial  functions  on  values  associated  with  a  gridcell  and  its 
neighboring  cells  within  a  layer,  such  as  calculating  proximity 
and  diversity  measures,  to  generate  new  layers.   Multiple  layers 
are  combined  by  the  logic  statements  and  mathematical  equations 
of  a  model,  and  the  resulting  information  forms  another  layer. 
This  process  can  be  repeated  to  simulate  changes  over  time. 

GIS  models  have  been  most  frequently  used  to  simulate 
single  events,  with  simpler  calculations  and/or  fewer  iterations 
than  applications  of  the  individual  cell  model  method  described 
above.   However,  GIS  software  programs  have  been  rapidly 
evolving  and  improving,  and  now  many  have  the  capability  to 
iteratively  model  a  landscape  over  time.   The  spatial  model  of 
Water  Conservation  Area  1  developed  as  part  of  this  work 
(described  in  Chapter  4)  is  a  GIS  model. 

The  Landscape  of  Water  Conservation  Area  1 

Water  Conservation  Area  1  lies  in  Palm  Beach  County, 
Florida,  and  consists  of  57,234  hectares  of  peat  wetlands  (see 
Figure  1-3).   The  majority  of  these  wetlands  are  distributed  in 
a  complex  mosaic  of  sloughs,  wet  prairie  and  sawgrass  (Cladium 
jamaicense)    marshes,  shrub  patches,  and  forested  islands. 
Disturbances,  particularly  floods  and  fires,  have  shaped  this 
mosaic  of  vegetation  communities  (Gunderson  1994,  DeAngelis 


14 


1994).   The  juxtaposition  of  wetlands  has  made  Water 
Conservation  Area  1  a  high  quality  wildlife  habitat, 
particularly  for  small  fish,  amphibians,  alligators  {Alligator 
misslssippiensis) ,  and  wading  birds  such  as  herons,  egrets,  ibis 
and  wood  storks  (Mycteria    americana)     (Thompson  1972).   Hoffman 
et  al.  (1994)  identified  this  wetland  mosaic  as  the  "Loxahatchee 
tree  island  marsh",  and  found  that  it  supported  the  highest 
density  of  wading  birds  of  all  the  water  conservation  areas  of 
the  Everglades . 

The  pattern  of  the  landscape  mosaic  in  Water  Conservation 
Area  1  resembles  other  peatlands  in  the  Everglades  and 
elsewhere,  but  is  unique  in  its  abundance  of  small  forested 
islands  known  as  bayheads  or  tree  islands.   Tree  islands 
contrast  with  the  marsh  vegetation  communities  around  them,  so 
they  highlight  the  pattern  of  the  landscape  mosaic  well. 
Because  they  are  easily  delineated  on  aerial  photography  and 
satellite  imagery,  tree  islands  are  useful  indicators  with  which 
to  quantify  the  landscape  pattern. 

The  following  sections  review  the  environmental  factors 
that  created  and  continue  to  influence  the  landscape  of  Water 
Conservation  Area  1  and  its  mosaic  pattern.   These  factors  are 
both  natural  and  human  in  origin. 

The  History  of  Water  Conservation  Area  1 

A  century  ago,  Water  Conservation  Area  1  was  known  as  the 
"Hillsborough  Lake",  although  it  was  not  a  lake.   It  was  part  of 
an  uninterrupted  expanse  of  wetlands  extending  from  Lake 
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Okeechobee  to  the  southern  tip  of  Florida  (see  Figure  1-4) . 
These  wetlands  — the  Everglades— were  geologically  young,  and  had 
been  in  a  state  of  flux  since  their  formation.   They  originated 
approximately  5000  years  ago,  after  glaciers  retreated  and  sea 
level  rose  at  the  beginning  of  the  Holocene  epoch  (Gleason  and 
Stone  1994).   The  southern  Florida  peninsula  had  probably  been 
an  oak  savanna  at  the  end  of  the  Pleistocene  epoch,  but  as  sea 
level  rose,  fresh  water  backed  up  into  low  lying  areas. 
Hillsborough  Lake  was  one  such  area.   Wetland  vegetation  grew  in 
these  new  wetlands,  and  peat  formed,  displacing  water  and 
increasing  the  flooded  area.   Sea  level  continued  to  rise,  and 
eventually  the  wetlands  spread  to  their  current  extent  (Wanless 
et  al.  1994,  Gleason  and  Stone  1994). 

The  Hillsborough  Lake  area  was  first  affected  by  human 
hydrologic  manipulations  in  the  1880' s,  when  the  Caloosahatchee 
River  canal  linked  Lake  Okeechobee  to  the  Gulf  coast,  reducing 
the  volume  of  water  flowing  south  to  the  Everglades  (see  Figure 
1-3).   Flooding  duration  (hydroperiod)  in  the  Hillsborough  Lake 
area  was  reduced  (Thompson  1970) .   In  1915  the  West  Palm  Beach 
and  Hillsboro  canals  were  dug  to  the  north  and  south  of 
Hillsborough  Lake  (Figure  1-3).   These  canals  further  altered 
the  volume  and  flow  pattern  of  water  in  the  area. 

Drastic  hydrologic  changes  came  to  Hillsborough  Lake  in 
the  1950 's,  when  flood  control  work  began  to  enclose  the  area 
with  levees,  canals  and  pump  stations.   Impoundment  was 
completed  in  1960,  creating  Water  Conservation  Area  1.   The 
South  Florida  Water  Management  District  began  to  use  the  area  to 
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store  runoff  from  adjacent  agricultural  lands.   Since  1951,  the 
U.S.  Fish  and  Wildlife  Service  also  has  administered  Water 
Conservation  Area  1  as  a  migratory  bird  refuge  (Light  and  Dineen 
1994,  VanArman  et  al.  1984)  . 


Figure  1-3.   Map  of  South  Florida  showing  the  Water  Conservation 
Areas 
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Figure  1-4 .   Map  of  the  historic  Everglades  system 
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Today  flood  control  structures  restrict  the  hydrologic 
connectivity  of  Water  Conservation  Area  1  to  the  rest  of  the 
Everglades  ecosystem.   Flood  control  and  water  management  have 
allowed  rapid  growth  of  the  agriculture  and  urban  development 
now  surrounding  the  area.   It  remains  a  northern  refuge  for 
Everglades  species  threatened  by  habitat  destruction. 

The  Substrate  of  Water  Conservation  Area  1 

Bedrock 

Water  Conservation  Area  1  lies  over  limestone  bedrock 
known  as  the  Fort  Thompson  Formation.   Originally  formed  on  the 
flat  floor  of  a  shallow  sea,  the  Fort  Thompson  Formation  has 
since  been  eroding  and  dissolving.   At  a  fine  scale  the  bedrock 
surface  is  uneven,  etched  with  a  topography  of  depressions  and 
ridges.   At  a  larger  scale,  a  wide  trough  in  the  bedrock 
surface,  known  as  the  Loxahatchee  Channel,  runs  roughly  north- 
south  from  Lake  Okeechobee  to  the  southern  end  of  Water 
Conservation  Area  1  (Jones  1948) .   Water  overflowing  from  Lake 
Okeechobee  followed  the  trough,  which  is  still  evident  in  the 
grain  of  the  landscape  pattern  in  Water  Conservation  Area  1. 
Before  water  management,  the  trough  was  relatively  wetter  than 
surrounding  areas,  inspiring  the  name  "Hillsborough  Lake."   The 
earliest  Everglades  peat  was  formed  in  the  Loxahatchee  trough 
(Gleason  and  Stone  1994). 
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Figure    1-5.       Bedrock   topography   of    the   Northern   Everglades 
Source:     Jones    1948 
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Peat 

On  top  of  the  bedrock  structure  lies  a  layer  of  peat, 
varying  between  124.5  and  457.2  cm  deep  (Richardson  and 
Kitchens,  unpublished  data) .   The  peat  layer  fills  in  and  masks 
the  microtopography  of  the  bedrock  below  it.   Its  surface  slopes 
gradually  from  north  to  south,  dropping  less  than  four  feet  with 
respect  to  sea  level  over  approximately  30  miles  (Richardson  et 
al.  1990) . 

Peat  is  built  up  by  the  deposition  of  dead  vegetation  at  a 
faster  rate  than  decomposition  occurs.   Estimates  of  peat 
accumulation  in  the  Everglades  ranged  from  7.62  to  10.77  cm  per 
century  (Bond  et  al.  1986) .   The  oldest  peat  samples  dated  from 
Water  Conservation  Area  1  are  approximately  4800  years  old 
(Gleason  and  Stone  1994) .   The  peat  found  at  Water  Conservation 
Area  1  is  derived  mainly  from  water  lily,  Nymphaea    odorata , 
(classified  as  Loxahatchee  peat),  sawgrass,  (classified  as 
Everglades  peat),  and  trees,  Ilex    cassine   and  Persea   borbonia, 
(classified  as  Gandy  peat)  (Gleason  et  al.  1980). 

While  the  peat  surface  of  Water  Conservation  Area  1 
appears  flat  when  observed  from  the  air,  its  micro-topography  is 
quite  uneven.   An  observer  on  the  ground  might  stand  high  and 
dry  on  a  tree  island,  but  wade  knee-deep  in  water  in  an  adjacent 
slough.   While  this  elevational  difference  may  not  be 
significant  in  an  upland  ecosystem,  its  ecological  effect  in  a 
wetland  is  amplified  by  flooding.   The  distribution  of 
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vegetation  communities  at  Water  Conservation  Area  1  is 
correlated  with  this  peat  surface  topography  (Pope  1989). 

Several  forces  shape  the  micro-topography  of  peat  in  Water 
Conservation  Area  1.   These  include:  the  buildup  of  peat  by 
vegetation;  the  disturbance  and  transportation  of  peat  by 
flowing  water  and  alligators;  and  the  removal  of  peat  by  fires. 
These  processes  are  discussed  below. 

Water  in  Water  Conservation  Area  1 

The  natural  hydrology 

Before  the  era  of  flood  control,  Everglades  hydrology  was 
driven  by  rainfall  (Parker  1984) .   Periods  of  high  rainfall  in 
the  Kissimmee  watershed  of  Lake  Okeechobee  caused  flooding  in 
the  Kissimmee  River  Basin.   As  a  result,  Lake  Okeechobee 
expanded  westward  into  wetlands  along  Fisheating  Creek.   Water 
from  Okeechobee  also  rose  over  a  natural  berm  along  the  lake's 
south  shore,  and  flowed  south  into  what  is  now  the  Everglades 
Agricultural  Area. 

During  the  annual  wet  season  (June  -  November) ,  water  from 
the  Kissimmee/Lake  Okeechobee  watershed  connected  with  the 
Everglades  water  body  to  create  a  continuous  inclined  hydrologic 
gradient  extending  from  Lake  Okeechobee  to  Florida  Bay  (Parker 
1984,  Fennema  et  al.  1994).   The  entire  Everglades,  with  the 
exception  of  tree  islands,  were  inundated.   Water  followed  this 
gradient  in  a  sheet  flow  across  the  vegetated  landscape, 
inspiring  the  names  "Pa-hay-okee"  (grassy  waters)  and  "river  of 
grass"  (Douglas  1988).   Unlike  a  river,  however,  water  flow  was 
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so  slow  and  evapotranspiration  so  high  that  it  is  unlikely  water 
from  Lake  Okeechobee  ever  reached  Florida  Bay.   In  Water 
Conservation  Area  1,  rainwater  fell,  flowed  slowly  south  along 
the  slight  elevational  gradient,  was  evaporated  or  transpired, 
and  fell  again. 

In  dry  years,  during  the  dry  season  (December  -  May),  the 
hydrologic  gradient  in  the  Everglades  was  broken.   In  Water 
Conservation  Area  1,  water  evaporated  and  transpired  without 
replenishment,  until  the  only  water  remaining  above  the  peat 
surface  was  found  in  sloughs  and  alligator-excavated  ponds. 
Below  the  peat  surface,  water  was  held  like  a  sponge. 

Although  a  relatively  small  amount  of  water  actually 
reached  Water  Conservation  Area  1  via  overland  flow  from  the 
north,  the  role  played  by  "upstream"  water  was  important. 
Because  peat  is  loosely  consolidated  and  can  be  temporarily 
suspended  in  water,  the  sheet  flow  of  water  probably  transported 
bits  of  disturbed  peat  and  vegetation.   Evidence  of  this 
phenomenon  lies  in  the  shape  of  Water  Conservation  Area  l's 
larger  tree  islands  (Gleason  and  Stone  1994).   All  large  tree 
islands  have  an  oblong  shape,  aligned  along  the  north/south 
direction  of  flow.   Sheet  flow  may  have  eroded  the  heads  and 
flanks  of  these  islands,  and/or  deposited  peat  at  the  downstream 
ends.   Smaller  tree  islands  in  Water  Conservation  Area  1  do  not 
show  this  effect  (see  Figure  1-1)  . 

Hurricanes,  tropical  storms  and  prolonged  high  rainfall 
probably  caused  deeper  flooding  in  the  Everglades,  which 
completely  inundated  tree  islands  and  flowed  at  faster  rates 
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than  the  usual  sheet  flow.   The  effects  of  high-energy 
infrequent  flooding  events  like  hurricanes  on  the  landscape  of 
Water  Conservation  Area  1  have  not  been  documented. 
Flood  control  hydrology 

The  advent  of  flood  control  technology  transformed  the 
hydrology  of  Water  Conservation  Area  1  (Richardson  et  al.  1990) . 
Today  the  hydrologic  regime  is  determined  by  a  regulation 
schedule  set  by  the  South  Florida  Water  Management  District  (see 
Figure  1-6).   The  schedule's  goal  is  to  allow  the  storage  of 
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Figure  1-6.   Rainfall  and  water  regulation  schedule  of  Water 

Conservation  Area  1 

Source:  Richardson  et  al .  1990 
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drainage  water  from  agricultural  areas  during  the  wet  season, 
particularly  in  the  event  of  a  hurricane.   Stored  water  is  held 
during  the  dry  season,  and  then  released  at  the  end  of  the  dry 
season  to  allow  room  for  next  year's  storage.   This  pattern  of 
flooding  and  dry-down  opposes  the  region's  natural  hydrology. 

The  construction  of  canals  and  dikes  halted  the  southward 
sheet  flow  of  water  across  the  Everglades.   Since  the 
impoundment  of  Water  Conservation  Area  1,  water  has  moved  from 
the  canals  at  the  perimeter  of  the  area  inward  to  the  interior 
and  been  held  there.   Rather  than  lying  evenly  like  a  blanket 
over  the  landscape,  as  it  did  when  sheet  flow  existed,  water  now 
pools  behind  the  dike  at  the  low  southern  end  of  the  area. 
Water  at  the  northern  end  of  the  area  drains  south,  but  with  no 
water  flowing  from  the  north  to  replace  it,  the  northern  end 
dries.   Therefore  hydroperiods  in  the  south  have  lengthened 
while  hydroperiods  in  the  north  have  shortened  relative  to 
natural  regimes.   The  high  energy  floods  caused  by  hurricanes  no 
longer  occur,  due  to  flood  control  in  and  around  Lake 
Okeechobee . 

The  alteration  of  the  hydrologic  pattern  of  Water 
Conservation  Area  1,  together  with  lower  water  quality,  is 
changing  the  landscape  mosaic.   Studies  of  vegetation  (Thompson 
1972,  Alexander  and  Crook  1975,  Richardson  et  al.  1990)  have 
documented  the  loss  of  sawgrass  and  tree  islands  at  the  southern 
edge  of  the  area,  the  invasion  of  dense  cattails  {Typha    sp.) 
along  the  western  edge,  and  consolidation  of  shrub  patches  at 
the  area's  northern  tip. 
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The  Tree  Islands  of  Water  Conservation  Area  1 

Tree  islands  at  Water  Conservation  Area  1  are  not  found  on 
top  of  limestone  bedrock  pinnacles  as  are  the  tree  islands  of 
the  southern  Everglades.   Rather  they  are  "bayhead"  islands: 
topographic  high  points  made  only  of  the  peat  itself.   These 
islands  are  formed  by  a  more  complex  process  than  the  pinnacle 
islands  farther  south,  an  effect  of  the  deep  peat  substrate 
(Gleason  and  Stone  1994) .   Because  peat  is  a  product  of 
vegetation,  and  the  vegetation  living  on  bayhead  tree  islands 
contributes  to  island  formation,  Gleason  et  al.  (1980)  termed 
the  islands  "unique  geobotanical  landf orms . " 
Tree  island  formation 

The  shapes  of  the  tree  islands  in  Water  Conservation  Area 
1  bear  a  strong  resemblance  to  tree  islands  found  in  boreal 
patterned  peatlands  in  the  Glacial  Lake  Agassiz  region  of 
central  Canada  and  Minnesota  (Glaser  1987).   There,  peat-based 
tree  islands  are  found  in  sedge-dominated  fens.   Like  Water 
Conservation  Area  l's  islands,  they  are  oblong  and  oriented 
parallel  to  the  direction  of  water  flow,  with  rounded  heads  and 
tapering  tails  like  islands  in  the  southern  Everglades. 
However,  these  islands  were  formed  differently  from  those  in 
Water  Conservation  Area  1.   The  boreal  peat  islands  are  forest 
fragments  left  after  minerotrophic  seeps  emerged  on  the  slopes 
of  raised  forested  bogs.   The  seeps  generated  extensive  "water 
tracks",  wide  treeless  fen  channels  dominated  by  sedges  and 
Sphagnum   species  (Glaser  1987).   The  water  tracks  formed  a 
network  similar  to  a  braided  stream,  then  expanded  and  merged, 
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leaving  islands  of  trees  in  the  interstices.   Some  water  tracks 
are  also  patterned  with  ridges  ("strings")  and  sloughs 
("f larks")  oriented  perpendicular  to  the  direction  of  water 
flow . 

Tree  islands  in  Water  Conservation  Area  1  are  formed  by  an 
additive  process  rather  than  the  subtractive  process  of  the 
boreal  peatlands  mentioned  above.   Gleason  and  others  (1980) 
found  radiocarbon  evidence  in  two  peat  cores  of  islands  in  Water 
Conservation  Area  1  that  suggested  the  islands  originated  as 
peat  "batteries"  (floating  rafts  of  peat  and  live  vegetation) . 
In  peat  prairies  of  the  Okefenokee  Swamp,  an  ecosystem  similar 
to  Water  Conservation  Area  1  but  without  the  influence  of  sheet 
flow,  a  connection  between  peat  batteries  and  tree  islands  has 
also  been  documented.   Research  at  Okefenokee  (Cypert  1972;  Rich 
1979;  Duever  and  Riopelle  1983,  Glasser  1986)  indicated  that 
islands  there  are  formed  under  wet  conditions  when  large 
batteries  of  peat  rise  to  the  surface  of  the  water.   Plant 
succession  occurs  on  the  floating  batteries,  further 
consolidating  the  peat  structure  until  woody  vegetation  is 
established . 

While  the  mechanisms  of  peat  battery  formation  at  Water 
Conservation  Area  1  have  not  been  documented,  newly  formed 
batteries  have  been  observed  during  high  water  following  a  very 
dry  period  (W.  M.  Kitchens,  personal  communication).   Trapped 
gasses  (Rich  1979)  or  buoyant  plant  rhizomes  (Gleason  and  Stone 
1994)  may  be  involved  in  floating  the  batteries. 
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Three  types  of  peat  batteries  have  been  described  in  the 
Okefenokee  Swamp:  mats,  bulges  and  aggregates  (Rich  1979).   Mats 
are  sheets  of  peat  as  large  as  20  feet  in  diameter,  which 
dislodge  from  the  substrate  and  rise  to  the  water  surface.   They 
may  or  may  not  completely  detach  from  the  substrate.   If  a  mat 
battery  completely  detaches,  water  can  carry  it  to  a  new 
location.   When  the  water  level  falls,  the  mat  settles  back  down 
to  the  peat  surface,  creating  a  topographic  high  which  becomes 
an  island. 

The  displacement  of  mat  batteries  is  not  common  in  the 
Okefenokee  Swamp,  where  sheet  flow  is  rare.   Thus,  Okefenokee 
tree  islands  are  round  or  irregularly  shaped,  and  are  not 
oriented  parallel  to  water  flow.   They  are  similar  to  small 
islands  in  Water  Conservation  Area  1.   Transportation  of  mat 
batteries  has  been  repeatedly  observed  in  Water  Conservation 
Area  1,  where  the  shapes  of  large  tree  islands  show  the 
influence  of  water  flow  (Baker  1952;   M.  Maffei,  personal 
communication,  Gleason  and  Stone  1994). 

Bulge  batteries  are  sheets  of  peat  which  rise  to  the 
surface  but  remain  entirely  attached  to  the  substrate.   These 
are  common  at  Okefenokee,  but  not  in  Water  Conservation  Area  1. 

Aggregate  batteries  are  accumulations  of  small  floating 
bits  of  peat,  which  are  gathered  together  by  wind  or  water  flow. 
Unlike  mats  and  bulges,  aggregates  appear  gradually  rather  than 
suddenly.   They  occur  in  a  wide  range  of  sizes  and  are  most 
often  found  at  the  edges  of  alligator-excavated  ponds 
("alligator  holes").   They  may  play  an  important  role  in  the 
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expansion  of  existing  islands  in  addition  to  the  creation  of  new 
islands . 

Glasser  (1986)  traced  all  batteries  and  islands  visible  in 
Chase  Prairie,  Okefenokee  Swamp,  using  aerial  photographs 
covering  a  40-year  period.   She  documented  the  formation  of  82 
islands  during  this  period,  in  an  area  of  approximately  30 
square  kilometers.   Of  these  82  islands,  64  derived  from 
aggregate  batteries  at  the  margins  of  alligator  holes.   During 
the  period  of  study,  all  aggregate  battery  islands  persisted  and 
developed  woody  vegetation  in  10-30  years,  even  if  the  original 
alligator  hole  had  disappeared.   After  40  years  the  islands 
supported  trees.   The  remaining  18  islands  Glasser  observed 
derived  from  bulges  and  mats,  but  only  10%  of  these  persisted 
long  enough  to  develop  woody  vegetation.   Most  bulges  and  mats 
subsided  below  the  water  surface. 

The  relationship  between  alligator  holes,  aggregate  mats 
and  tree  islands  has  not  been  investigated  at  Water  Conservation 
Area  1.   Alligators  are  very  numerous  there,  and  alligator  holes 
are  frequently  found  adjacent  to  tree  islands.   The  role  of 
alligators  in  shaping  peat  topography  is  clearly  important  and 
deserves  study,  although  it  is  not  explicitly  addressed  in  this 
work. 
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Fire  in  Water  Conservation  Area  1 

Fire  and  the  Everglades  ecosystem 

A  major  factor  shaping  the  Water  Conservation  Area  1 
landscape  is  fire.  Fire  is  an  omnipresent  process  in  the 
Everglades  ecosystem,  since 

the  climate  and  the  vegetation  of  south  Florida, 
coupled  with  the  area's  elaborate  system  of  drainage 
ditches  and  canals,  have  created  one  of  the  highest 
fire  potentials  in  the  United  States.   Abundant 
summer  rainfall  (even  during  extremely  dry  years)  and 
warm  temperatures  from  May  through  October  promote  a 
rank  buildup  of  fuel  that  becomes  explosive  during 
the  6-  to  8-month  dry  season.  .  .  .  There  are  more 
days  with  lightning  recorded  in  south  Florida  than 
elsewhere  in  the  Nation.   Over  6,000  lightning 
strikes  were  recorded  in  inland  south  Florida  during 
one  6-hour  period  in  the  summer  of  1976.  (Wade  et  al. 
1980:13) 

Seasonal  probability  of  fires 

The  probability  a  fire  will  be  ignited  in  the  Everglades 

varies  seasonally.   In  Everglades  National  Park  most  fires  are 

ignited  by  humans,  in  winter,  spring  and  fall  (Taylor  1981) . 

Humans  have  been  present  and  lighting  fires  in  the  Everglades 

since  its  genesis,  although  human-ignited  fires  were  probably 

less  frequent  before  Europeans  arrived  (Robertson  1953).   Fires 

not  caused  by  humans  are  ignited  by  lightning,  which  is  most 

common  during  the  summer  storm  season.   While  the  two  types  of 

fires  in  Everglades  National  Park  generally  occur  at  different 

times  of  the  year,  both  types  burn  the  greatest  area  in  the 

month  of  May  (Taylor  1981)  because  of  factors  determining  fire 

propagation  and  intensity. 
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Fire  behavior  factors 

Propagation  and  intensity  of  fires  in  the  Everglades,  as 
in  all  wildland  fires,  are  influenced  by  fuel  density  and 
moisture  content,  atmospheric  humidity,  wind  speed  and 
direction,  and  soil  moisture/water  level  (Anderson  1982, 
Rothermel  1983,  Burgan  and  Rothermel  1984,  Andrews  1986,  Andrews 
and  Chases  1989) .   Topographic  factors  such  as  slope  do  not 
influence  Everglades  fires  as  they  do  in  more  variable  terrain, 
but  canals  and  open  water  may  stop  or  direct  fire  movement. 

The  most  important  factor  in  fire  behavior  is  the  amount 
and  moisture  content  of  fine  fuel  found  on  a  site  (Rothermel 
1983) .   Rothermel  defined  fine  fuel  as  dead  or  living  plant 
material  less  than  1/4  inch  in  diameter  that  can  change  moisture 
content  by  a  specified  amount  within  1  hour.   Needles,  leaves, 
grass  and  dead  herbaceous  stems  are  fine  fuel.   Because  the 
moisture  content  of  fine  fuel  changes  rapidly  in  response  to 
meteorological  and  soil  moisture  conditions,  these  conditions 
are  also  important  in  determining  the  probability  and  behavior 
of  a  fire. 

Fire  potential  variables  in  Water  Conservation  Area  1  vary 
daily,  seasonally  and  between  wet  and  dry  years.   During  the  dry 
season  (December  -  May),  evapotranspiration  lowers  the  water 
level,  exposing  peat  in  some  areas.   Moisture  decreases  in  the 
exposed  peat,  in  live  vegetation,  and  in  dead  fuel  of  all  sizes. 
Some  plant  species  senesce  during  the  winter,  increasing  the 
amount  of  dead  fuel.   Rare  winter  frosts  also  kill  vegetation. 
The  optimal  conditions  for  fires  occur  at  the  end  of  the  dry 
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season,  when  water  levels  and  soil  moisture  are  at  their  lowest, 
dead  fuel  is  abundant  and  dry,  and  the  first  summer  storms 
provide  lightning  as  a  source  of  ignition.   Fires  at  this  time 
of  year  propagate  well  and  may  be  moderate  to  severe. 

During  the  wet  season  (June  -  November)  lightning  is 
common  and  ignites  fires,  but  in  most  years  Water  Conservation 
Area  1  is  inundated  with  water  at  this  time.   Fuel  is  moist,  so 
fires  burn  only  a  limited  area  around  the  strike  location  at  a 
low  intensity. 

During  drought  years  both  the  wet  and  dry  seasons  are  dry, 
creating  conditions  for  severe  fires.   Minimal  rain  falls  during 
summer  storms  and  winter  weather  fronts.   Water  evaporates  and 
transpires  without  replenishment  until  it  stands  only  in  sloughs 
and  alligator  holes.   Vegetation  desiccates,  and  the  peat 
surface  itself  may  become  dry  enough  to  burn.   The  absence  of 
standing  water  permits  annual  plants  adapted  to  shorter 
hydroperiods  to  grow  in  the  wet  prairies.   The  resulting  dense 
vegetative  cover  increases  fuel  density.   Fires  ignited  under 
these  dense,  dry  fuel  conditions  are  intense,  and  spread  quickly 
over  large  areas  (Wade  et  al.  1980,  Duever  1984). 
Fire  severity 

The  terms  "fire  intensity"  and  "fire  severity"  are 
sometimes  used  interchangeably.   More  specifically,  intensity 
refers  to  the  energy  expended  by  a  fire,  while  severity  refers 
to  the  effect  of  a  fire  on  the  landscape.   The  severity  of  fires 
in  the  Everglades  landscape  varies,  depending  upon  the  intensity 
of  the  fire.   The  mortality  of  foliage  caused  by  fire  is  a 
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function  of  plant  moisture  content,  the  temperature  of  the  fire 
and  the  time  the  plant  is  exposed  to  the  fire  (Wright  and  Bailey 
1982) .   The  mortality  of  plant  individuals  after  a  fire  is  a 
function  of  many  factors,  including  their  species,  their  age, 
the  season  of  the  fire,  plant  competition  and  the  post-fire 
hydrology. 

Fires  are  low-intensity  during  the  Everglades  wet  season, 
since  plant  moisture  is  high  and  damage  to  living  plant  tissues 
is  reduced.   Live  foliage  above  the  water  surface  may  burn  along 
with  dead  dry  fuel,  but  the  submerged  crowns  of  plants  are  left 
intact.   Everglades  vegetation  is  adapted  to  frequent,  low- 
intensity  fires.   Unless  a  low-intensity  fire  is  followed  by  a 
long  period  of  deep  flooding,  most  species  resprout  quickly 
(Wade  et  al .  1980) .   Low-intensity  fires  hasten  the  release  of 
nutrients  from  dead  biomass  and  slow  the  rate  of  peat 
accumulation  by  consuming  fine  material  which  would  otherwise 
become  peat.   They  maintain,  rather  than  change,  vegetation 
patterns  at  the  landscape  level. 

Moderately  intense  fires  during  the  dry  season  also 
maintain  the  existing  landscape  pattern  by  slowing  vegetation 
succession  and  reducing  peat  buildup.   Moderate  dry-season  fires 
burn  a  greater  proportion  of  aboveground  live  vegetation  than 
wet-season  fires.   The  destruction  of  living  foliage  requires 
plants  to  expend  energy  reserves  to  resprout.   While  sawgrass 
resprouts  quickly  after  a  moderate  fire  (Forthman  1973),  small 
woody  seedlings  with  little  energy  stored  in  their  root  systems 
may  not  survive.   This  serves  to  set  back  the  vegetation 
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succession  of  woody  species.   By  burning  dead  fuel  which  would 
otherwise  be  incorporated  into  the  peat,  moderate  fires  maintain 
the  lower  elevations  and  longer  hydroperiod  conditions  preferred 
by  wet  prairie.   In  the  absence  of  fire,  peat  will  accrete  and 
create  conditions  more  favorable  to  sawgrass  and  brush  (Loveless 
1959,  Wade  et  al.  1980) . 

Extremely  intense  fires,  occurring  during  drought  years 
or  where  water  levels  have  been  artificially  lowered,  cause 
long-term  changes  in  the  vegetation  pattern  of  the  landscape 
(Alexander  and  Crook  1975) .   Intense  fires  burn  all  desiccated 
aboveground  vegetation  susceptible  to  ignition  and  mortality. 
The  high  heat  generated  by  the  consumption  of  the  ample  dry  fuel 
penetrates  into  the  soil,  killing  the  roots  of  some  plants. 
Intense  fires  can  kill  even  sawgrass  and  mature  trees.   They  can 
burn  into  the  peat,  lowering  elevation.   The  lower  elevation 
increases  hydroperiod  once  the  drought  ends,  leading  to  a  long- 
term  shift  in  vegetation  community  (Craighead  1971,  Loope  and 
Urban,  1980) . 
Patterns  of  fire  events 

The  frequency  of  fires  is  important,  because  fires  can 
cause  plant  mortality  if  they  occur  too  frequently,  even  if 
their  intensity  is  low.   Gunderson  and  Snyder  (1994)  analyzed 
fire  records  from  Everglades  National  Park  from  1948  to  1990, 
and  from  Water  Conservation  Areas  2  and  3  from  1980  to  1990  to 
search  for  patterns  of  fire  occurrence.   They  performed  a 
Fourier  analysis  on  the  log  of  fire  sizes  over  time  and  found 
that  the  dominant  temporal  cycle  was  an  annual  one.   A  10  to  14 
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year  cycle  accounted  for  the  next  largest  amount  of  variation. 
This  longer  cycle  "occurs  at  similar  time  ranges  as  other 
environmental  variables  in  the  system;  water  levels,  water 
flows,  and  pan  evaporation  all  fluctuate  on  cycles  of  similar 
time  frames  (Gunderson  1992).   These  variables  indicate  longer 
term  fluctuations  in  wet  or  dry  conditions  and,  hence, 
susceptibility  to  burning."  (Gunderson  and  Snyder  1994:302) 

The  cycles  found  by  the  Fourier  analysis  do  not  indicate 
how  often  a  point  on  the  ground  might  be  expected  to  burn,  and 
neither  do  fire  frequency  (the  mean  number  of  fires  per  time 
period) ,  and  fire  return  interval  (the  mean  time  period  between 
fires) .   The  frequency  at  which  vegetation  at  any  given  location 
is  subjected  to  a  fire  is  more  pertinent  to  the  formation  of 
landscape  pattern  than  the  number  of  fire  events  in  a  region. 
While  fire  frequency  is  often  tabulated  within  a  known  area, 
such  as  a  park  or  management  area,  this  statistic  does  not 
incorporate  the  area  burned  by  each  fire.   Fire  rotation  time  is 
calculated  by  dividing  the  total  area  from  which  fire  data  are 
collected  by  the  mean  area  burned  annually.   Rotation  time 
estimates  how  often  any  point  on  the  ground  is  expected  to  burn, 
or  the  period  required  for  the  entire  area  to  be  burned  if  no 
area  burns  twice.   The  rotation  time  calculated  from  fire 
records  for  Everglades  National  Park  from  1948  to  1979  (Taylor 
1981)  is  33.6  years.   This  includes  wildfires  and  controlled 
burns . 
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Fire  in  Water  Conservation  Area  1 

Detailed  fire  records  for  Water  Conservation  Area  1  do  not 
exist.   Until  recently,  only  very  large  fires  were  noted  by 
wildlife  refuge  managers,  and  only  rough  acreage  estimates  were 
possible.   Some  large  recorded  fires  are  listed  in  Table  1-1 
The  most  severe  fires  burned  in  1962  and  1989.   These  fires 
ignited  in  drought  years  in  which  large  fires  also  burned  in 
Everglades  National  Park  (Gunderson  and  Snyder  1994). 

Table  1-1.   Large  fires  recorded  in  Water  Conservation  Area  1 


Year 

Hectares  Burned 

1955 

4,  047 

1962 

>  40, 470 

1981 

2,  631 

1989 

18, 211 

The  most  recent  severe  fire  in  Water  Conservation  Area  1 
burned  in  April  of  1989,  at  the  end  of  a  dry  season  following 
eight  months  of  drought.   A  winter  freeze  had  heightened  the 
fire  potential  by  increasing  dead  vegetation.   The  fire  ignited 
when  a  controlled  burn  just  outside  the  area  jumped  across  the 
canal  and  started  a  wildfire.   Over  18,000  hectares  (40,000 
acres)  burned  in  the  northern  and  western  portions  of  the  Water 
Conservation  Area  1  (M.  Maffei,  personal  communication) . 

Aerial  photographs  taken  soon  after  the  1989  fire  showed 
that  large  tree  island  edges,  small  and  medium  sized  tree 
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islands,  brush,  sawgrass  and  wet  prairie  were  severely  burned 
(Loxahatchee  National  Wildlife  Refuge,  unpublished  data) .   No 
green  vegetation  was  visible  in  the  burned  area,  with  one 
exception — the  photographs  showed  that  tree  crowns  in  the 
centers  of  large  tree  islands  were  green  and  undamaged. 

It  appears  that  the  1989  fire  in  Water  Conservation  Area  1 
swept  quickly  through  the  prairie  and  sawgrass.   When  it  met 
with  an  island,  the  fire  swept  around  the  island  and  continued 
on,  but  ignited  the  island's  perimeter  in  the  process.   Fire 
then  burned  inwards  on  all  sides  of  the  island.   This  behavior 
is  evident  in  a  panchromatic  SPOT  satellite  image  of  Water 
Conservation  Area  1  taken  while  the  fire  was  burning  (see  Figure 
1-7) .   If  a  tree  island  was  small,  its  vegetation  was  consumed 
by  the  fire.   If  it  was  large,  the  interior  of  the  island  did 
not  burn,  perhaps  because  the  wetter,  sparser  and/or  calmer 
interior  environment  extinguished  the  fire.   This  resulted  in 
the  fire's  energy  being  released  differently  on  large  and  small 
islands.   If  this  phenomenon  occurs  in  all  severe  fires,  then 
fire  is  an  important  factor  determining  the  sizes  of  tree 
islands,  and  thus  the  landscape  pattern  in  Water  Conservation 
Area  1.   (This  hypothesis  is  explored  in  detail  in  Chapter  4) . 
The  interiors  of  large  islands  may  have  escaped  the  fire 
because  their  understory  vegetation  and  microclimate  are  quite 
different  from  other  vegetation  communities  in  Water 
Conservation  Area  1.   The  leaf  area  index  of  a  tree  island 
overstory  is  significantly  higher  than  sawgrass  or  wet  prairie 
vegetation  (Gunderson  1992).   Therefore  tree  island  canopies 
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provide  dense  shade,  lowering  air  temperature  and  conserving 
litter  moisture.   Insolation,  air  temperature  and  dead  fuel 
moisture  are  factors  in  fire  potential  (Andrews  and  Chases  1989, 
Rothermel  1983,  Andrews  1986) .   Tree  canopies  also  act  as  a 
windbreaks  to  reduce  wind  speed  (Fons  1940,  Marston  1956), 
another  important  factor  in  fire  potential  (Schroeder  and  Buck 
1970)  . 
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Figure  1-7 
Area  1 


Satellite  image  of  1989  fire  in  Water  Conservation 


Transpiration  on  tree  islands  is  also  higher  than  in 
sawgrass  or  wet  prairies  (Gunderson  1992).   Tree  roots  reach 


lower  into  the  water  table  than  the  roots  of  grasses,  sedges, 
and  herbs,  providing  the  trees  with  moisture  during  drought 
after  shallow-rooted  marsh  vegetation  has  begun  to  desiccate. 
Higher  transpiration  on  the  islands  elevates  relative  humidity 
under  the  shaded  canopy  relative  to  the  open  marsh,  further 
conserving  the  moisture  of  the  understory  and  litter  on  tree 
islands.   While  the  understory  at  the  perimeters  of  tree  islands 
consists  of  dense  sawgrass,  vines  and  brush,  providing  ample 
fire  fuel,  the  understory  in  the  interior  of  large  islands  is 
sparse,  consisting  mainly  of  ferns  (Pope  1989,  Wade  et  al. 
1980)  . 


Vegetation 

Vegetation  community  types 

Plant  species  in  the  Everglades  can  be  organized  into 
associations,  or  communities,  which  share  similar  environmental 
requirements.   An  early  classification  of  Everglades  species 
into  community  types  was  based  upon  groups  of  species  observed 
growing  together,  but  did  not  use  hydrologic  regime  as  an 
organizing  criterion  (Davis  1943).   Later,  Loveless  (1959) 
characterized  a  set  of  dominant  communities  in  the  Everglades, 
tree  islands,  sawgrass  marshes,  wet  prairies  and  aquatic 
sloughs,  and  discussed  the  influence  of  hydrologic  regime  in 
structuring  communities. 
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More  detailed  classifications  of  Everglades  plant 
communities  were  made  later  (Craighead  1971,  Wade  et  al.  1980, 
Olmsted  and  Loope  1984).   Gunderson  (1994:327)  divided 
Everglades  freshwater  wetlands  into  forested  communities, 
graminoid  associations  and  areas  with  "little  or  no  emergent 
vegetation."   Forested  communities  include  bayheads,  willow 
heads  and  cypress  forests.   Graminoid  associations  include 
sawgrass  marshes,  peat  wet  prairies  and  marl  wet  prairies. 
Little  or  no  emergent  vegetation  includes  ponds,  creeks  and 
sloughs.   Using  this  community  classification,  the  common  plant 
species  of  Water  Conservation  Area  1  are  listed  in  Table  1-2 . 


Table  1-2.   Common  plant  species  of  Water  Conservation  Area  1, 
by  vegetation  community 


Vegetation 
Community 

Species 

Common  Name 

Bayhead 
(Tree  Island) 

Persea   borbonia 

Red  Bay 

Ilex    cassine 

Dahoon  Holly 

Blechnum   serrulatum 

Swamp  Fern 

Acrostichum   danaei folium 

Leather  Fern 

Osmunda    regalis 

Royal  Fern 

Lygodium   japonicum* 

Japanese  Climbing 
Fern* 

Smilax    lauri folia 

Bamboo  Vine 

Myrica    cerifera 

Wax  Myrtle 

Chrysobalanus    icaco 

Coco  Plum 

Cephalanthus 
occidentalis 

Buttonbush 

Willow  Head 

Salix    caroliniana 

Willow 

Ludwigia    octovalis 

Water  Primrose 

Mikania    scandens 

Climbing  Hemp-vine 

Polygonum 
hydropiperoi des 

Smartweed 

Phragmites    communis 

Common  (Giant)  Reed 

Sawgrass  Marsh 

Cladium   jamaicense 

Sawgrass 

*  indicates  exotic  species 

Sources :  Richardson  et  al.  1990,  Silveira 
data) ,  Thompson  1970,  Alexander  and  Crook 
al.  1994 


1988  (unpublished 
1975,  and  Browder  et 


40 


Table  l-2--continued 


Vegetation 
Community 

Species 

Common  Name 

Peat  Wet 
Prairie 

Rhynchospora    traceyi 

Tracy's  Beakrush 

Eleocharis   montevidensis 

Spikerush 

Eleocharis    quadrangulata 

Spikerush 

Panicum   hemitomon 

Maidencane 

Eriocaulon    compressum 

Hat-pins 

Fuirena    scirpoides 

Rush  Fuirena 

Xyri s    spp . 

Yellow-eyed  Grass 

Hyperi c um    spp . 

St.  John's  Wort 

Micro-algae : 
Myxophyceae, 
Bacillophyceae  and 
Chlorophyceae  groups 

Periphyton  Mats 

Pontederia    cordata 

Pickerelweed 

Sagittaria    lati folia 

Arrowhead 

Slough 

Utricularia    spp. 

Bladderworts 

Potamogeton    illinoensis? 

Pondweed 

Najas    guadalupensis 

Southern  Naiad 

Nuphar   advena (luteum?) 

Spatterdock 

Nymphoides    aquaticum? 

Floating  Heart 

Nymphaea    odorata 

White  Water  Lily 

Eichhornia    crassipes* 

Water  Hyacinth* 

Hydrocotyle    umbellata 

Water  Pennywort 

Salvinia    spp. 

Water  Fern 

Invasive 
Species 

Typha    domingensis 

Cattail 

Melaleuca    quinquenervia* 

Punk  Tree,  Cajeput* 

*  indicates  exotic  species 

Sources :  Richardson  et  al.  1990,  Silveira 
data),  Thompson  1970,  Alexander  and  Crook 
al.  1994 


1988  (unpublished 
1975,  and  Browder  et 


Environmental  factors  in  vegetation  community  distribution 

The  distribution  of  plant  communities  in  the  Everglades 
is  a  function  of  hydrologic  and  nutrient  factors,  along  with 
fire  history,  animal  disturbance,  hurricanes,  and  previous 
vegetation  type  (Davis  1943,  Loveless  1959,  Craighead  1971, 
McPherson  1973,  Davis  1994). 
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Hydrologic  factors 

Plant  species  in  Water  Conservation  Area  1  have  varying 
tolerances  for  depth,  duration  and  variability  of  flooding  (Pope 
1989,  Richardson  et  al .  1990).   The  spatial  distribution  of 
Everglades  vegetative  communities  is  largely  a  function  of 
hydroperiod  (McPherson  1973).   In  semi-permanently  flooded 
areas,  hydroperiod  is  largely  a  function  of  elevation.   In  the 
Everglades  "very  minor  changes  in  elevation  can  cause  major 
shifts  in  the  location  of  plant  communities"  (Alexander  and 
Crook  1974  :  63)  . 

Submerged  and  floating  aquatic  plants  found  in  the  slough 
community  tolerate  continuous  and  deep  flooding.   The  species  of 
the  wet  prairie  require  shorter  hydroperiods  and  shallower  water 
than  slough,  but  tolerate  longer  hydroperiods  than  sawgrass. 
Tree  islands  are  higher  in  elevation  than  the  surrounding 
communities,  and  have  the  shortest  hydroperiod. 

Hydroperiod  at  a  site  is  determined  by  rainfall, 
evapotranspiration,  the  elevation  of  the  peat  surface, 
surrounding  topography,  and  water  flow  patterns.   While  these 
environmental  parameters  influence  the  vegetation  community,  at 
the  same  time  vegetation  also  influences  its  environment. 
Vegetation  increases  elevation  by  adding  dead  biomass  to  the 
peat  and  intercepting  floating  peat  and  litter.   It  reduces 
standing  water  via  transpiration,  impedes  water  flow,  and 
provides  fire  fuel.   Interactions  between  plants  and  the 
environment  drive  vegetation  succession  (Loveless  1959,  Vogl 
1969)  . 
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Nutrient  concentrations 

Before  drainage  and  impoundment,  Water  Conservation  Area 
l's  waters  were  oligotrophic :  nutrients  such  as  nitrogen  and 
phosphorus  entered  the  system  via  rainfall,  and  therefore  were 
limited.   Low  nutrient  levels  gave  oligotrophy-adapted  species 
such  as  sawgrass  a  competitive  advantage  (Steward  and  Ornes 
1975).   Historically,  nutrient  levels  had  little  influence  on 
the  distribution  of  plant  communities,  except  perhaps  in  the 
vicinity  of  faunal  activity  sites  such  as  alligator  holes  and 
bird  rookeries  (Richardson  et  al.  1990).   Cattail  {Typha 
domingensis)    and  willow  {Salix    caroliniana)    grow  at  these  sites 
as  a  result  of  the  nutrient  concentrations  or  physical 
disturbance  caused  by  the  animals.   Since  the  construction  of 
the  canals  on  the  perimeter  of  Water  Conservation  Area  1  and  the 
resulting  influx  of  nutrient-rich  agricultural  runoff,  cattail 
and  willow  head  species  have  also  become  dominant  near  the 
canals  (Richardson  et  al.  1990)  . 

Vegetation  Succession 

The  first  theory  of  vegetation  succession  was  put  forward 
by  Frederick  Clements,  based  on  his  knowledge  of  forests  in  the 
Great  Lakes  region  (Clements  1898).   After  observing  the  species 
composition  and  structure  of  old  growth  forest  and  the  regrowth 
of  cleared  areas,  Clements  proposed  that  land  in  its  natural 
state  passes  through  a  succession  of  vegetation  community  types 
over  time  (serai  stages),  until  it  reaches  a  climax  stage.   Each 
serai  stage,  as  it  matures,  alters  the  conditions  of  its 
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environment  until  the  environment  is  more  favorable  to  other 
species.   Colonization  and  conversion  to  the  next  serai  stage 
then  follow.   Early  serai  stages  are  composed  of  pioneer  species 
proficient  at  rapid  reproduction  and  resource  exploitation.   As 
time  goes  on,  succeeding  stages  have  more  slower-growing, 
resource-conserving  species.   Climax  species  are  adapted  for 
equilibrium  conditions.   Each  geographic  region  has  its  own 
climax  community. 

Clements'  theory  was  widely  accepted,  but  some  later 
ecologists,  notably  Henry  A.  Gleason,  objected  to  the  rigid 
concept  of  a  fixed  serai  progression  and  climax.   They  proposed 
that  species  composition  at  a  site  was  determined  only  by  the 
summation  of  a  variety  of  stochastic  factors,  such  as  seed 
dispersal,  weather  and  fire  history  (Gleason  1926).   Gleason 
denied  the  existence  of  a  linear  set  of  serai  stages  leading  to 
a  regional  climax.   Disturbance  was  a  peripheral  part  of 
Clements'  theory,  but  as  the  field  of  ecology  grew,  ecosystems 
were  described  in  which  disturbance  played  a  critical  role.   In 
these  systems  a  regional  climax  was  never  achieved.   Terms  like 
"fire  climax"  were  coined  to  keep  these  systems  under  the 
umbrella  of  succession  theory. 

The  concept  of  succession  is  useful  for  understanding  the 
spatial  pattern  of  vegetation  in  Water  Conservation  Area  1 . 
Vegetation  community  change  in  accordance  with  Clements'  theory 
can  be  observed  in  peat  wetland  systems  (Vogl  1969)  .   In  Water 
Conservation  Area  1,  water  lilies  in  sloughs  and  alligator  holes 
build  up  peat,  raising  elevation  and  shortening  hydroperiod 


until  the  grasses  and  rushes  of  the  wet  prairie  can  invade. 
Sawgrass  builds  up  peat  until  the  hydroperiod  is  short  enough  to 
favor  shrubs;  shrubs  alter  the  environment  to  allow  the  invasion 
of  trees  (Gleason  and  Stone  1994,  Loveless  1959).   We  do  observe 
a  loosely  predictable  process  of  vegetation  community  change 
after  disturbance,  and  it  is  convenient  to  organize  these 
changes  into  types  or  serai  stages  (i.e.  prairie,  shrub, 
forest) .   Based  on  the  peat  record  and  plant-environment 
relationships,  a  successional  sequence  from  slough  to  wet 
prairie,  to  sawgrass,  to  shrubs  and  then  trees  has  been 
documented  in  the  Everglades  (see  Figure  1-8)  (Loveless  1959, 
Pesnell  and  Brown  1977,  Olmsted  and  Loope  1984,  Duever  1984, 
Kushlan  1990) . 
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Figure  1-8.   Successional  sequence  in  Water  Conservation  Area  1 
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However,  Clements'  emphasis  on  a  linear  progression 
through  the  serai  sequence,  ending  in  equilibrium  at  a  stable 
climax  is  not  applicable  to  the  Everglades  ecosystem. 
Disturbances  and  other  factors  perceived  as  stochastic 
(droughts,  floods,  fires,  alligator  activity,  seed  dispersal  and 
nutrient  concentration  by  birds)  all  influence  plant  community 
distribution  in  the  Everglades.   For  a  young,  disturbance- 
adapted  ecosystem  like  the  Everglades,  the  concept  of  a  stable 
climax  at  equilibrium  is  not  appropriate.   Tropical  hardwood 
hammock  forest  has  been  identified  as  the  climax  community  of 
the  Everglades  (Olmsted  and  Loope  1984),  but  it  is  not  found 
anywhere  in  Water  Conservation  Area  1 . 

A  more  useful  non-linear  concept  of  succession  and 
disturbance  is  found  in  C.S.  Holling' s  cyclical  theory  of 
ecosystem  dynamics  (Holling  1992a) .   In  this  theory,  disturbance 
is  part  of  a  continuum  with  succession.   An  ecosystem  progresses 
from  an  exploitation  phase  (corresponding  to  early  serai  stages) 
to  a  conservation  phase  (late  serai  stages),  as  in  the 
succession  concept.   However,  while  the  system  becomes  more 
stable  and  organized  as  it  approaches  climax,  it  also  becomes 
less  resilient.   Although  more  able  to  absorb  small 
perturbations,  it  is  more  vulnerable  to  a  major  disturbance. 
Inevitably  disturbance  releases  the  stored  energy  and  material 
in  the  system  and  disorganizes  it.   The  system  goes  through  a 
release  phase  and  then  a  reorganization  phase,  from  which  it  may 
follow  the  same  exploitation  and  conservation  series  it  followed 
before.   If  conditions  have  changed,  however,  it  may  follow  a 
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different  path.   This  last  concept  is  especially  applicable  to 
the  Everglades  ecosystem.   As  humans  change  the  hydrology  and 
introduce  exotic  species  to  the  Everglades,  vegetation 
communities  have  a  greater  potential  to  change  in  ways  not 
predicted  by  the  our  current  understanding  of  succession  in  the 
Everglades  landscape. 


CHAPTER  2 
LANDSCAPE  DATA  COLLECTION 


Introduction 


The  foundation  of  an  understanding  of  landscape  pattern  is 
a  thorough  familiarity  with  the  landscape.   Many  sources  of 
information  were  used  in  this  work  to  gain  knowledge  of  the 
landscape  of  Water  Conservation  Area  1.   This  chapter  will 
describe  the  acquisition  and  processing  of  this  information  (see 
Figure  2-1).   Data  sources  included  field  observations  made  by 
airboat  during  1988-1991  and  aerial  photographs  taken  in  1992. 
Aerial  photographs  from  1952  were  compared  with  satellite  images 
from  1987  and  1990  to  assess  recent  landscape  changes. 
Satellite  images  were  studied  in  both  their  raw  form  and  as 
vegetation  maps.   Historical  accounts  of  the  vegetation  of  the 
area  were  used  as  clues  to  the  landscape  in  the  past  (Davis 
1943,  Thompson  1972).   Literature  on  the  natural  history  of  the 
Everglades,  other  Florida  wetlands  and  the  Okefenokee  swamp  also 
was  used  to  provide  additional  information  on  the  species  found 
in  the  area  (Dressier  et  al.  1987,  Craighead  1971,  Glasser  1986, 
Kushlan  1990,  Olmsted  and  Loope  1984,  Pesnell  and  Brown  1977, 
Wade  et  al.  1980) . 

Using  the  knowledge  gained  from  the  above  information,  a 
map  of  Water  Conservation  Area  1  vegetation  was  made  from  a  1987 
satellite  image.   The  image  processing  required  to  create  this 
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map  is  described  in  the  latter  part  of  this  chapter.   Vegetation 
was  computer-classified  from  the  satellite  data  using  field 
observations,  aerial  photographs  and  a  vegetation  classification 
of  the  same  image  by  J.R.  Richardson  (Richardson  et  al.  1990) . 
Vegetation  classes  were  based  on  a  successional  sequence  and 
will  be  referred  to  as  "successional  classes"  to  distinguish 
them  from  the  classes  in  the  Richardson  et  al.  classification. 
The  creation  of  this  successional  vegetation  classification 
required  a  synthesis  of  all  the  landscape  and  community  ecology 
information  gathered  in  this  work. 


1952  Aerial 
Photographs 


1990  Multispectral 
&  1990  Panchromatic 
Satellite  Images 


1989  Raw 

Panchromatic 

Image 


1990 

Merged 

Image 


1987  Multispectral 
&  1987  Panchromatic 
Satellite  Images 


Field  Data  and 
Ecological  Literature 


18-Class  Vegetation 
Map  (J.R.  Richardson) 


1992  Aerial 
Photographs 


|  24-Class  Vegetation 
Wj  Succession  Map 


=  Data  Source 

=  Classified  Image 


Figure  2-1.   Diagram  of  landscape  data  sources 
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Ground  Observations 

Because  travel  on  the  deep  peat  of  Water  Conservation  Area 
1  is  difficult,  few  vegetation  studies  of  the  area  exist,  and 
most  are  based  on  inadequately  sized  samples  (Richardson  et  al. 
1990  and  Pope  1989  are  exceptions) .   To  gain  a  familiarity  with 
the  landscape  as  a  whole  and  with  the  vegetation  at  a  finer 
scale  than  remote  sensing  allows,  observations  of  the  vegetation 
communities  of  Water  Conservation  Area  1  were  made  in  the  field. 
Georef erenced  observations  also  were  made  to  be  used  in 
classifying  satellite  imagery. 

Drought  conditions  and  a  major  fire  interrupted  and 
delayed  these  observations,  but  these  events  also  provided  the 
opportunity  to  observe  fire  effects  and  vegetation  recovery. 

Pre-fire  Field  Observations 

Methods 

Between  1988  and  1991,  Water  Conservation  Area  1  was 
explored  on  the  ground  by  airboat.   A  Loran-C  navigational  unit 
and  paper  prints  of  an  April  1987  satellite  image  were  used  to 
track  the  path  of  the  airboat.   The  Loran-C  unit  used  signals 
from  microwave  stations  to  calculate  latitude  and  longitude 
coordinates  within  an  error  range  of  100  m.   The  Loran-C 
coordinates  were  recorded  at  different  sites,  along  with  the 
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vegetation  types  and  estimates  of  the  coverage  of  those  types 
within  a  50  feet  radius  of  the  airboat. 

Ground  exploration  began  in  the  spring  of  1988.   At  that 
time,  the  water  level  in  Water  Conservation  Area  1  was  high. 
Vegetation  communities  and  corresponding  Loran  coordinates  were 
documented  at  numerous  locations,  and  located  on  the  1987 
satellite  image.   During  1988  however,  drought  conditions 
lowered  the  water  level  in  the  area.   Once  the  water  level 
dropped  below  the  peat  surface,  airboat  travel  was  limited  to 
the  perimeter  canals  and  areas  near  them.   A  high  concentration 
of  alligators  seeking  water  in  the  canals  prompted  the  wildlife 
refuge  managers  to  close  Water  Conservation  Area  1  to  visitors. 
Ground  data  collection  became  impossible. 

Drought  persisted  through  the  winter  of  1988-89,  and  a 
freeze  further  heightened  the  fire  potential  by  increasing  the 
dead  fuel  load.   As  described  in  the  introductory  chapter,  on 
April  11,  1989,  at  the  end  of  the  dry  season,  firebrands  from  a 
controlled  burn  just  outside  the  northwestern  boundary  of  Water 
Conservation  Area  1  blew  across  the  canal  and  started  a  large 
wildfire.   Over  18,000  ha  (40,000  acres)  burned  in  the  northern 
and  western  portions  of  the  area  (see  Figure  2-2) .   The  fire  was 
controlled  by  firefighters  at  the  very  northern  tip  of  Water 
Conservation  Area  1  to  prevent  its  spread  to  private  property, 
but  was  allowed  to  burn  freely  elsewhere  (M.  Maffei,  personal 
communication) .   Because  access  to  the  area  was  still 
impossible,  ground  observation  of  the  immediate  effects  of  the 
fire  could  not  be  made. 
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Figure  2-2.   Map  of  the  1989  fire,  with  the  location  of  the 
successional  vegetation  map  used  in  the  model  in  Chapter  4 
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Post-fire  Field  Observations 


Methods 


Although  the  initial  effects  of  the  fire  could  not  be 
observed,  vegetation  recovery  could  be  observed.   The  drought 
did  not  end  until  1991.   However,  in  November  1989  the  water 
level  temporarily  rose  enough  in  Water  Conservation  Area  1  to 
allow  limited  airboat  travel.   The  status  of  vegetation  recovery 
after  the  fire  was  observed  from  an  airboat  along  a  wide  east- 
west  corridor  running  through  the  center  of  Water  Conservation 
Area  1  (varying  between  26°  30'  and  26°  32'  latitude).   Even 
eight  months  after  the  fire,  the  boundaries  and  effects  of  the 
fire  were  obvious.   Recovery  was  limited,  because  there  had  not 
been  a  drought-free  growing  season  since  the  fire. 

In  February  1991,  after  vegetation  recovery  was  well 
underway,  a  larger  percentage  of  the  burned  area  was  visited  by 
airboat.   Three  burned  and  three  unburned  tree  islands  were 
explored  on  foot,  measured  and  described  in  detail.   A  plan  to 
describe  18  burned  and  18  unburned  islands  was  not  feasible  and 
had  to  be  abandoned.   Brushy  resprouting  trees,  fallen  dead 
trees,  and  dense  growth  of  the  exotic  Japanese  climbing  fern 
(Lygodium   japonicum)    that  had  invaded  burned  islands  made 
traversing  burned  islands  by  foot  extremely  difficult  and  time- 
consuming  . 

The  six  islands  described  were  selected  at  random  using 
pairs  of  latitude/longitude  coordinates  chosen  at  random  within 
the  burned  area,  and  randomly  chosen  cardinal  directions. 
Coordinates  and  cardinal  directions  were  assigned  before 
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sampling  began.   A  sampling  team  traveled  to  a  pair  of  ground 
coordinates  by  airboat,  then  from  there  traveled  in  the  assigned 
cardinal  direction  until  an  island  of  the  desired  type  was 
reached.   Two  transects  were  set  up  on  the  island,  one  running 
south  to  north,  the  other  west  to  east.   The  samplers  followed 
these  transects  on  foot,  recording  whether  or  not  the  vegetation 
showed  signs  of  burning.   The  distance  along  the  transect  was 
recorded  where  the  vegetation  changed  from  burned  to  unburned 
and  vice  versa.   The  numbers  of  tree  individuals  touching  the 
transect  were  recorded  and  tallied  by:  burned  and  living,  burned 
and  dead,  unburned  and  living,  unburned  and  dead.  It  was  assumed 
that  any  burned  trees  that  had  not  resprouted  (22  months  post- 
fire)  were  dead. 
Results 

While  the  information  gathered  on  the  ground  after  the 
fire  was  limited,  it  was  nonetheless  essential  as  ground 
training  for  the  interpretation  of  aerial  photographs  taken 
after  the  fire.   Depictions  of  selected  surveyed  islands  are 
shown  in  Figure  2-3. 

On  the  small  islands  visited  in  the  unburned  area  of  Water 
Conservation  Area  1,  a  central  core  of  trees  with  an  understory 
of  ferns  was  surrounded  by  a  border  of  shrubs,  sawgrass  and 
vines  3  to  4  m  wide.   On  the  small  islands  visited  in  the  burned 
area,  all  vegetation  had  been  burned.   However,  all  trees  along 
the  survey  transects  were  sprouting  from  their  bases.   A  ring  of 
sawgrass  less  than  1  m  wide  surrounded  the  islands. 
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On  the  large  burned  island  surveyed,  the  core  of  the 
island,  made  up  of  trees  with  a  fern  understory,  was  not  touched 
by  fire.   The  margins  of  this  tree  core  were  burned.   The  width 
of  this  burned  zone  varied  from  7  m  to  77  m.   The  wider  width 
was  found  where  small  peninsulas  of  trees  extended  out  from  the 
general  ovoid  shape  of  the  island  (see  Figure  2-3) .   Seventeen 
percent  of  the  burned  trees  along  the  transects  were  dead,  and 
the  rest  were  living.   It  was  not  possible  to  tell  whether  the 
dead  trees  had  been  killed  by  the  fire.   Although  wax  myrtle 
{Myrica    cerifera)    is  not  typically  killed  by  fire,  the  wax 
myrtle  shrubs  making  up  the  southern  "tail"  of  the  island  were 
burned,  with  complete  mortality.   An  extensive  band  of  sawgrass 
around  the  island  had  been  burned,  but  had  recovered  with  no 
sign  of  mortality. 

Aerial  Photographs 

Aerial  photographs  were  used  as  a  link  between  the  fine 
scale  of  ground  observations  and  the  course  resolution  of 
satellite  images.   When  laid  out  in  order  with  their  edges 
matching,  the  photographs  provided  a  landscape  view,  while  still 
showing  details  lost  to  satellite  images  such  as  the  texture  of 
vegetation  and  alligator  trails.   Photographs  were  used  to 
provide  visual  information  on  the  landscape,  to  aid  in  the 
identification  of  vegetation  at  different  sites  for  satellite 
classification,  and  to  detect  vegetation  change.   A  set  of  1952 
photographs  provided  an  approximation  of  vegetation  before 
impoundment  and  was  compared  with  recent  data  to  detect  change. 


55 


540 
m 


36  m 


25  m 


21  m 


Small  burned 
island 


Small  unburned 
island 

Approx.  scale:  1:1700 


N 


Approx.  direction  of 
fire  movement 


Trees,  not  burned 
Trees,  burned 

Shrubs,  burned 

Sawgrass:  burned  on  large  island 


240  m 
Large  burned  island 

Approx.  scale:  1:8600 


unburned,  combined  with  shrubs,  on  small  island 


Figure  2-3.   Map  of  three  islands  observed  on  the  ground  in 
February,  1991 


Photographs  taken  in  1991  and  1992  were  taken  to  observe 
disturbance  and  recovery  from  fire.   The  disadvantage  of  the 
aerial  photographs  used  was  that  they  could  not  be 
georef erenced .   Due  to  the  lack  of  man-made  features  in  Water 
Conservation  Area  1  and  the  repetitive  nature  of  the  landscape 
pattern,  the  matching  of  the  photographs  to  map  coordinates  or 


other  data  could  not  be  verified.   While  an  approximate  scale  of 
the  1952  photography  was  known,  the  actual  scale  of  each 
photograph  could  not  be  calculated. 

Aerial  Photographs,  1952 

A  set  of  black  and  white  aerial  photographs  of  Water 
Conservation  Area  1  was  taken  in  1952.   The  aerial  photograph 
set  contained  29  photographs  at  a  scale  of  approximately 
1:16,500,  covering  approximately  90%  of  Water  Conservation  Area 
1 .   The  L-7  canal  on  the  northwest  border  of  Water  Conservation 
Area  1  was  still  under  construction  at  the  time  the  photographs 
were  taken,  but  the  remaining  perimeter  canals  were  complete. 
Methods 

The  1952  photographs  were  examined  and  visually  compared 
with  a  satellite  image  dated  April  5,  1990  to  detect  vegetation 
change.   Obvious  change  between  1952  and  1990  was  identified. 
All  islands  larger  than  1000  m2  in  the  photographs  were  searched 
for  on  the  satellite  image  to  detect  island  degradation  and 
disappearance . 

Twenty-five  of  the  1952  aerial  photographs  were  of 
sufficient  quality  to  sample  for  tree  islands.   A  quadrat  1  km2 
was  randomly  placed  on  each  photo,  using  a  random  number  table 
to  assign  x  and  y  distances  from  the  upper  left  corner  of  the 
photo.   Within  the  quadrat,  tree  islands  were  delineated  and 
assigned  a  category:  small  (less  than  160  m2 )  ,  medium  to  large 
(between  160  and  2500  m2  )  and  very  large  (greater  than  2500  m2  )  . 
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Results 

The  results  of  the  sample  of  tree  islands  in  the  1952 
photographs  are  listed  in  Table  2-1. 

Table  2-1.   Results  of  tree  island  sample,  1952  photography 


Island  density 
(islands/km2) 

Small  islands 
(less  than  160  m  ) 

93.81 

Medium  to  large  islands 
(between  160  and  2500  m  ) 

5.  62 

Very  large  islands 
(greater  than  2500  m  > 

0.48 

All  islands 

99.  91 

A  careful  comparison  between  the  1952  photographs  and  the 
1990  satellite  image  showed  that  all  tree  islands  greater  than 
1000  m^  in  the  1952  photographs,  could  be  identified  by  shape  in 
the  1990  image.   In  some  cases  the  island's  trees  had  been 
replaced  by  other  vegetation,  or  woody  vegetation  surrounded  the 
island,  but  the  outline  of  the  island  was  still  visible. 

Other  vegetation  communities  were  difficult  to  identify  in 
the  photographs,  due  to  the  lack  of  color  and  the  variation  in 
exposure  between  photographs.   In  particular,  sawgrass  {Cladium 
jamaicense)     could  not  be  distinguished  from  cattail  {Typha 
domingensis)  ,  so  the  extent  of  cattail  in  1952  could  not  be 
estimated.   (Cattail  was  found  to  increase  in  photoplots  and 
vegetation  transects  near  canals  in  Water  Conservation  Area  1 
between  the  1960's  and  1987  (Richardson  et  al.  1990,  Davis  et 
al.  1994) ) . 
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Two  types  of  vegetation  change  were  apparent  however.   The 
first  of  these  was  f looding-induced  change.   In  the  1952 
photographs,  water  was  impounded  north  of  the  Hillsborough 
canal,  along  the  southern  border  of  the  future  water 
conservation  area.   This  water  created  a  linear  pond  devoid  of 
emergent  vegetation.   A  network  of  wide  sloughs  extended  north 
from  the  pond  approximately  1  km.   The  sloughs  cut  through 
vegetation  community  patches,  and  were  clearly  a  recent 
phenomenon  altering  the  landscape  pattern.   Tree  island  canopies 
in  this  aquatic  zone  were  difficult  to  delineate.   This 
phenomenon  was  documented  later,  in  1970  (Thompson  1970) .   The 
satellite  image  from  1990  showed  a  similar  aquatic  zone,  not 
only  in  the  south,  but  along  all  the  canals  surrounding  Water 
Conservation  Area  1  (shown  in  Figure  2-4  as  zones  of  change) . 

The  second  type  of  change  was  related  to  reduced  flooding. 
The  1990  satellite  image  shows  dense  concentrations  of  shrubs  at 
the  northern  tip  and  in  the  south  central  portion  of  Water 
Conservation  Area  1  (also  part  of  the  zones  of  change  shown  in 
Figure  2-4) .   The  dominant  shrubs  in  these  areas  are  wax  myrtle 
and  buttonbush  {Cephalanthus    occidentalis)  .       The  same  shrub 
areas  appear  in  a  1987  satellite  image.   However,  neither  of 
these  shrub  communities  are  visible  in  the  1952  photographs. 
They  occurred  on  topographic  high  points  (Richardson  et  al. 
1990),  and  may  be  the  products  of  fire  frequency  change  or 
shorter  hydroperiods  since  the  loss  of  sheet  flow. 
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Figure  2-4.   Areas  of  vegetation  change  in  Water  Conservation 
Area  1  between  1952  and  1990,  and  locations  of  photoplots  from 
two  vegetation  studies 
Sources:  Richardson  et  al.  1990  and  Davis  et  al.  1992 


Discussion 

The  1952  photographs  were  analyzed  visually  for  change,  in 
a  general,  non-quantitative  fashion.   A  conscious  decision  was 
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made  to  avoid  standard  photogrammetric  change  detection  methods, 
based  on  several  factors.   First,  the  precise  matching  or 
registration  of  photographs  of  the  same  site  from  different 
years  is  critical  for  change  detection.   The  interior  of  Water 
Conservation  Area  1  is  devoid  of  man-made  feature  or  permanent 
geographic  landmarks  such  as  mountain  peaks .   The  vegetation, 
including  tree  islands,  is  in  flux.   Therefore,  accurate 
registration  of  photography  from  different  years  is  not 
possible.   The  problem  is  compounded  by  the  different  scales  of 
available  photography  and  imagery. 

Second,  the  consistent  delineation  of  vegetation  polygons 
from  remotely  sensed  data  is  difficult  in  Water  Conservation 
Area  1.   The  apparent  boundaries  of  vegetation  communities  are 
strongly  influenced  by  water  level.   Patches  of  emergent 
vegetation  appear  to  shrink  when  water  is  high.   Water  level  in 
Water  Conservation  Area  1  fluctuates  constantly,  introducing 
interpretation  differences  between  seasons  or  years.   In 
addition,  the  available  types  of  remotely-sensed  data  from 
different  years  may  vary  from  black  and  white  photographs,  to 
color  infrared  photographs,  to  satellite  imagery.   These 
different  data  types  may  yield  different  interpretations  due  to 
their  different  spectral  qualities. 

Finally,  the  photo-interpretation  of  the  entire  Water 
Conservation  Area  would  be  too  labor-intensive  to  undertake  in 
the  scope  of  this  work.   A  sampling  of  the  landscape  by 
photographic  plots  could  have  been  done  to  reduce  the  labor  of 
interpretation.   Change  detection  photo-sampling  compares 
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detailed  interpretations  of  vegetation  on  sample  plots  using 
photographs  or  satellite  images  over  time  (Silveira  1996). 
Interpretations  of  vegetation  communities  are  usually  recorded 
as  overlay  maps.   Areas  changed  between  dates  on  each  plot  are 
calculated  by  comparing  overlay  maps.   The  results  are  often 
used  to  draw  conclusions  regarding  vegetation  change  in  the 
ecosystem  as  a  whole.   However,  important  landscape  changes  may 
be  missed  when  an  area  is  sampled  rather  than  examined  in  its 
entirety,  particularly  if  the  sample  size  is  small. 

The  oversight  of  landscape  change  was  realized  in  previous 
photoplot  studies  of  Water  Conservation  Area  1  (Thompson  1972, 
Richardson  et  al.  1990,  Alexander  and  Crook  1975,  Davis  et  al. 
1994).   The  studies  failed  to  detect  the  large  area  of  sawgrass 
that  changed  to  brush  in  the  south-central  part  of  Water 
Conservation  Area  1,  because  no  plots  were  located  in  this 
region.   In  addition,  these  studies  found  conflicting  trends  in 
the  wet  prairie  and  sawgrass  communities  in  central  Water 
Conservation  Area  1. 

Thompson  (1972)  estimated  change  in  vegetation  communities 
in  Water  Conservation  Area  1  between  1948  and  1968  using  a  time 
series  of  7  40.5-ha  photoplots .   Richardson  et  al .  (1990) 
compared  these  same  plots  with  a  1987  classified  satellite  image 
(see  Figure  2-4).   Alexander  and  Crook  (1975)  estimated 
vegetation  change  between  1940-53  and  1965-71  on  100  259.2-ha 
photoplots  throughout  the  Everglades,  including  5  plots  in  Water 
Conservation  Area  1.   Davis  et  al.  (1994)  compared  25  of  these 
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plots  with  1987-1987  photographs,  including  all  the  plots  in 
Water  Conservation  Area  1  (see  Figure  2-4). 

Thompson  (1972)  and  Richardson  et  al.  (1990)  found  small 
vegetation  community  changes  in  their  7  photoplots,  except  in  2 
central  plots.   In  plot  4,  25%  of  the  plot  changed  from  sawgrass 
to  wet  prairie  between  the  years  1962  and  1968.   By  1987 
however,  sawgrass  and  wet  prairie  had  returned  to  their  1962 
distribution.   In  plot  3,  a  similar  switch  from  sawgrass  to  wet 
prairie  dominance  occurred  between  1968  and  1987  (see  Table  2- 
2).   Richardson  et  al.  hypothesized  that  the  changes  found  in 
plot  4  between  1962  and  1968  were  due  to  fire.   A  large  fire 
burned  in  Water  Conservation  Area  1  in  1962  (see  Table  1-1) . 
The  later  changes  in  plot  4  may  have  been  post-fire  vegetation 
recovery. 

Davis  et  al.  (1994)  analyzed  a  plot  in  the  center  of  Water 
Conservation  Area  1  (plot  75)  that  overlapped  the  above- 
mentioned  plot  3.   The  changes  in  plot  75  are  shown  in  Table  2- 
2,  and  are  the  opposite  of  trends  found  by  Richardson  et  al.  in 
plot  3. 


Table  2-2.  Percent  area  in  sawgrass  and  wet  prairie  of  two 
photographic  sample  plots  in  Water  Conservation  Area  1 


Year 

Plot  75 
Sawgrass 

Plot  75 
Wet  Prairie 

Plot  3 
Sawgrass 

Plot  3 

Wet  Prairie 

1968 

0.0% 

89.  6% 

58  .0% 

21.0% 

1984/87 

37.7% 

47.3% 

5.7% 

69.8% 

Source:  Richardson  et  al.  1990:134  and  Davis  et  al.  1994:441 


Two  factors  may  account  for  the  discrepancy  between  the 
results  of  the  Richardson  et  al.  and  Davis  et  al.  studies.   One 
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factor  is  the  interpretation  of  the  1968  photography  by  Thompson 
(1972)  and  Alexander  and  Crook  (1975).   Alexander  and  Crook 
found  no  sawgrass  in  plot  75  (see  Table  2-2),  while  Thompson 
interpreted  58%  of  his  plot  as  sawgrass.   Since  plot  75 
contained  plot  3,  it  should  have  held  some  sawgrass  if 
Thompson's  interpretation  of  plot  3  was  correct.   The  water 
level  in  Water  Conservation  Area  1  in  1968  was  high,  and  the 
area  was  flooded  nearly  all  year  (Thompson  1972).   When  the 
lower  portions  of  sawgrass  plants  are  submerged,  sawgrass  and 
wet  prairie  appear  similar  in  aerial  photographs.   The 
difference  in  plot  composition  must  have  been  due  to  the 
judgment  of  the  photo-interpreters. 

The  other  factor  contributing  to  the  conflicting  trends  is 
the  interpretation  of  the  plots  at  the  later  dates.   Some 
differences  between  these  interpretations  are  expected,  because 
Richardson  et  al.  used  satellite  imagery  instead  of  photography, 
the  interpretations  are  separated  by  3  years,  and  rainfall 
varied  between  these  years.   If  a  fire  burned  in  the  area  of  the 
plots  during  this  interval,  then  a  real  change  in  vegetation 
composition  could  have  occurred. 

Based  upon  the  data  from  their  plots  in  Water  Conservation 
Area  1  and  9  other  plots  in  the  Everglades,  Davis  et  al. 
concluded  that  a  significant  increase  in  sawgrass  and  a 
significant  decrease  in  wet  prairie  had  occurred  throughout  the 
Everglades'  "wet  prairie/slough-sawgrass-tree  island  mosaic," 
(Davis  et  al.  1994:432).   However,  the  different  results  of 
photoplot  studies  in  Water  Conservation  Area  1  illustrate  the 
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limitations  of  using  a  small  sample  of  photographic 
interpretations  to  characterize  vegetation  community  change 
within  an  Everglades  landscape. 

Aerial  Photographs,  1991-1992 

A  satellite  image  was  used  to  assess  the  effects  of  the 
1989  fire  on  the  landscape  of  Water  Conservation  Area  1.   The 
first  post-fire  cloud-free  image  available  for  this  purpose  was 
dated  April  1990.   To  assess  the  progress  of  vegetation  recovery 
after  1990,  a  low-cost  landscape  overview  was  needed. 
Photography  was  preferred  to  distinguish  resprouting  trees  from 
shrubs.   It  was  not  necessary  for  the  photographs  to  be 
vertical;  in  fact,  oblique  angles  were  useful  for  vegetation 
identification.   Because  professional  aerial  photography  was  too 
costly,  a  photographic  survey  was  made  using  a  35  mm.  camera 
held  out  the  window  of  a  small  plane. 
Methods 

On  December  1991,  20  months  after  the  1989  fire,  an  aerial 
photographic  survey  of  Water  Conservation  Area  1  was  conducted. 
Photographs  were  taken  from  an  Cessna  172  airplane,  using  a 
hand-held  35  mm  camera.   Both  near-vertical  and  oblique  color 
photographs  were  taken.   The  survey  flight  was  part  of  a  wading 
bird  survey  (H.  Jelks,  unpublished  data),  but  was  also  well 
suited  for  a  landscape  inventory.   East-west  transects  were 
flown  on  two  consecutive  mornings,  at  a  constant  speed  and  305  m 
altitude.   The  transects  followed  known  coordinates,  spaced  0.5 
minutes  latitude  apart.   Photographs  were  taken  from  the 
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aircraft  at  regular  intervals,  Loran  coordinates  were  noted  for 
distinctive  landmark  features,  and  additional  notes  were  taken 
between  transects.   Weather  conditions  were  clear. 

On  March  24,  1992,  a  return  flight  was  made  to  document 
the  continued  recovery  of  tree  islands  after  the  fire.   The  time 
of  day,  aircraft,  and  method  of  photographing  were  the  same  as 
the  first  survey.   Rather  than  following  a  systematic  flight 
path  however,  individual  islands  were  revisited  and  photographed 
at  altitudes  ranging  from  300  to  750  m. 
Results 

Three  hundred  and  eighty-eight  photographs  were  taken 
during  the  1991  survey.   The  quality  of  the  photographs  was 
good,  and  the  lighting  and  weather  conditions  were  excellent  for 
recording  the  texture  of  the  vegetation.   In  addition  to  showing 
the  condition  of  the  vegetation,  the  photographs  were  useful  in 
identifying  vegetation  types  for  the  classification  of  a  1987 
satellite  image,  as  described  in  the  next  section. 

During  the  second  survey  60  photographs  were  taken.   The 
weather  was  hazy,  so  the  color  quality  of  the  photographs  was 
poor.   The  purpose  of  the  survey  was  to  view  the  extent  of 
resprouting  by  the  tree  stumps  and  snags  left  by  the  fire.   The 
lack  of  color  definition,  combined  with  the  dense  growth  of  the 
exotic  Japanese  climbing  fern  {Lygodium   japonicum)    over  tree 
snags  prevented  this. 
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Satellite  Images 

Five  satellite  images  of  Water  Conservation  Area  1  were 
used  for  a  variety  of  purposes  in  this  work.   The  images  were 
used  in  their  raw  form  for  visual  analysis  in  the  same  way 
aerial  photographs  might  be  used:  as  navigation  aids,  to 
identify  vegetation,  and  to  identify  the  extent  and  behavior  of 
the  1989  fire.   Images  also  were  classified  into  two  vegetation 
maps.   The  maps  were  used  in  a  geographic  information  system 
(GIS)  to  analyze  tree  island  distributions,  (see  Chapter  3),  and 
to  model  landscape  dynamics,  (see  Chapter  4) .   The  satellite 
data  were  processed  on  the  computer  to  make  the  most  of  their 
spatial  and  spectral  resolution  and  to  create  the  vegetation 
maps.   This  image  processing  is  described  in  detail  below. 

SPOT  Satellite  Image  Information 

The  satellite  data  were  acquired  from  the  French  SPOT 
Image  Corporation.   SPOT  Image  Corporation  currently  has  two 
satellites  deployed  in  sun-synchronous  orbits  823  km  above  the 
earth.   The  swath  width  covered  by  the  satellite  sensors  is  60 
km,  so  each  satellite  image  covers  a  parallelogram  60  km  on  a 
side  . 

The  SPOT  satellites  have  two  sensors:  one  panchromatic  and 
one  multispectral .   The  panchromatic  sensor  measures  reflectance 
of  visible  light  in  the  green  to  red  portion  of  the  spectrum 
(.51-. 73  micrometers).   It  has  a  10  m  pixel  resolution,  so  it 
senses  and  records  one  panchromatic  reflectance  value  for  each 
10-m  by  10-m  square  of  earth  it  scans.      The  multispectral 
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sensor  has  a  20  m  pixel  resolution,  and  collects  three  data 
values  for  each  pixel.   The  three  data  values  are  recorded  in 
separate  files  or  bands .   Band  1  is  the  reflectance  of  green 
visible  light  (.50-. 59  micrometers).   Band  2  is  the  reflectance 
of  red  visible  light  (.61-. 68  micrometers) .   Band  3  is  the 
reflectance  of  near  infrared  radiation  (.79-. 89  micrometers). 
The  images  used  in  this  work  are  listed  in  Table  2-3. 

Table  2-3.   SPOT  satellite  images  of  Water  Conservation  Area  1 


Scene  Date 

Image  Type 

SPOT  Scene  ID  Number 

April  4,  1987 

multispectral 

12X623297870404155912 

April  4,  1987 

panchromatic 

11P623297870404155910 

April  11,  1989 

panchromatic 

11P62 32 97 890411155 903 

April  5,  1990 

multispectral 

11X62  32  97  900  4  051603  35 

April  5,  1990 

panchromatic 

12P62 32 97 900405160333 

A  SPOT  image  dating  from  April  11,  1989  was  taken  during 
the  severe  fire  that  burned  over  18,000  ha  of  Water  Conservation 
Area  1.   Smoke  obscured  part  of  Water  Conservation  Area  1,  and 
created  spectral  variation  in  the  image  data,  preventing 
vegetation  interpretation  or  classification  of  the  multispectral 
image.   The  panchromatic  image  was  nonetheless  useful  for  visual 
analysis.   It  was  used  as  a  documentation  of  fire  behavior 
during  the  1989  fire,  particularly  as  it  affected  tree  islands. 

The  multispectral  SPOT  image  dating  from  April  4,  1987  was 
merged  with  the  panchromatic  image  from  the  same  date  to  improve 
the  utility  of  the  imagery.   Likewise,  the  April  5,  1990 
multispectral  image  was  merged  with  the  panchromatic  image  from 
April  5,  1990.   The  result  of  merging  the  panchromatic  and 
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multispectral  images  was  an  enhanced  data  set  with  the  3-band 
spectral  resolution  and  10-m  spatial  resolution. 

The  1987  merged  SPOT  image  was  used  in  several  ways.   The 
image  was  rectified  to  State  Plane  Coordinates,  and  classified 
into  an  18-class  vegetation  map  by  J.  R.  Richardson,  as  part  of 
an  evaluation  of  vegetation-water  relationships  in  Water 
Conservation  Area  1  (Richardson  et  al.  1990) .   Paper  printouts 
of  this  vegetation  map  were  used  to  navigate  Water  Conservation 
Area  1  by  airboat . 

Both  the  unclassified  1987  merged  image  and  the  vegetation 
map  were  used  to  gain  familiarity  with  the  appearance  of  the 
landscape  in  the  imagery  medium  and  to  identify  the  spectral 
signatures  of  different  vegetation  communities.   The  sites 
visited  by  airboat  in  Water  Conservation  Area  1  were  located  on 
the  vegetation  map  and  the  unclassified  merged  image  using  the 
Loran-C  coordinates  of  the  sites.   The  spectral  characteristics 
of  the  site  were  compared  with  field  notes  taken  at  the  site. 
The  vegetation  map  was  also  used  to  analyze  the  size 
distribution  of  tree  islands  in  Water  Conservation  Area  1,  using 
a  Geographic  Information  System  (see  Chapter  3) . 

In  a  separate  procedure,  the  multispectral  SPOT  image  from 
April  4,  1987  was  merged  with  the  panchromatic  image  from  the 
same  date,  and  a  subset  of  the  merged  image  was  classified  into 
24  classes  based  on  vegetation  type  and  successional  stage  (see 
Table  2-7).   The  Richardson  et  al.  vegetation  map,  the  1992 
aerial  photographs,  field  notes  and  ecological  information  were 
used  to  define  the  24  classes  (see  Figure  2-1) .   One  copy  of 
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this  successional  class  map  was  rectified  to  State  Plane 
coordinates,  and  a  second  copy  was  rectified  to  Universal 
Transverse  Mercator  (UTM)  coordinates.   This  allowed  the 
successional  class  map  to  be  referenced  to  data  in  both  State 
Plane  and  UTM  coordinates.   The  successional  class  map  was  used 
for  modeling  landscape  change  in  response  to  fire  (see  Chapter 
4),  and  was  compared  with  the  1990  SPOT  images. 

The  1990  SPOT  images  were  merged  and  rectified  to  UTM 
coordinates.   The  1990  merged  image  was  compared  with  1952 
aerial  photographs  as  described  above  in  this  chapter.   It  was 
also  used  to  assess  the  performance  of  the  landscape  model  by 
providing  an  example  of  post-fire  vegetation  recovery  (see 
Chapter  4 ) . 

Image  Merging 

The  image  processing  techniques  used  in  this  work  included 
image  merging,  vegetation  index  calculation,  rectification  and 
classification.   All  image  processing  tasks  and  GIS  analyses  in 
this  work  were  performed  using  ERDAS  Inc.  computer  software, 
versions  7.4  and  7.5  (ERDAS  1991).   The  image  merging  and 
vegetation  index  calculations  were  performed  similarly  for  the 
1987  panchromatic-multispectral  image  pair  and  the  1990 
panchromatic-multispectral  image  pair,  and  so  are  described  only 
once  below. 

Image  merging  was  a  two-step  process.  The  first  step  was 
the  resampling  of  the  multispectral  data  at  10-m  resolution  by 
dividing  the  20-m  pixels  into  quarters.   The  resampled 
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multispectral  data  were  then  registered  to  the  panchromatic 
data.   Once  the  two  images  were  matched  together  with  10-m 
resolution,  the  second  stage  of  the  merging  process  was  an 
incorporation  of  the  panchromatic  data  into  the  multispectral 
data  through  a  color  transformation  procedure.   A  normalized 
difference  vegetation  index,  (NDVI),  was  also  calculated  from 
the  multispectral  image,  and  added  to  the  merged  image  as  a 
fourth  data  band. 
Resampling 

Each  pixel  in  a  SPOT  image  has  a  data  value  for  each  band, 
and  a  row  number  (y  coordinate)  and  column  number  (x  coordinate) 
to  show  its  location  in  space.   Resampling  the  multispectral 
data  reassigned  these  spatial  coordinates,  as  shown  in  Figure  2- 
5.   This  increased  the  size  of  the  data  file  by  a  factor  of  4. 


A  20-meter  pixel 


column  (x)  =  1 
row  (y)  =1 
data  (d)=  5 


The  same  20 -meter  pixel 
resampled  to  10-meter  pixels 


x=l 

x=2 

y=i 

y=i 

d=5 

d=5 

x=l 

x=2 

y=2 

y=2 

d=5 

d=5 

Figure  2-5.   Pixel  coordinate  values  before  and  after  resampling 


The  multispectral  data  were  resampled  at  the  same  time  the 
multispectral  data  were  registered  to  the  panchromatic  data. 
Registration  is  the  process  of  matching  the  x  and  y  coordinates 
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of  the  multispectral  and  panchromatic  images.   Twenty-five 
corresponding  points  were  identified  on  the  two  images  (control 
points) .   These  points  were  man-made  features  that  could  be 
clearly  identified  on  both  images.   For  each  control  point,  x 
and  y  coordinates  on  each  of  the  images  were  recorded  and  paired 
together.   The  pairs  of  coordinates  for  all  the  control  points 
were  run  through  a  computer  algorithm  that  identified  the  mean 
geometric  offset  of  the  points,  and  created  an  algebraic  matrix 
to  transform  the  coordinates  of  the  multispectral  pixels  to  the 
coordinates  of  the  panchromatic. 

The  computer  program  that  generated  the  transformation 
matrix  also  calculated  the  root  mean  square  of  the  distance 
between  the  actual  coordinates  of  the  control  points  and  the 
coordinates  predicted  for  each  point  by  the  transformation 
matrix.   The  transformation  matrix  for  both  the  1987  and  1990 
image  pairs  had  a  root  mean  square  distance  less  than  one  half 
of  a  pixel  (5m). 

While  the  multispectral  pixels  were  spatially  redefined 
during  the  resampling  and  registering  process,  the  multispectral 
data  values  were  not  altered.   "Nearest  neighbor"  resampling  was 
used,  so  the  spectral  data  values  of  the  original  pixels  were 
simply  reassigned  to  the  nearest  spatially  redefined  pixel. 
Image  transformation 

Once  the  two  images  were  matched  together  with  10-m 
resolution,  the  panchromatic  data  were  incorporated  into  the 
multispectral  data  via  a  color  transformation  procedure 
(Richardson  et  al.  1990). 
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By  convention,  the  data  values  in  each  of  the  three  SPOT 
multispectral  image  bands  are  designated  as  color  values  of  red, 
green  and  blue  to  emulate  a  color  infrared  photograph  for 
display  purposes.   Multispectral  values  range  from  0  to  255. 
Band  1  data  are  assigned  as  blue  values,  band  2  as  green,  and 
band  3  as  red.   When  the  red,  green  and  blue  color  values  of 
each  pixel's  three  data  bands  are  displayed  on  a  monitor,  one  of 
16,777,216  possible  colors  appears.   This  color  can  be  plotted 
in  3-dimensional  color  space,  with  axes  for  red,  green  and  blue 
and  values  ranging  from  0  to  255.   In  addition  to  red-green-blue 
(RGB)  color  space,  this  pixel  color  can  be  defined  by  an 
alternate  colorspace-  a  polar  coordinate  system  known  as 
intensity,  hue  and  saturation  (IHS). 

To  merge  the  panchromatic  and  multispectral  data,  the 
multispectral  pixel  values  in  RGB  colorspace  were  transformed  to 
IHS  values.   Next  the  intensity  value  for  each  multispectral 
pixel  was  deleted  and  replaced  by  the  value  found  at  that  pixel 
location  on  the  panchromatic  image.   Finally,  the  IHS 
coordinates  of  the  pixels  were  converted  back  to  RGB  colorspace. 
The  band  3  (red),  band  2  (green),  and  band  1  (blue)  values  had 
been  altered  by  the  inclusion  of  the  panchromatic  information. 
This  process,  using  the  spectral  values  of  a  sample  pixel,  is 
shown  in  Figure  2-6. 

The  advantage  of  a  merged  image  is  the  addition  of  data 
with  higher  spatial  resolution  to  the  multispectral  data.   It  is 
not  possible  to  increase  the  spatial  resolution  of  multispectral 
data  itself.   The  addition  of  panchromatic  data  to  the 
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multispectral  image  greatly  enhances  the  definition  of  small 
vegetation  patches,  ponds,  sloughs,  roads,  and  buildings. 
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Figure  2-6.   Diagram  of  merging  of  images  using  color 
transformation  process 


The  spectral  resolution  of  the  multispectral  data  was  not 
increased  by  the  merging.   The  panchromatic  sensor  records  data 
in  a  single  band  that  lumps  together  reflectance  from  the  green 
and  red  portions  of  the  spectrum  already  sensed  by  bands  1  and  2 
of  the  multispectral  sensor.   Some  spectral  information  was 
actually  lost;  the  intensity  values  from  the  infrared  portion  of 
the  spectrum  (band  3)  were  replaced  by  panchromatic  data  that 
did  not  include  infrared  reflectance.   However,  the  infrared  hue 
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and  saturation  data  which  contain  the  majority  of  information 
describing  infrared  reflectance  were  retained. 

While  all  band  values  decreased  as  a  result  of  the 
transformation  process  (see  Figure  2-7),  they  decreased 
consistently.   The  relationships  between  the  band  values  changed 
little,  as  shown  by  the  band  ratios  in  Figure  2-8. 
Vegetation  index 

A  normalized  difference  vegetation  index,  or  NDVI,  is  a 
ratio  of  the  red  and  infrared  bands  which  highlights  vegetation 
with  high  biomass  and/or  photosynthetic  productivity  (Jensen 
1986) .   An  NDVI  data  band,  (band  4),  was  added  to  the  images  by 
calculating  the  following  formula  for  each  pixel: 

(band  3  -  band  2)/  (band  3  +  band  2+0.5)  *  100 

To  avoid  values  of  0,  0.5  was  added  to  the  denominator,  and  a 
multiplier  of  100  was  used  as  a  scaling  factor  to  put  the  NDVI 
in  the  same  range  as  the  other  bands. 

The  NDVI  band  was  used  in  addition  to  the  3  other  spectral 
bands,  and  its  calculation  did  not  affect  the  data  in  the  other 
bands . 


75 


24 

25 

24 

24 

24 

26 

26 

23 

26 

27 

27 

25 

26 

27 

27 

26 

Panchromatic  values  (band  1)  of  pixels  in  mesic  forest 


66 
17 
30 

66 
17 
30 

64 
18 
30 

64 
18 
30 

66 
17 
30 

66 

17 
30 

64 
18 
30 

64 
18 
30 

67 

17 
31 

67 

17 
31 

67 
17 
30 

67 
17 
30 

67 
17 
31 

67 
17 

31 

67 
17 
30 

67 
17 
30 

Multispectral  red,  green  and  blue  values  (bands  3,  2,  and  1)  of 
pixels  registered  to  same  the  locations 


38 

9 

17 

39 
10 
18 

37 
10 
17 

37 
10 
17 

38 

9 

17 

41 
10 

18 

40 
11 
19 

35 
10 
16 

41 
10 
19 

43 

10 
19 

43 
10 
19 

39 
10 
17 

41 
10 
19 

43 
10 
19 

43 
10 
19 

41 
10 
18 

Multispectral  red,  green  and  blue  values  (bands  3,  2,  and  1) 

after  merging 
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Rectification 

Rectification  of  a  satellite  image  is  the  process  of 
linking  a  map  coordinate  system  to  the  pixels  of  the  image.   The 
1990  merged  image  was  rectified  to  Universal  Transverse  Mercator 
(UTM)  coordinates  without  altering  the  pixel  size.   The 
vegetation  map  made  from  the  1987  merged  image  of  Water 
Conservation  Area  1  (Richardson  et  al .  1990)  was  rectified  to 
State  Plane  coordinates.   During  the  rectification,  the  pixel 
size  was  redefined  from  10  m  to  30  ft.   The  unclassified  1987 
merged  image  was  rectified  to  UTM  coordinates  by  registering  it 
to  the  rectified  1990  merged  image,  using  the  procedure 
described  in  the  Registration  section  above.   A  copy  of  the  1987 
merged  image  was  also  rectified  to  State  Plane  coordinates  using 
the  same  algebraic  transformation  matrix  used  by  Richardson  et 
al.  to  rectify  the  vegetation  map.   This  allowed  the 
unclassified  1987  merged  image  to  be  compared  on  a  pixel-by- 
pixel  basis  with  both  the  1990  merged  image  and  Richardson  et 
al.  1987  vegetation  map. 

To  rectify  the  1990  merged  image,  25  to  30  control  points 
were  used.   These  points  were  man-made  features,  such  as  canals, 
road  intersections  and  bridges,  which  could  be  clearly 
identified  on  both  the  image  and  1:24,000  scale  USGS  topographic 
maps.   The  control  points  were  distributed  as  evenly  as  possible 
throughout  the  image.   There  were  no  suitable  points  in  the 


interior  of  Water  Conservation  Area  1,  but  many  points  were 
located  along  the  perimeter  canals. 

To  identify  the  coordinates  of  control  points,  the 
locations  of  the  image  control  points  were  digitized  on  the 
computer  monitor.   The  locations  of  the  map  control  points  were 
digitized  from  a  digitizing  table.   To  achieve  this,  a  portion 
of  the  image  was  displayed  on  the  monitor  while  a  corresponding 
map  was  "locked  down"  to  the  digitizing  table  by  digitizing 
several  UTM  tic  marks  along  the  edge  of  the  map.   The  accuracy 
of  the  lock-down  was  tested  by  digitizing  points  with  known 
coordinates  and  comparing  the  digitizer's  estimate  of  the 
points'  coordinates  with  their  actual  coordinates.   If  the 
digitized  coordinates  varied  more  than  10  m  from  the  actual 
coordinates,  the  map  lock-down  procedure  was  repeated.   If 
repeated  lock-downs  did  not  reduce  the  discrepancy  to  less  than 
10  m,  the  map  was  assumed  to  be  distorted  and  was  not  used. 

Thirty  control  points  were  digitized  on  the  image  and  the 
map  to  create  a  set  of  corresponding  image  file  coordinates  and 
map  coordinates,  in  the  same  way  that  the  panchromatic  and 
multispectral  images  were  registered  together.   The  pairs  of 
coordinates  for  all  the  control  points  were  run  through  a 
computer  algorithm  which  identified  the  mean  geometric  offset  of 
the  points,  and  created  an  algebraic  matrix  to  rectify  the  image 
to  UTM  coordinates .   The  program  that  generated  the 
transformation  matrix  also  calculated  the  root  mean  sguare  of 
the  distance  between  the  actual  coordinates  of  the  control 
points  and  the  coordinates  predicted  for  each  point  by  the 
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transformation  matrix.   The  error  contributed  by  each  control 
point  was  also  identified.   If  the  root  mean  square  distance  was 
greater  than  1  pixel  (10  m)  ,  the  control  point  contributing  the 
greatest  error  was  discarded,  and  the  transformation  matrix  was 
recalculated.   This  process  was  repeated  until  the  root  mean 
square  distance  was  less  than  1  pixel. 

The  1990  merged  image  was  rectified  using  a  linear 
rectification,  i.e.  the  algebraic  transformation  matrix  was  a 
linear  function.   Using  the  transformation  matrix,  the  image 
file  coordinates  were  converted  to  map  coordinates.   Spectral 
values  were  not  altered,  but  were  assigned  to  the  nearest 
rectified  pixel. 

It  was  unnecessary  to  adjust  for  elevation  in  this 
rectification,  due  to  the  height  of  the  satellite  and  the  flat 
terrain.   The  SPOT  satellite  has  a  fixed  array  of  sensors 
aligned  in  the  x  direction,  and  the  data  are  corrected  for  the 
curvature  of  the  earth.   Therefore,  SPOT  images  do  not  have  the 
displacement  and  distortion  found  in  aerial  photographs. 
Distortion  in  the  image  can  occur  in  the  y  direction,  if  the 
speed  of  the  satellite  varies  while  the  image  is  being  sensed. 

A  rectification  using  control  point  coordinates  identified 
from  a  ground  survey  or  Global  Positioning  System  is  preferable 
to  one  using  control  points  taken  from  a  map.   However,  the 
accuracy  of  a  map-based  rectification  is  appropriate  to  the  10-m 
resolution  of  the  data  when  the  control  points  are  numerous, 
well  distributed  and  checked  for  error. 
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Image  Classification 

During  satellite  image  classifications,  pixels  are 
assigned  to  classes  based  upon  their  spectral  values. 
Classifications  may  be  unsupervised,  supervised,  or  a 
combination  of  both.   In  an  unsupervised  classification,  the 
classes  to  which  pixels  are  assigned  are  derived  from  the 
spectral  attributes  of  the  data  and  are  defined  by  a  clustering 
algorithm.   Only  the  number  of  classes  in  the  classification  are 
specified  beforehand.   The  pixels  are  assigned  to  these  classes, 
and  the  analyst  studies  the  locations  of  the  different  classes 
on  the  image  and  identifies  their  vegetation  or  land  cover  type. 

In  a  supervised  classification,  the  analyst  first  defines 
known  areas  of  vegetation  on  the  image.   Statistics  on  the 
spectral  data  of  the  pixels  in  each  area  are  used  to  define 
spectral  signatures  for  each  class.   An  algorithm  such  as 
maximum  likelihood  is  used  to  assign  the  image  pixels  to  the 
classes  (ERDAS  1991,  Jensen  1986) . 

The  Richardson  et  al .  18-class  1987  vegetation  map  was 
classified  using  a  supervised  procedure  in  which  approximately 
25  classes  were  initially  defined,  and  then  narrowed  down  to  18 
classes  after  testing  and  evaluation  (Richardson  et  al.  1990). 
The  classification  process  gave  special  attention  to  the 
identification  of  cattail  {Typha   domingensis)    and  other 
vegetation  in  the  areas  of  greatest  change . 

To  create  a  map  suitable  for  modeling  the  least-changed 
vegetation  of  Water  Conservation  Area  1,  a  successional  class 
map  was  classified  from  a  subset  of  the  1987  merged  image  using 
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an  unsupervised  procedure.   The  unsupervised  classification 
focused  on  identifying  vegetation  in  the  interior  of  Water 
Conservation  Area  1,  and  on  defining  classes  which  elucidated 
the  transitions  between  vegetation  communities . 
Unsupervised  Classification 

A  subset  was  extracted  from  the  1987  merged  image  which 
contained  a  5000  ha  portion  of  the  interior  of  Water 
Conservation  Area  1  (see  Figure  2-2) .   This  subset  was  run 
through  an  unsupervised  clustering  routine  to  produce  a  28-class 
classification  (ERDAS  1991).   The  clustering  was  based  on  the 
values  in  the  four  spectral  data  bands  of  each  pixel:  green 
reflectance  (band  1),  red  reflectance  (band  2),  near  infrared 
reflectance  (band  3)  and  the  NDVI,   (a  vegetation  index  of  the 
relationship  between  bands  2  and  3) .   The  main  factors 
influencing  spectral  reflectance  in  Water  Conservation  Area  1 
are  the  vegetation  species,  the  density  of  vegetation,  and  the 
amount  of  water  present.   Conveniently,  these  factors  also 
define  vegetation  communities  and  the  gradations  between  them. 
It  was  expected,  for  example,  that  the  community  of  wet  prairie 
species  would  be  divided  between  several  classes  separated  by 
their  proportions  of  open  water,  wet  prairie  vegetation,  and 
sawgrass.   These  classes  represent  stages  along  a  successional 
progression  from  slough  to  wet  prairie  to  sawgrass. 
Class  Identifications 

The  28  unsupervised  classes  had  to  be  identified  and 
ranked  into  a  serai  progression.   This  was  accomplished  by 
comparing  the  classes  with  georef erenced  ground  observations, 
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with  the  aerial  photographs  of  Water  Conservation  Area  1  taken 
in  December  1991,  and  with  the  locations  of  classes  in  the  1987 
vegetation  map. 

The  18-class  1987  vegetation  map  was  made  as  part  of  an 
effort  to  investigate  relationships  between  water  quality  and 
quantity  and  vegetation  change  in  Water  Conservation  Area 
l(Richardson  et  al.  1990).   The  study  required  detailed  mapping 
of  vegetation  near  the  canals,  particularly  sawgrass  and 
cattail.   However,  to  model  the  landscape  of  the  interior  of  the 
Water  Conservation  Area  1,  a  map  of  the  canal  zone  was 
unnecessary.   A  map  showing  different  densities  of  wet  prairie, 
sawgrass,  brush  and  trees  was  needed  instead.   The  1987 
vegetation  map  was  unsuitable  for  this  purpose,  and  thus  for  use 
in  the  landscape  model.   However  it  was  a  valuable  tool  used  in 
identifying  the  classes  of  the  unsupervised  classification. 

Seven  of  the  18  classes  in  the  1987  vegetation  map  did  not 
appear  in  the  5000-ha  area  of  Water  Conservation  Area  1  used  in 
the  unsupervised  classification  (see  Table  2-4) .   These  classes 
were  not  used  in  this  task.   The  11  classes  found  in  the  5000  ha 
subset  are  described  in  Table  2-5. 

The  unsupervised  classification  was  cross-referenced  with 
the  1987  vegetation  map  using  ERDAS  software's  SUMMARY  program 
(ERDAS  1991).   For  each  unsupervised  class,  the  program 
summarized  what  the  pixels  in  that  class  were  classified  as  on 
the  1987  vegetation  map.  By  analyzing  the  locations  of  the 
unsupervised  classes,  field  notes,  the  aerial  photography,  and  a 
nearest  neighbor  analysis  of  the  vegetation  map  (Richardson  et 
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al.  1990:46),  a  set  of  rules  for  identifying  the  unsupervised 
classes  based  on  the  summary  results  was  defined.   The  rules 
specified  the  successional  class  of  the  unsupervised  classes, 
based  on  the  percentage  of  each  class  in  the  vegetation  map 
cross-referenced  to  it  (see  Table  2-6) . 

The  28  classes  of  the  unsupervised  classification  were 
identified  according  to  the  rules,  renamed,  and  ranked  along  a 
successional  continuum  (see  Table  2-7).   This  created  a 
successional  class  map  suitable  for  use  as  a  base  map  for  the 
spatial  model  of  landscape  pattern  (see  Chapter  4) . 


Table  2-4.   Classes  in  the  1987  vegetation  map  of  Water 
Conservation  Area  1  found  in  less  than  1%  of  the  successional 
class  map  area. 


Class 
Number 

Description 

1 

High  density  sawgrass  with  some  fern  tussocks 

2 

Sawgrass  with  invasion  of  cattail  near  canal 

10 

Cattail  community  close  to  canal 

11 

Open  deep  water  mostly  along  Hillsborough  canal 

12 

Slough  deeper  or  less  vegetated  than  class  9 

15 

Sawgrass  with  invasion  of  wax  myrtle  close  to  canal 

18 

Cattail  community 

Source:  Richardson  et  al.  1990:40-41 
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Table  2-5.   Description  of  classes  in  the  1987  vegetation  map  of 
Water  Conservation  Area  1  found  in  the  successional  class  map 
area 


Class 
Number 

Description 

3 

Sawgrass-  ...primary  sawgrass  class  occurring  on 
all  parts  of  the  refuge  including  the  vast 
sawgrass  areas  on  the  west  side  of  the  refuge 

4 

Brush/sawgrass-  ...primarily  sawgrass  with  large 
amounts  of  wax  myrtle.  Some  tree  islands  which  may 
have  been  burned  out  previously  are  made  up 
entirely  of  the  class,  particularly  in  the 
southern  part  of  the  refuge 

5 

Tree  island-  ...lower  stature  tree  island 
community  made  up  of  a  mix  of  wax  myrtle,  dahoon 
holly  and  red  bay.  Occurs  along  edges  of  large 
tree  islands  with  some  small  islands  composed 
entirely  of  this  class 

6 

Wet  prairie-  ...largest  area  of  the  refuge.   This 
class  is  the  denser  wet  prairies  occurring  over 
all  of  the  refuge  but  the  primary  community  type 
of  the  central  portion  of  the  refuge.  This  class 
may  contain  small  tree  islands  smaller  than  30 
feet  across  and  areas  with  small  sawgrass  strands 

7 

Tree  island-  ...core  of  larger  tree  islands, 
larger  stature  trees  made  up  primarily  of  dahoon 
holly  and  red  bay  with  lesser  amounts  of  way 
myrtle  than  class  5 

8 

Brush-  ...smaller  brush  clumps  primarily  in  wet 
prairie 

9 

Wet  prairie-  ...sparser  wet  prairie  community 
often  with  sparse  sawgrass 

13 

Willow/brush-  ...predominantly  willow  but  with 
some  mixed  classification  whit  wax  myrtle  brush 
areas 

14 

Brush/tree  island-  ...class  of  many  of  the  smaller 
tree  islands  with  mostly  wax  myrtle,  some 
sawgrass,  occasionally  dahoon  holly  and  red  bay 

16 

Sawgrass-  ...slightly  higher  elevation  sawgrass 
than  class  3.   Core  of  sawgrass  ridges 

17 

Willow- ...  willow  along  canal  edge,  some 
misclassif ied  floating  aquatics  along  Hillsborough 
canal 

Source:  Richardson  et  al.  1990:40-41 
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Table  2-6.   Rules  fo 
on  the  1987  vegetati 


r  successional  class  identifications  based 
on  map  of  Water  Conservation  Area  1 


Vegetation  Type 

Rule  for  Identification  Based  on  Cross- 
Reference  with  Vegetation  Map  by  Richardson 
et  al . 

1.  Slough 

Type  is  >=50%  sparse  wet  prairie  (class  9), 
and  >=99%  sparse  wet  prairie  and  wet 
prairie  classes  combined  (classes  6  +  9) . 
Classes  are  ranked  in  descending  order  by 
percentage  of  sparse  wet  prairie  (class  9). 

2 .  Sparse  wet 
prairie 

Type  is  <50%  sparse  wet  prairie  (class  9), 
but  >=85%    sparse  wet  prairie  and  wet 
prairie  classes  combined  (classes  6+9). 
Classes  are  ranked  in  descending  order  by 
percentage  of  sparse  wet  prairie  (class  9) . 

3.  Dense  wet 
prairie  with 
sawgrass  and/or 
sparse  shrubs 

Type  is  >=50%  but  <85%  sparse  wet  prairie 
and  wet  prairie  classes  combined  (classes  6 
+  9) .  Sparse  wet  prairie  (class  9)  is  <  5%. 
Sawgrass  (classes  3  +  4  +  16)  and  brush 
(classes  8  +  14)  are  <  30%. 

Classes  are  ranked  in  descending  order  by 
percentage  of  wet  prairie  (class  6) . 

4 .  Sawgrass  and 
shrubs  with  some 
wet  prairie 

Type  is  >=25%  but  <50%  wet  prairie  (class 
6)  . 

Sawgrass  (classes  3  +  4  +  16)  and  brush 
(classes  8  +  14)are  present.   No  trees 
(classes  5  +  7  +  13  +  17)  are  present. 
Classes  are  ranked  in  descending  order  by 
percentage  of  wet  prairie  (class  6) . 

5 .  Shrubs 

Type  is  <25%  wet  prairie  (class  6) .  Brush 
is  >  =  40%  and  trees  are  present  but  <10% 
(classes  5  +  7  +  13  +  17) . 

Classes  are  ranked  in  descending  order  by 
percentage  of  wet  prairie  (class  6) . 

6.  Dense  shrubs  and 
trees 

Type  is  >=10%  and  <50%  trees  (classes  5+7 
+  13  +  17)  . 

Classes  are  ranked  in  descending  order  by 
percentage  of  wet  prairie  (class  6) . 

7 .  Low  stature 
trees 

Type  is  >=50%  non-core  trees  (classes  5  + 
13  +  17)  . 

8 .  Island  core 
trees 

Type  is  >=50%  core  trees  (class  7) .  Ranked 
in  ascending  order  by  percentage  of  core 
trees  ( class  7 )  . 

Source:  Richardson  et  al.  1990:40-41 
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Table  2-7.   Classes  in  the  successional  map  of  Water 
Conservation  Area  1 


Classes 

Vegetation  type 

1-3 

Slough 

4-9 

Sparse  wet  prairie 

10-14 

Dense  wet  prairie  with  sawgrass  and/or 
sparse  shrubs 

15-19 

Sawgrass  and  shrubs  with  some  wet  prairie 

20-23 

Shrubs 

24-25 

Dense  shrubs  and  trees 

26 

Low  stature  trees 

27-28 

Island  core  trees 

CHAPTER  3 
ANALYSIS  OF  TREE  ISLAND  DISTRIBUTION 


Introduction 

This  chapter  describes  a  set  of  analyses  of  the 
distribution  of  tree  islands,  which  led  to  the  development  of  a 
theory  of  landscape  pattern  dynamics  in  Water  Conservation  Area 
1.   Rather  than  attempt  an  analysis  of  the  pattern  of  all 
vegetation  types  in  Water  Conservation  Area  1  at  once,  tree 
islands  were  chosen  as  indicators  of  the  landscape  pattern. 
There  were  several  reasons  for  this  choice.   First,  trees  appear 
clearly  in  imagery  and  photography.   The  probability  of 
confusing  trees  with  slough,  wet  prairie,  or  sawgrass  is  low. 
Second,  because  tree  islands  have  abrupt  boundaries  and  the 
highest  elevation  relative  to  other  vegetation,  they  are  less 
subject  to  change  in  size  due  to  water  level  variation. 
Finally,  the  striking  visual  effect  of  the  pattern  of  tree 
islands  in  Water  Conservation  Area  1  was  the  inspiration  for 
this  body  of  work. 

The  analyses  took  two  different  approaches  (see  Figure  3- 
1).   In  the  first  approach,  the  distribution  of  different-sized 
tree  islands  found  in  a  satellite  image  of  the  landscape  of 
Water  Conservation  Area  1  was  examined.   The  number  and  sizes  of 
tree  islands  in  an  18-class  vegetation  map  classified  from  a 
1987  satellite  image,  (described  in  Chapter  2),  were  tallied 
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Figure  3-1.   Diagram  of  GIS  landscape  analysis 


using  a  geographic  information  system  (GIS) .   The  size 
distribution  of  the  islands  was  evaluated  in  light  of  expected 
distributions  according  to  hierarchy  and  other  ecological 
theories.   In  the  second  approach,  several  simple  population 
models  written  in  Basic  programming  language  were  used  to 
generate  size  distributions  using  different  general  processes. 
The  model-generated  distributions  were  compared  with  the  tree 
island  size  distribution  to  identify  which  mathematical 
functions  generated  a  similar  pattern. 
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The  information  gained  from  the  GIS  analysis  and  basic 
models  was  synthesized  with  the  landscape  data  described  in 
Chapter  2,  to  develop  a  conceptual  model  of  vegetation 
succession  and  tree  island  growth.   This  conceptual  model  is 
described  in  Chapter  4 . 

GIS  Analysis  of  Tree  Island  Distribution 

Analysis  of  Size  Distribution 

As  described  in  Chapter  2,  a  merged  1987  SPOT  satellite 
image  of  Water  Conservation  Area  1  was  studied  at  length,  both 
in  the  field  and  the  laboratory.   A  print  of  the  image  was  taken 
into  the  field  and  used  as  a  navigation  tool  and  a  guide  for 
making  ground  observations.   A  vegetation  map  classified  from 
this  image,  (Richardson  et  al.  1990),  was  analyzed  using  a 
geographic  information  system  (ERDAS  1991)  to  quantify  the 
distribution  of  tree  islands  in  Water  Conservation  Area  1.   The 
vegetation  map  was  composed  of  18  vegetation  classes,  and  was 
rectified  to  State  Plane  coordinates  using  a  30  ft  by  30  ft 
pixel  resolution. 
Methods 

A  4,440  ha  area  of  the  north-central  portion  of  the 
vegetation  map  was  analyzed  for  tree  island  distribution.   The 
area  was  chosen  because  it  contained  a  high  concentration  of 
islands  and  was  free  from  the  influence  of  the  perimeter  canals. 
Three  of  the  18  vegetation  classes  in  the  map  were  used  to 
delineate  tree  islands:  class  5  ("low  stature  tree")  class  7 
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("high  stature  tree")  and  class  8  ("brush").   Within  this  area 
all  contiguous  tree  island  pixels  in  classes  5,  7  and  8  were 
grouped  as  a  single  class  and  clumped  together  into  individual 
patches  using  ERDAS  software's  (version  7.5)  CLUMP  program 

(ERDAS  1991) .   Pixel  counts  were  made  for  each  patch,  and  these 
data  were  imported  into  a  spread  sheet. 

Because  the  islands  numbered  in  the  thousands,  every 
island  could  not  be  graphed.   Islands  were  grouped  by  size,  and 
the  data  were  graphed  according  to  size  class.   Because  of  a 
scaling  phenomenon,  small-size  classes  contained  many  islands, 
while  large-size  classes  contained  just  one  island  each.   Small- 
size  classes  were  one  pixel  different  from  each  other,  while  the 
largest  size  classes  were  hundreds  of  pixels  apart  in  size.   For 
this  reason,  the  data  were  graphed  in  rank  order  by  size,  rather 
as  a  histograms  with  pre-defined  size  intervals  on  the  x  axis. 
Histograms  can  obscure  the  pattern  of  data,  smoothing  out  or 
artificially  generating  discontinuities,  depending  on  the 
interval' s  size  in  relation  to  the  resolution  of  the  data 

(Holling  1992,  Gunderson  1992).   However,  it  is  important  to 
note  that  all  satellite  data  are  grouped  into  pre-defined  size 
intervals,  since  they  are  measured  in  discrete  pixels.   For 
example,  it  is  very  unlikely  any  10-pixel  islands  would  measure 
exactly  1000  m2  if  surveyed  using  more  precise  methods  on  the 
ground.   In  effect  the  rank  order  graphs  in  this  work  are  still 
histograms,  with  a  size  category  increment  of  one  pixel,  and 
with  empty  categories  not  shown. 


Results 

The  frequency  of  the  counted  islands  are  graphically 
displayed  by  size  in  Figure  3-2.   The  frequency  distrubution  is 
shown  in  more  detail,  with  only  islands  greater  than  10  pixels 
(836  m2)  in  Figure  3-3.   A  total  of  3623  islands  were  counted. 
An  unknown  number  of  islands  smaller  than  the  resolution  of  the 
satellite  image  could  not  be  counted.   The  island  density  was  81 
islands/km2,  less  than  the  99  islands/km2  sampled  in  the  1952 
aerial  photographs  (see  Chapter  2) .   The  smallest  islands 
detected  were  a  single  pixel  in  size  (90  ft^  or  83.6  m^ ) ,  while 
the  largest  island  detected  was  7522  pixels  (nearly  63 
hectares ) . 

The  sizes  of  islands  found  within  the  counted  area,  ranked 
from  smallest  to  largest,  are  shown  in  Figure  3-4.   The  mean 
island  area  was  43  pixels  (3595  m2)  .   The  median  island  area  was 
5  pixels  (418  m2)  .   At  the  lower  end  of  the  range  of  island 
sizes,  islands  were  found  at  every  possible  size  increment.   In 
other  words,  size  classes  of  islands  were  separated  by  one 
pixel.   However,  above  the  366-pixel  size  class  (306  hectares), 
the  gaps  between  the  size  classes  became  larger  than  one  pixel, 
increasing  along  with  island  size.   Above  the  2238-pixel  size 
class  (1872  hectares),  the  gaps  between  size  classes  became  much 
larger.   A  portion  of  this  distribution  is  shown  in  more  detail 
in  Figure  3-5 . 

The  total  area  of  all  the  islands  found  within  each  size 
class  is  shown  in  Figure  3-6.   This  graph  combines  both 
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frequency  and  size  information.   For  example,  the  465  islands  in 
the  2-pixel  size  class  have  a  total  area  of  930  pixels. 
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Figure    3-2.       Tree    island   distribution   by   size 
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Figure    3-3.      Number    of    islands    in   each    size    class    (detail) 
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Figure    3-4.       Island    size    class    distribution 
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Figure  3-5.  Island  size  classes  (detail) 
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Total  area  of  all  islands  in  each  size  class 
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Figure  3-6.   Total  area  of  islands  within  each  size  class 


Discussion 

The  number  of  tree  islands  in  this  portion  of  Water 
Conservation  Area  1  declined  exponentially  with  increasing  size 
(see  Figure  3-2) .   Approximately  half  of  the  islands  were  4 
pixels  or  less,  and  one  third  of  the  islands  were  2  pixels  or 
less.   Numbers  stabilized  with  one  island  per  size  class  at  70 
pixels  and  larger  (see  Figure  3-3) . 

The  frequency  distribution  of  islands  by  size  class  did 
not  resemble  a  Poisson  distribution.   A  Poisson  distribution 
would  indicate  a  random  process  determining  island  size,  as 
would  be  found  in  raindrops  splattered  on  a  surface.   It  instead 
resembled  a  frequency  distribution  typically  found  in  a  biotic 
population  with  high  juvenile  mortality. 

Discontinuities  in  the  size  distribution  of  entities  in  an 
ecosystem  may  indicate  the  influence  of  different  processes  on 


the  landscape  (Holling  1992,  Gunderson  1992) .   Changes  in  the 
slope  or  mathematical  function  of  the  distribution  at  some  size 
threshold (s)  may  mark  the  domain  boundaries  of  different 
processes.   Similarly,  the  age  or  size  distribution  of  a  biotic 
population  may  show  discontinuities  at  points  in  the  life  cycle 
where  different  processes  (such  as  predation  or  reproduction) 
begin  or  cease  to  be  influential.   If  small  islands  and  large 
islands  in  Water  Conservation  Area  1  are  formed  by  different 
processes,  (as  suggested  by  Loveless  1959  and  Gleason  and  Stone, 
1991:185),  we  would  expect  the  frequency  or  size  distribution  of 
small  and  large  islands  to  exhibit  different  characteristics. 

Such  a  discontinuity  was  suggested  in  the  size  class 
distribution  of  tree  islands  (see  Figure  3-3) .   This 
discontinuity  is  not  necessarily  evidence  of  two  separate 
processes  of  island  genesis;  other  processes  influencing  island 
structure  which  are  operating  at  different  scales  or  were  active 
during  different  periods  in  the  past  could  have  had  same  effect. 
Some  processes  of  island  expansion  or  reduction  that  might 
affect  islands  of  different  sizes  differently  are  vegetation 
succession,  peat  accretion,  erosion  by  sheet  flow  of  water, 
faunal  activity,  and  vegetation  disturbance  by  wind  and  fire. 
Some  of  these  processes  operate  most  actively  at  the  edges  of 
islands  (succession,  peat  accretion,  and  erosion) .   Because  the 
ratio  of  edge  to  island  area  decreases  with  increasing  island 
size,  the  effects  of  these  processes  may  change  with  islands 
size.   However,  this  change  should  be  gradual  and  not  cause  a 
discontinuity  in  the  size  distribution.   Other  processes,  (wind, 
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fire,  and  faunal  activity) ,  may  have  different  effects  above  and 
below  some  threshold  related  to  island  size.   A  threshold- 
triggered  process  such  as  fire  could  produce  a  discontinuity  in 
the  distribution  of  island  sizes. 

One  theory  that  addresses  the  relationship  between 
processes  and  spatial  extent  is  hierarchy  theory.   A  key 
principle  of  hierarchy  theory  (O'Neill  et  al.  1986,  Alan  and 
Starr  1982)  states  that  higher  order  processes  operate  at  larger 
spatial  scales  and  over  longer  time  scales  than  lower  level 
processes.   Higher  order  processes  control  or  organize  lower 
level  subsystems. 

While  ecologists  have  defined  hierarchical  levels  in  terms 
of  temporal  and  spatial  scales,  H.  T.  Odum  (1983)  described 
hierarchies  based  on  ecosystem  energetics.   Hierarchical  levels 
are  defined  by  energy  quality.   The  quality  of  energy  is 
measured  by  the  amount  of  solar  energy  embodied  in  it.   For 
example,  in  a  simple  hierarchy  of  plant,  herbivore  and 
carnivore,  energy  quality  increases  at  each  higher  trophic 
level.   A  single  calorie  of  plant  material  requires  10  calories 
of  solar  energy  for  its  manufacture.   To  make  a  single  calorie 
within  an  herbivore  requires  10  calories  of  plant  energy,  or  the 
equivalent  of  100  calories  of  solar  energy.   The  herbivore 
therefore  has  a  higher  quality  energy  than  a  plant. 

Odum' s  exponential  theory  of  energy  hierarchy  predicts 
that  while  energy  quality  increases  at  each  higher  hierarchical 
level,  the  amount  of  calories  present  (not  embodied)  in  each 
level  decreases  (See  Figure  3-7). 
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Since  each  energy  transformation  step  requires  utilization 
and  dispersion  of  potential  energy,  the  total  energy  flows 
decrease  as  one  passes  work  from  step  to  step. 
Conversely,  the  embodied  energy  per  unit  of  actual  energy 
increases  down  the  chain.  .  .  .   The  less  actual  energy 
remains,  the  more  embodied  energy  has  been  involved. 
(Odum  1983:270) 

Energy  quality  and  hierarchical  level  also  relate  to  the 

spatial  scale  of  processes  or  entities: 

Flows  of  low-quality  energy  are  abundant  and  widely 
dispersed,  and  individual  units  [receiving  that  energy] 
are  small  in  size.   Higher-quality  units  and  their  flows, 
although  less  in  total  energy  flow,  are  more  concentrated 
and  each  unit  is  larger  in  size  with  a  larger  territory 
from  which  it  receives  energy  and  feeds  back  its  actions. 
(Odum  1983:15) 

Examples  of  this  phenomenon  are  the  trophic  structure  of 

ecosystems,  where  primary  producers  outnumber  consumers;  the  age 

structure  of  stable  populations,  where  small,  young  individuals 

outnumber  larger  adults;  and  the  occurrence  of  earthquakes, 

where  weak  quakes  outnumber  strong  ones.   Tree  islands  in  Water 

Conservation  Area  1  are  distributed  exponentially  by  size,  with 

small  islands  outnumbering  large  islands.   Does  this  indicate  an 

energy  hierarchy  in  the  island  population? 

An  energy  hierarchy  would  not  apply  to  island  patches 

unless  larger  islands  received  inputs  of  higher  quality  energy 

than  smaller  islands.   Unlike  mice  versus  bobcats,  or  nestlings 

versus  adults,  large  islands  do  not  have  a  different  strategy 

for  obtaining  energy  than  small  islands.   If  the  amount  and 

quality  of  energy  required  to  create  one  unit  area  of  island  is 

the  same  regardless  of  island  size,  then  all  islands  are  at  the 

same  hierarchical  level.   There  is  no  evidence  that  large 

islands  are  receiving  higher  quality  energy  than  smaller  islands 
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at  the  present  time.  It  is  more  likely  that  islands  are  at  the 
same  hierarchical  level,  and  a  higher  order  disturbance  process 
is  organizing  them  into  the  sizes  we  observe. 
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Figure  3-7.   Diagram  of  energy  hierarchy 
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Island  Population  Model  Analysis 


To  better  understand  the  observed  size  distribution  of 
tree  islands  in  Water  Conservation  Area  1,  simple  population 
models  were  written  to  generate  island  distributions  using 
different  growth  functions.   The  models  were  based  on  the 
assumption  that  islands  appear  and  grow  in  size  over  time  and  do 
not  "die."   The  models  were  used  to  illustrate  the  curves 
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produced  by  different  growth  functions  when  the  data  were 
graphed,  and  to  explore  the  interaction  of  the  scales  of  island 
growth  and  the  unit  used  to  measure  growth. 
Methods 

Two  simple  programs  were  written  in  Basic  programming 
language.   The  first  program  generated  a  population  distribution 
based  on  either  a  constant  growth  model  or  an  area-dependent 
growth  model,  (see  Figure  3-8).   The  program  generated  a  size 
distribution  by  "growing"  an  island  and  recording  its  size  at 
each  iteration  (time  step) .   It  assumed  that  conditions  for 
growth  remain  constant  over  time;  i.e.  a  small  island  in  the 
present  grows  the  same  way  a  large  island  did  in  the  past  when 
it  was  small.   It  also  was  assumed  that  new  islands  were  created 
by  peat  batteries  as  described  in  Chapter  1 .   The  number  of  new 
islands  introduced  to  the  population  could  be  altered  by  the 
program,  but  most  model  runs  were  made  with  a  constant  number  of 
islands  added  to  the  population  each  time  period,  to  narrow  the 
number  of  factors  influencing  the  distribution.   New  islands 
were  1  m2  in  size. 

The  second  program  generated  a  population  distribution 
using  a  growth  function  based  on  island  edge  (see  Figure  3-9) . 
Otherwise,  the  program  was  similar  to  the  one  described  above. 
The  model  assumed  islands  were  circular,  and  calculated  island 
edge  based  on  the  circumference  of  a  circle. 

Both  programs  rounded  the  output  data  into  integer  values, 
to  simulate  the  effect  of  measuring  island  areas  in  pixels. 
Pixels  in  a  satellite  image  are  classified  based  on  their 
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spectral  values.   The  spectral  values  of  a  pixel  integrate  the 
reflectance  of  all  the  vegetation  in  the  scanned  unit.   The 
dominant  vegetation  therefore  has  the  greatest  influence  on 
pixel  spectral  values.   If  the  amount  of  woody  vegetation  within 
a  scanned  unit  on  the  ground  reaches  a  certain  threshold,  a 
pixel  will  be  classified  as  an  island  class  rather  than  a  non- 
island  class.   The  area  of  woody  vegetation  cannot  be  measured 
in  fractions  of  a  pixel.   In  the  same  way,  the  area  values  for 
islands  in  these  basic  models  are  rounded  to  the  nearest 
integer . 

Models  were  run  using  different  coefficient  values, 
generally  for  1000  iterations.   Frequency  distributions  and 
graphs  were  generated  from  the  output  data  using  spreadsheet 
software . 
Results 

The  constant  growth  model  was  first  run  with  the  simplest 
scenario:  constant  recruitment  of  one  new  island  per  iteration 
(n=l)f  and  constant  incremental  island  growth  of  one  pixel  per 
time  period  (g=l).   This  function  is  analogous  to  the  following: 
every  island  has  one  alligator  hole  next  to  it,  regardless  of 
island  size.   The  activity  at  every  alligator  hole  generates  the 
same  amount  of  disturbed  peat,  which  drifts  to  the  adjacent 
island,  adding  to  its  area. 
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REM 

MODEL  TO  GENERATE  ISLAND  DISTRIBUTION-  INCREMENTAL  OR 

REM 

AREA-DEPENDENT  GROWTH 

REM 

i  =  iteration  (or  age) 

REM 

n  =  number  of  new  islands  per  iteration 

REM 

g  -  fixed  area  of  growth  added  per  iteration 

REM 

L  =  growth  coefficient  as  a  function  of  island  area 

REM 

s  =  size  class 

REM 

a  =  total  area  of  islands  in  size  class 

CLS 

OPEN  ' 

'c:\treebas\tree.dat"'  FOR  OUTPUT  AS  #1 

INPUT 

"enter  max  iteration  and  number  of  islands";  k,  n 

INPUT 

"enter  island  growth  increment  (g)  for  arithmetic  growth, 

or  enter  0";  g 

INPUT 

"enter  island  growth  coefficient  (L)  for  size-dependent 

growth, 

or  enter  1";  L 

PRINT 

#1,  "max  iteration  =  ",  k 

PRINT 

#1,  "number  of  new  islands  per  iteration  (n)  =",  n 

PRINT 

#1,  "island  growth  increment  (g)  =",  g 

PRINT 

#1,  "island  growth  factor  (L)  =",  L 

PRINT 

#1,  "Iteration       s              n               a   " 

IF  g  = 

=  0  THEN  s  =  1 

FOR  i 

=  1  TO  k 

s  =  (s  *  L)  +  g 

PRINT  #1,  CINT(i),  CINT(s),  CINT(n),   (CINT(n)  *  CINT(s)) 

PRINT  ,  CINT(i),  CINT(s),  CINT(n),   (CINT(n)  *  CINT(s)) 

NEXT  i 

CLOSE 

#1 

Figure  3-8.   Basic  model  code  to  generate  island  distributions: 
incremental  and  area-dependent  growth 
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REM    Model  to  generate  island  distribution-  edge-dependent 
growth 

REM  i  =  iteration  (or  age) 

REM  n  =  number  of  new  islands  per  iteration 

REM  s  =  size  class 

REM  c  =  circumference, 

REM  a  =  total  area  of  islands  in  size  class 

CLS 

OPEN  "c:\treebas\edgear.txt"  FOR  OUTPUT  AS  #1 

INPUT  "what  is  the  max  iteration";  k 

INPUT  "what  is  the  number  of  new  islands  per  iteration";  n 

INPUT  "what  is  the  initial  size  class";  s 

INPUT  "growth  factor  (x)";  x 

PRINT  #1,  "max  iteration  =  ",  k 

PRINT  #1,  "number  of  new  islands  per  iteration  (n)  =",  n 

PRINT  #1,  "island  growth  factor  (x)",  x 

PRINT  #1,  "Iteration       s  n  a   " 

PRINT  "  i",  "  s",  "  c",  "n",  "(n  *  s);  ";  "" 

c  =  (SQR(s  /  3.141592))  *  2    *  3.141592 
FOR  i  =  1  TO  k 

PRINT  #1,  CINT(i),  CINT(s),  CINT(n),  CINT(n  *  s) 

PRINT  i,  s,  c,  n,  (n  *  s) 

s  =  s  +  (x  *  c) 

c  =  (SQR(s  /  3.141592))  *    2    *    3.141592 
NEXT  i 
CLOSE  #1 


Figure  3-9.   Basic  model  code  to  generate  island  distributions 
edge-dependent  growth 
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Figure  3-10.   Distribution  from  incremental  growth:  island  size 
classes 
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Figure  3-11.   Distribution  from  incremental  growth  with  high 
constant  island  recruitment 
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This  function  generated  a  linear  size  distribution  (see 
Figure  3-10) .   There  was  just  one  island  per  size  class,  so  if 
the  total  area  per  size  class  were  graphed,  the  graph  would  be 
identical  to  Figure  3-10.   If  the  growth  increment  (g)  was 
increased,  the  function  generated  a  line  with  a  steeper  slope. 

When  more  islands  were  added  per  iteration  (n=10),  the 
size  distribution  remained  the  same,  but  the  total  area  within 
size  class  increased  at  a  steeper  slope  (Figure  3-11) . 

When  growth  is  incremental,  but  new  islands  are  added  in 
changing  numbers  rather  than  at  a  constant  rate,  the  total  area 
in  each  size  class  is  affected.   The  model  was  run  to  simulate 
island  recruitment  that  was  high  when  older  islands  were  formed 
and  decreased  in  a  linear  fashion  over  time  (n  =  n+1) .   Such  a 
situation  might  occur  if  battery  formation  decreased  over  time 
as  the  proliferating  islands  occupied  more  and  more  of  the  wet 
prairie  area  where  batteries  form.   Figure  3-12  shows  the 
results  of  this  model  run.   The  total  area  per  size  class  was  an 
exponential  function  (see  Figure  3-13). 

The  model  was  next  run  to  simulate  island  recruitment  that 
was  low  when  older  islands  were  formed  and  increased  for  younger 
age  classes,  (n  =  n-1) .   Such  a  situation  might  occur  if  long- 
term  peat  buildup,  changing  hydroperiods,  or  some  other  factor 
increased  the  rate  of  battery  formation  over  time.   Figure  3-14 
shows  the  results.   In  this  case  the  total  area  per  size  class 
was  quadratic  function  (Figure  3-15) . 
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Figure  3-12 .   Distribution  from  incremental  growth  with 
increasing  island  recruitment:  island  size  classes  and  numbers 
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Figure  3-13.   Distribution  from  incremental  growth  with 
increasing  island  recruitment:  total  area  in  each  size  class 
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Figure  3-14.   Distribution  from  incremental  growth  with 
decreasing  island  recruitment:  island  size  classes  and  numbers 
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Figure  3-15.   Distribution  from  incremental  growth  with 
decreasing  island  recruitment:  total  area  in  each  size  class 
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The  model  was  next  run  using  size-dependent  growth  instead 
of  constant  growth.   This  growth  function  would  occur  if  the 
rate  of  island  growth  increased  as  the  island  grew  larger. 
Size-dependent  growth  could  occur  if  the  probability  of  an 
island  intercepting  floating  peat  increased  with  increasing 
size,  or  a  greater  proportion  of  core  tree  area  speeded  the  rate 
of  island  growth. 

The  model  was  run  with  the  size-dependent  growth  function, 
constant  island  recruitment  (n=l),  and  size-dependent  growth 
equal  to  1%  of  the  island's  area  per  iteration  (L=1.01).   With 
the  size-dependent  growth  function,  a  scaling  effect  emerged,  as 
the  model  generated  islands  with  areas  that  were  not  an  integer 
value.   These  values  were  rounded  off  and  put  into  integer  bins. 
While  the  resolution  of  the  analysis  was  fixed  at  integer 
intervals,  the  scale  of  island  growth  changed  between  large  and 
small  islands.   When  islands  were  small,  the  amount  of  growth 
added  each  iteration  was  small  enough  that  several  iterations  in 
a  row  produced  islands  of  the  same  integer  size.   When  plotted, 
islands  from  multiple  iterations  were  grouped  into  the  same 
integer  size  class. 

This  scaling  effect  is  apparent  in  Figure  3-16.   Large 
numbers  of  islands  are  grouped  into  the  smallest  size  classes 
even  though  only  one  new  island  was  added  to  each  iteration. 
While  the  growth  function  of  this  model  produced  an  exponential 
curve  (Figure  3-17),  the  curve  was  forced  into  a  linear 
distribution  at  the  small  scale  end  of  the  scale  due  to  the 
integer  constraint  (Figure  3-18).   When  these  effects  are 
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Figure  3-17.   Distribution  from  area-dependent  growth:  island 
size  classes 
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Figure  3-18.   Distribution  from  area-dependent  growth:  iland 
size  classes  (detail) 
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Figure  3-19.   Distribution  from  area-dependent  growth:  total 
area  in  each  size  class 
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Figure  3-20.   Distribution  from  area-dependent  growth:  total 
area  in  each  size  class  (detail) 


combined  in  the  calculation  of  the  total  area  per  size  class, 
the  curve  is  clearly  perturbed  at  the  small  scale  (Figures  3-19 
and  3-20) . 

Finally,  a  model  using  growth  based  on  the  edge  length  of 
the  island  was  run,  using  constant  island  recruitment.   This 
growth  function  could  occur  if  the  peat-building  and  vegetative 
reproduction  of  sawgrass  at  the  perimeter  of  islands  were  the 
main  factor  in  island  growth  rather  than  the  accumulation  of 
floating  peat  and  the  litter  contributed  by  trees  and  other 
species  in  the  interior  of  islands.   The  growth  coefficient  (x) 
was  the  unit  area  added  per  unit  of  edge  length.   The  model 
assumed  islands  were  circular,  and  estimated  edge  length  each 
iteration  by  calculating  the  circumference  of  a  circle  with  the 
same  area  as  the  island. 

A  run  of  the  model  using  .14  units  of  area  added  per  unit 
of  edge  length  (x  =  .14)  generated  Figure  3-21.   The  curve  of 


Ill 


the  distribution  is  exponential,  but  more  gradual  than  the  area- 
dependent  curves.   The  scaling  effect  is  less  using  this 
function;  islands  were  grouped  only  in  the  five  smallest 
classes . 
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Figure  3-21.   Distribution  from  edge-dependent  growth:  island 
size  classes 


Discussion 

Of  all  the  distributions  generated  by  the  basic  models, 
the  size-dependent  growth  model  most  resembles  the  data  from  the 
GIS  analysis  described  earlier  in  this  chapter.   The  scaling 
phenomenon  illustrated  by  this  size-dependent  growth  model  also 
appears  to  be  affecting  the  size  class  distribution  in 
Richardson  et  al.'s  map,  (see  Figures  3-3  and  3-6) .   In  both 
cases  there  are  high  numbers  of  islands  in  the  smaller  size 
classes,  (Figures  3-3  and  3-16),  only  one  island  apiece  in  the 
larger  classes,  and  the  difference  between  the  island  size 
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classes  increases  as  the  islands  get  larger  (Figures  3-4  and  3- 
17)  . 

This  scaling  phenomenon  is  important  to  consider  when 
using  the  sizes  and  numbers  of  entities  as  indicators  of 
ecosystem  process  domains,  especially  if  satellite  imagery  is 
used.   The  basic  model  of  area-dependent  growth  illustrates  how 
a  single  spatially-dependent  function  producing  a  flat  frequency 
distribution  (one  island  per  iteration) ,  appears  to  produce  an 
exponential  frequency  distribution  when  measured  at  too  coarse  a 
resolution.   This  frequency  distribution,  when  combined  with  the 
island  size  distribution,  exhibited  noisy  variation  at  the  low 
end  of  the  scale.   This  variation  could  be  misinterpreted  as  a 
discontinuity  caused  by  the  influence  of  a  second  process 
(Figures  3-19  and  3-20) .   The  distribution  of  the  islands  in 
Richardson  et  al.'s  map  shows  similar  variation  at  the  low  end 
of  the  scale  (Figure  3-6) .   This  is  also  likely  to  be  a  scaling 
phenomenon. 


CHAPTER  4 
LANDSCAPE  MODEL 


Introduction 


This  chapter  describes  two  models  of  the  landscape  of 
Water  Conservation  Area  1:  a  conceptual  model  of  island 
formation  and  vegetation  community  change,  and  a  spatial  model 
of  landscape  pattern.   The  conceptual  model  is  a  synthesis  of 
the  information  presented  in  the  first  three  chapters  of  this 
work,  made  to  identify  the  most  important  factors  in  island 
formation  and  vegetation  community  change.   From  the  conceptual 
model,  a  spatial  model  with  a  narrower  scope  was  developed.   The 
spatial  model  was  used  to  project  the  relationship  between 
landscape  pattern  and  disturbance  by  fire  in  Water  Conservation 
Area  1  over  time.   The  model  simulated  vegetation  succession, 
the  effect  of  fire  on  vegetation,  and  post-disturbance 
vegetation  succession.   The  landscape  patterns  generated  by  the 
model  expressed  the  interaction  of  these  processes. 

Conceptual  Model  of  Tree  Island  Growth  and  Landscape  Pattern 

A  conceptual  model  of  island  growth  and  vegetation 
community  change  in  Water  Conservation  Area  1  was  developed  as  a 
precursor  to  a  spatial  model  of  landscape  pattern.   The 
conceptual  model  united  the  many  pieces  of  information  gathered 
in  the  first  three  chapters  of  this  work:  ecological  data 
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gathered  and  analyzed  by  other  researchers,  landscape  data 
collected  both  remotely  and  in  the  field,  and  data  analyses 
aimed  at  a  theoretical  understanding  of  landscape  pattern  (see 
Figure  4-1) .   The  synthesis  of  this  information  led  to  the 
conceptualization  of  a  system  that  could  generate  and  maintain 
the  landscape  pattern  we  observe  in  Water  Conservation  Area  1. 
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Figure  4-1.  Information  used  to  develop  the  conceptual  model  of 
vegetation  change  and  island  dynamics  in  Water  Conservation  Area 
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This  conceptual  model  is  described  below.   The  description 
summarizes  information  already  presented  and  referenced  earlier 
in  this  work.   For  more  complete  information,  the  reader  is 
referred  to  these  earlier  sections. 

Processes  Contributing  to  Island  Growth 

The  conceptual  model  was  limited  to  the  main  processes 
contributing  to  island  growth  and  to  the  pattern  of  all 
vegetation  communities  in  the  interior  of  Water  Conservation 
Area  1.   These  processes  are  fire,  peat  accumulation,  and 
vegetation  succession.   These  three  processes  are  in  turn 
affected  by  a  number  of  factors,  the  majority  of  which  are 
driven  by  climate  and  water  management  (see  Figure  4-2) . 
Vegetation  succession 

Vegetation  succession  is  a  positive  force  in  the  formation 
and  growth  of  tree  islands.   Vegetation  succession  is  the 
process  by  which  land  passes  through  a  succession  of  vegetation 
community  types  over  time  (serai  stages)  (Chapter  1:42) .   Each 
serai  stage,  as  it  matures,  alters  the  conditions  of  its 
environment  until  the  environment  is  more  favorable  to  the  next 
serai  stage.   Since  tree  islands  in  Water  Conservation  Area  1 
are  defined  by  the  presence  of  the  bayhead  vegetation  community 
(Chapter  1:25),  vegetation  succession  is  critical  to  the 
dynamics  of  tree  islands.   In  Water  Conservation  Area  1,  the 
environmental  conditions  driving  succession  are  water  level  and 
hydroperiod  (Chapter  1:41).   Water  level  and  hydroperiod  have 
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Diagram  of  processes  contributing  to  tree  island 


positive  or  negative  effects  on  succession.   Lower  water  level 
and  shorter  hydroperiod  move  succession  forward  from  herbaceous 
vegetation  to  woody  vegetation.   Higher  water  level  and  longer 
hydroperiod  move  succession  backward  from  woody  vegetation  to 
herbaceous  vegetation.   Vegetation  itself  influences  water  level 
and  hydroperiod,  mainly  by  building  peat  (Chapter  1:20).   Peat 
accumulation  moves  succession  forward  by  raising  elevation, 
which  reduces  water  depth  and  hydroperiod.   Hurricanes  and  fires 
move  succession  backwards  by  destroying  or  damaging  vegetation. 

Plant  reproduction  is  essential  to  the  process  of 
succession,  whether  by  seed  or  vegetative  reproduction.   Changes 
in  seed  sources,  nutrient  levels,  and  the  ability  of  plants  to 
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reproduce  vegetatively  can  have  a  positive  or  negative  effect  on 
succession,  changing  the  species  composition  of  plant 
communities  and  the  competitive  balance  between  species.   These 
changes  may  redefine  the  series  of  vegetation  communities  in  a 
successional  sequence. 
Peat  accumulation 

Peat  accumulation  is  a  positive  force  in  the  formation  and 
growth  of  tree  islands.   Peat  accumulates  on  islands  via  litter 
production  or  drifts  in  from  outside  the  island.   When  litter 
from  island  vegetation  accumulates  faster  than  it  decomposes,  it 
builds  peat.   Alligator  activity  disturbs  peat,  which  can  drift 
and  collect  at  island  margins  (Chapter  1:28) .   Peat  batteries 
may  float  and  resettle,  creating  mounds  that  may  become  islands 

(Chapter  1:26) .   Sheet  flow  of  water  across  the  landscape  has  a 
positive  or  negative  effect  on  island  growth:  eroding, 
transporting  or  depositing  peat  (Chapter  1:22). 
Fire 

Fire  is  a  negative  force  in  the  formation  and  growth  of 
tree  islands.   Fire  destroys  or  damages  vegetation  and  can 
reduce  elevation  by  burning  peat,  moving  succession  backwards 

(Chapter  1:33).   The  probability  of  a  fire  and  the  intensity  of 
fire  is  determined  by  fire  potential.   Fire  potential  is 
determined  by  a  variety  of  factors,  including  the  density  and 
moisture  content  of  the  fuel  load,  atmospheric  humidity,  wind 
speed,  and  soil  moisture  (Chapter  1:30).   Fuel  moisture  and  soil 
moisture  are  a  function  of  water  level,  which  is  determined  by 
water  management  and  the  drought  cycle.   Atmospheric  humidity 


118 


and  wind  speed  are  a  function  of  current  weather  conditions  and 
the  microclimate  of  a  site  (Chapter  1:36) . 

Fire  potential  combines  with  a  source  of  ignition  to 
produce  a  fire.   Ignition  may  come  from  a  human  source  or  from 
lightning.   If  fire  potential  is  low,  ignition  will  produce  only 
a  limited  fire  affecting  less  than  a  hectare.   If  fire  potential 
is  very  high,  fires  may  propagate  quickly,  affecting  tens  of 
thousands  of  hectares  (Chapter  1:31)  . 

Fire  potential  in  Water  Conservation  Area  1  is  high  when 
the  water  level  is  low,  particularly  at  the  end  of  the  dry 
season.   During  droughts,  the  conditions  in  Water  Conservation 
Area  1  are  especially  favorable  for  the  initiation  and 
propagation  of  fire.   With  high  evapotranspiration  and  no 
rainfall,  the  water  level  falls  below  the  peat  surface.   This 
allows  a  shift  in  the  species  composition  and  density  of  the  wet 
prairies.  Panicum   hemitomon    and  annual  herbs  sprout  among  the 
sparse  cover  of  Eleocharis   and  Rhynchospora    species,  increasing 
the  fuel  load  to  a  level  capable  of  carrying  a  fire.   When 
ignition  occurs,  fire  can  then  sweep  across  the  landscape  (Wade 
et  al.  1980) . 

Vegetation  and  microclimatic  variables  that  increase  fire 
potential —  fuel  moisture  content,  understory  vegetation 
density,  wind  speed,  temperature  and  humidity--  are  lower  under 
tree  island  canopies  than  in  adjacent  vegetation  communities 
(Chapter  1:31)  .   As  observed  during  a  large  severe  fire  in  Water 
Conservation  Area  1,  fire  can  sweep  across  the  landscape,  and 
slow  when  it  reaches  the  lower-fire-potential  large  tree 
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islands.   The  fire  front  moves  around  rather  than  through  the 
island.   While  the  front  continues  on,  islands  burn  from  the 
margins  inward  until  the  moister  interior  conditions  extinguish 
the  fire.   While  small  islands  do  not  have  enough  interior  area 
to  generate  a  wetter  microclimate,  the  cores  of  large  islands 
are  protected  from  fire. 
Tree  islands 

The  tree  islands  in  Water  Conservation  Area  1  are  composed 
of  peat,  and  are  built  up  in  the  process  of  vegetation 
succession   (Chapter  1:26) .   Tree  islands  may  form  on  peat 
batteries  displaced  in  wet  prairie,  or  they  may  be  created  more 
slowly  by  litter  accumulation,  peat  buildup,  and  succession  from 
lower  serai  stage  vegetation. 

If  the  battery-island  theory  is  correct,  initially  a  peat 
battery,  alligator,  or  other  disturbance  creates  a  small 
topographic  high  in  the  wet  prairie.   The  reduced  hydroperiod  at 
this  higher  elevation  allows  colonization  by  sawgrass  {Cladium 
jamaicense) .   Senescent  biomass  from  the  sawgrass  accumulates 
and  builds  up  the  peat.   In  the  absence  of  disturbance  by  fire, 
colonization  by  woody  plants  follows.   First,  shrubs  {Myrica 
cerifera    and  Cephalanthus    occidentalis)     create  a  brush  island. 
Then  trees,  (dominated  by  Persea   borbonia    and  Ilex    cassine) , 
become  established  and  create  a  tree  island. 

Tree  islands  could  also  develop  without  a  battery  as  a 
catalyst.   With  time  and  the  absence  of  disturbance,  vegetation 
continuously  builds  peat.   If  elevation  was  lowered  by 
disturbance  in  one  area,  adjacent  undisturbed  areas  would  be 
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left  relatively  higher  and  drier.   The  process  of  succession 
would  then  progress  as  on  a  displaced  peat  battery,  forming  a 
tree  island. 

As  succession  proceeds  over  time,  tree  islands  spread 
outward,  increasing  in  area.   On  established  tree  islands,  dead 
biomass  accumulates  from  litter  originating  on  the  island  and 
from  floating  debris .   This  accumulation  raises  elevation  at  the 
island's  perimeter.   Narrow  bands  of  sawgrass  and  brush  ring  the 
perimeter  of  tree  islands,  where  sawgrass  has  spread  outward 
into  wet  prairie,  and  brush  has  colonized  the  sawgrass  (Figure 
2-3).   This  configuration  suggests  that  tree  islands  grow  in 
area  as  a  function  of  time. 

Conceptual  Model 

The  many  processes  contributing  to  tree  island  growth 
listed  in  Figure  4-2  were  narrowed  into  a  simpler  and  more 
specific  conceptual  model  (see  Figure  4-3) .   The  model  is  an 
abstraction  of  the  dynamics  of  vegetation  communities  in  the 
interior  of  Water  Conservation  Area  1.   This  conceptual  system 
cycles  through  three  phases:  1)  vegetation  succession  and  the 
growth  of  tree  islands;  2)  fire  behavior;  and  3)  post-fire 
vegetation  recovery.   These  phases  are  described  below.   The 
landscape  processes  associated  with  the  canals  at  the  perimeter 
of  Water  Conservation  Area  1  (the  increase  of  nutrients,  water 
depth  and  hydroperiod)  are  not  addressed  in  the  model. 
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Figure  4-3.   Conceptual  model  of  vegetation  and  island  dynamics 
within  the  interior  of  Water  Conservation  Area  1 


Phase  1:  vegetation  succession 

In  the  absence  of  disturbance,  vegetation  succession 
follows  a  sequence  of  serai  stages  driven  by  hydrologic 
conditions,  as  indicated  by  the  plus-sign  track  in  Figure  4-3. 
Succession  moves  vegetation  to  the  right  along  this  track  (from 
wet  prairie  to  sawgrass,  to  brush,  and  to  trees),  and 
disturbance  by  fire  moves  it  to  the  left.   The  transition  from 
wet  prairie  to  sawgrass  and  brush  may  be  accelerated  by  the 
formation  of  peat  batteries. 


122 

Phase  2 :  fire 

When  conditions  are  favorable  for  fire  and  ignition 
occurs,  fires  will  burn  (indicated  by  the  dashed  line  in  Figure 
4-3) .   The  damage  to  vegetation  communities  by  a  fire  depends  on 
the  fire's  severity  and  the  vulnerability  of  plant  species  in 
each  community.   Damage  is  portrayed  in  Figure  4-3  by  a  series 
of  questions  and  pathways. 

As  shown  on  the  right-hand  side  of  Figure  4-3,  if  a  tree 
island  is  above  some  size  threshold,  the  trees  in  its  interior 
escape  burning  and  remain  intact.   On  tree  islands  below  this 
size  threshold,  and  at  the  margins  large  islands,  trees  are 
burned.   If  the  burn  kills  only  the  above-ground  portion  of  the 
trees,  the  trees  will  resprout  in  a  brushy  form,  creating  a 
brush  island.   With  the  absence  of  a  tree  canopy,  shrub  species 
will  also  germinate  on  the  brush  island.   If  the  burn  is  severe 
enough  to  cause  tree  mortality,  but  the  peat  itself  does  not 
burn  below  a  thin  layer  at  the  surface,  then  a  tree  island  will 
become  a  patch  of  sawgrass  and  ferns.   If  the  peat  burns  deeply, 
reducing  the  peat  elevation,  the  island  will  become  wet  prairie. 
This  is  a  rare  event  in  the  interior  of  Water  Conservation  Area 
1. 

Patches  of  shrubs  (brush  islands)  do  not  have  the  moist 
microclimate  of  large  tree  islands,  so  entire  shrub  patches  will 
burn  in  a  fire.   If  only  the  above-ground  portion  of  the  shrubs 
are  killed,  then  they  will  resprout,  keeping  the  patch  in  a 
shrub  community.   If  the  fire  causes  shrub  mortality  but  the 
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peat  does  not  burn,  then  the  patch  will  become  sawgrass.   If  the 
peat  burns  deeply,  the  patch  will  become  wet  prairie. 

Likewise,  if  a  sawgrass  patch  burns  but  the  fire  does  not 
burn  into  the  peat,  the  sawgrass  will  resprout.   If  it  does  burn 
down  into  the  peat,  causing  sawgrass  mortality,  then  the  patch 
will  become  wet  prairie. 
Phase  3:  post-fire  vegetation  recovery 

Vegetation  recovers  after  a  fire,  either  by  the  regrowth 
of  surviving  individuals  or  by  the  germination  of  new 
individuals.   In  the  interior  of  Water  Conservation  Area  1, 
where  environmental  conditions  have  changed  less  than  at  the 
perimeter  of  the  area,  the  recovery  of  vegetation  follows  the 
same  succession  track  in  shown  Figure  4-3. 

At  the  perimeter  of  Water  Conservation  Area  1,  changes  in 
the  environment  including  the  invasion  of  exotic  species, 
nutrient  enrichment,  and  alterations  of  hydroperiod,  have 
changed  the  successional  pathway  followed  after  a  fire  or  other 
disturbance.   This  conceptual  model  does  not  apply  to  the 
perimeter  area. 
Fire  and  landscape  pattern 

If  the  process  of  fire  disturbance  is  influenced  by  the 
amount  of  core  area  of  tree  islands  as  this  model  describes, 
then  fire  acts  to  organize  landscape  pattern.   If  tree  islands 
above  some  size  threshold  are  spared  the  burning  of  their 
interiors,  their  undisturbed  cores  will  promote  more  rapid 
successional  recovery  and  continued  expansion  than  that 
occurring  on  completely  burned  islands.   If  they  expand  at  a 
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rate  greater  than  their  setback  by  fire,  their  vulnerability  to 
fire  will  further  diminish  in  a  positive  feedback.   In  time, 
repeated  fires  will  create  a  population  of  tree  islands  divided 
between  islands  whose  growth  is  controlled  by  fire,  and  larger 
islands  with  a  core  impervious  to  fire  and  a  higher  growth  rate. 
Such  a  situation  could  produce  a  discontinuity  in  the 
distribution  of  tree  islands  similar  to  that  found  in  the  map  of 
Richardson  et  al.  (see  Figure  3-4) . 

This  theory  cannot  be  tested  in  the  field  in  Water 
Conservation  Area  1,  because  the  temporal  and  spatial  scales  of 
the  variables  are  too  great  to  measure  within  a  human  life-span 
and  a  reasonable  budget.   However,  a  spatial  model  can  be  used 
to  show  whether  this  dynamic  system  can  be  reproduced  in 
abstract . 

Spatial  Model  of  Vegetation  Succession  and  Disturbance  by  Fire 

A  spatial  model  of  the  landscape  dynamics  of  the  interior 
of  Water  Conservation  Area  1  was  developed  using  a  cartographic 
approach  and  a  geographic  information  system  (see  Chapter  1:12). 
Although  cartographic  models  have  usually  been  used  to  predict 
an  outcome  in  a  single  step,  this  model  used  a  geographic 
information  system  iteratively,  to  show  landscape  change  over  an 
extended  time  under  different  fire  regimes .   As  shown  in  Figure 
4-4,  the  foundation  of  the  model  was  a  map  of  the  vegetation  of 
Water  Conservation  Area  1  by  successional  class  (see  Chapter 
2:81-86) .   The  functions  of  the  model  were  based  on  the 
conceptual  model  in  Figure  4-3.   The  model  was  used  to  generate 


125 


landscapes    subjected    to    fires    at    different    intervals,    and    at 
different    levels    of    severity.       The    composition    and    spatial 
heterogeneity   of    the    landscapes    generated   by    the   model    under 
these    different    scenarios    were    compared.       An    unsuccessful 
attempt    was    also    made    to    compare    landscapes    generated    by    the 
model   with    independent    data    from   a    1990    SPOT    satellite    image 
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Figure  4-4.   Flow  chart  of  spatial  modeling  and  analysis 


Model  Design 


The  spatial  model  simulated  succession,  fire,  and 
vegetation  recovery  using  a  geographic  information  system  (GIS) 
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A  succession  component  of  the  model  moved  vegetation  pixels  on  a 
map  through  a  sequence  of  serai  stages  by  increasing  their 
designated  successional  class.   The  spreading  of  sawgrass  into 
wet  prairie  via  vegetative  reproduction  also  was  simulated, 
based  on  the  juxtaposition  of  pixels  of  these  two  serai  stages. 
A  burn  component  in  the  model  changed  pixels  from  one 
successional  class  to  another  based  on  their  current  class  and 
their  distance  inward  from  the  edge  of  an  island.   Vegetation 
recovery  was  simulated  when  the  succession  component  of  the 
model  ran  after  a  fire.   Hydrologic  influences  on  both  fire  and 
succession  were  not  explicit  in  the  model.   They  were  held 
constant  to  isolate  the  effects  of  fire.   The  generation  of  new 
islands  by  peat  batteries  or  alligators  was  not  included  in  the 
model,  because  there  were  insufficient  data  documenting  these 
sources  of  island  formation  in  the  Everglades. 
GIS  software 

The  model  was  constructed  using  ERDAS  Inc.'s  VGA  ERDAS 
software,  version  7.5  (ERDAS  1991).   The  main  GIS  programs  used 
to  run  the  model  were  the  SEARCH  (spatial  analysis)  and  GISMO 
(GIS  Modeling)  programs. 

The  ERDAS  SEARCH  program  calculates  the  distance  from  each 
pixel  to  other  pixels  of  a  specified  class  or  classes.   The 
SEARCH  program  was  used  in  each  model  iteration  to  generate  maps 
of  spatial  information  from  the  successional  vegetation  base 
map. 

The  ERDAS  GISMO  program  uses  a  set  of  model  rules  to 
produce  model  output  based  on  conditions  found  in  different 
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maps.   The  user  specifies  logic  statements  or  equations  that 
define  the  model  rules.   The  GISMO  program  considers  each  pixel 
in  a  base  map,  and  queries  other  map  layers  to  find  data  values 
for  the  same  location.   These  values  are  then  run  through  the 
model  rules,  producing  a  new  value  for  the  pixel.   The  results 
are  saved  as  a  new  map. 

The  SEARCH  and  GISMO  programs  were  both  required  to 
produce  a  single  iteration  of  the  model.   Each  model  iteration 
represented  a  10-year  time  step.   When  the  programs  were  run 
repeatedly  in  batch  mode,  the  model  simulated  landscape  change 
over  an  extended  time. 
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Figure  4-5.   Diagram  of  model  structure 
Model  structure 

The  structure  of  the  model  is  shown  in  Figure  4-5.   At  the 
start  of  a  model  iteration  (the  top  of  Figure  4-5),  a  base 
vegetation  map  was  run  through  the  SEARCH  program  to  create  two 
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spatial  information  layers.   The  first  was  a  "borders"  map,  that 
defined  the  boundaries  between  wet  prairie  and  early 
successional  sawgrass.   In  the  borders  map,  all  pixels  in  the 
"slough/sparse  wet  prairie"  and  "wet  prairie"  classes  (classes 
1-14)  that  were  adjacent  to  a  higher  class  (classes  15  and 
higher)  were  coded  as  1.   Wet  prairie  pixels  were  rarely  found 
adjacent  to  shrub  or  tree  classes.  Pixels  in  classes  15  and 
higher  were  coded  as  0.   Pixels  in  classes  1-14  that  were  not 
adjacent  to  a  higher  class  were  coded  as  2. 

If  a  fire  was  set  to  burn  during  the  iteration,  a  second 
spatial  information  layer,  "distance, "  was  also  created  from  the 
base  map.   This  layer  showed  the  distance  inward  from  the  edges 
of  tree  islands.   Slough/wet  prairie,  wet  prairie,  and  sawgrass 
classes  (classes  1-20)  were  coded  as  0.   Pixels  in  the  shrub  and 
tree  classes  (classes  20  and  above)  were  then  re-coded  as  their 
distance  in  pixels  from  class  0.   The  result  was  a  map  of 
distance  inward  from  island  edges.   For  example,  a  pixel  coded 
as  2  was  20  meters  into  the  interior  of  an  island  from  the 
island  edge. 

A  fire  layer  was  also  specified  for  use  by  the  model. 
If  the  iteration  included  a  fire,  fire  was  switched  on  by 
setting  the  fire  layer  to  "burn"  (1).   If  it  did  not  include  a 
fire,  the  fire  layer  was  switched  to  "no  burn"  (0) . 

After  generating  the  spatial  information  maps  and  setting 
the  fire  switch  for  an  iteration,  the  GISMO  program  produced  the 
iteration's  outcome.   The  outcome  was  a  new  vegetation  map  which 
became  the  base  map  for  the  next  iteration.   The  succession 


129 


rules  within  the  GISMO  model  reassigned  the  successional  class 
of  a  pixel  in  the  base  map,  based  upon  its  present  successional 
class  and  its  location  on  the  "borders"  layer.   In  fire  years, 
the  model  rules  reassigned  a  pixel's  successional  class  to  a 
lower  class,  based  upon  its  present  class  and  its  location  on 
the  distance  layer. 
Model  base  map 

The  base  map  used  for  the  initial  conditions  of  the  model, 
referred  to  hereafter  as  the  "starting  base  map,"  was  produced 
from  a  merged  data  set  of  April  1987  SPOT  panchromatic  and 
multispectral  satellite  imagery.   The  processing  of  this  image 
is  described  in  detail  in  Chapter  2.   A  5000  hectare  subset  of 
the  satellite  data  set  (see  Figure  2-2)  was  classified  to 
produce  a  28-class  unsupervised  classification.   The  28  classes 
were  then  identified  as  vegetation  classes  and  ranked  in  order 
of  a  successional  progression  (see  Table  2-6).   Aerial 
photography  and  a  vegetation  map  classified  from  the  same  data 
using  1987  ground  truth  (Richardson  et.  al  1990)  were  used  as  a 
reference  for  the  class  identifications. 

The  staring  base  map  was  classified  into  28  successional 
classes  rather  than  just  the  five  basic  vegetation  types  or 
serai  stages  shown  in  Figure  1-8.   This  was  necessary  because 
the  conversion  of  an  entire  pixel  from  one  serai  stage  to  the 
next  higher  stage  is  not  completed  in  a  single  decade  (the  time 
step  of  the  model) .   Serai  stages  were  subdivided  into  a  number 
of  successional  classes  representing  gradations  between  stages. 
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The  number  of  classes  in  a  serai  stage  depended  on  the  estimated 
duration  of  that  stage  (see  Table  4-1) . 


Table  4-1.   Subdivisions  of  serai  stages  in  model's  starting 
base  map 


Serai  Stage 

Successional 
Classes 

Duration  of  Stage 

Slough/Sparse  wet 
prairie 

1-9 

90  years 

Wet  prairie 

10-14 

40  years 

Sawgrass 

15-19 

50  years 

Brush 

20-23 

40  years 

Trees 

24-28 

no  limit 

Model  rules  and  rates 

The  rates  built  into  the  model  were  estimated  from 
hydrologic,  vegetation  and  fire  data  in  the  ecological 
literature.   Assumptions  regarding  vegetation  succession  and  the 
spatial  distribution  of  vegetation  were  based  on  the  landscape 
information  collected  in  this  work  (see  Chapter  2),  and  on  other 
studies  of  Everglades  vegetation  (Loveless  1959,  Pesnell  and 
Brown  1977,  Olmsted  and  Loope  1984,  Duever  1984,  Kushlan  1990), 
particularly  two  studies  of  vegetation  community  distribution 
and  environmental  variables  in  Water  Conservation  Area  1  (Pope 
1989,  Richardson  et  al.  1990).   Rates  of  vegetation  succession 
were  based  on  studies  of  the  Everglades  and  the  Okefenokee  swamp 
(Craighead  1971,  Deuver  and  Riopelle  1983,  Glasser  1986) . 
Ranges  of  fire  frequency  were  based  on  fire  record  in  Everglades 
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National  Park  (Taylor  1981).   Rules  for  fire  behavior  were  based 
on  observations  of  a  fire  in  Water  Conservation  Area  1  (see 
Chapter  2)  and  on  fire  ecology  literature  (Wade  et  al.  1980, 
Anderson  1982,  Rothermel  1983,  Burgan  and  Rothermel  1984, 
Andrews  1986,  Andrews  and  Chases  1989) . 

The  initial  model  rules,  before  verification  and  fine- 
tuning  of  the  model,  are  shown  in  Table  4-2.   They  were  designed 
to  simulate  the  conceptual  model  shown  in  Figure  4-3.   If  there 
was  not  a  fire  during  an  iteration  (i.e.  burn  =  0),  pixels  in 
each  successional  class  advanced  to  the  next  higher  successional 
class.  Because  sawgrass  was  assumed  to  spread  vegetatively  into 
wet  prairie,  succession  was  accelerated  in  wet  prairie  pixels 
adjacent  to  sawgrass  pixels.   This  was  accomplished  by  advancing 
all  wet  prairie  pixels  defined  as  1  on  the  borders  map  to  the 
lowest  sawgrass  class  (class  15). 


Table  4-2.   Initial  landscape  mo 

del  rules 

Succession  without  fire 

Condition 

Result 

fire  is  0  and  border  is  0 

veg  =  veg  +  1 

fire  is  0  and  border  is  1 

veg  =  15 

fire  is  0  and  border  is  2 

veg  =  veg  +  1 

Succession  with  fire 

Condition 

Result 

fire  is  1  and  distance  >=  3 

veg  =  veg  +  1 

fire  is  1  and  distance  <  3  and 
veg  <  5) 

veg  =  1 

fire  is  1  and  distance  <  3  and 
veg  >=  5 

veg  =  veg  -  4 
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If  fire  was  set  to  burn  during  an  iteration,  succession 
occurred  only  in  core  island  pixels  over  30  meters  from  an 
island  edge  (i.e.  with  a  value  greater  than  three  on  the 
distance  map) .   All  other  pixels  were  reassigned  three  classes 
lower,  to  simulate  the  net  effect  of  succession  during  the 
iteration  and  damage  from  the  fire. 

Model  Verification 

The  initial  model  was  run  to  check  its  performance,  in 
order  to  adjust  the  rules  or  rates  so  that  it  reproduced  the 
assumptions  that  went  into  its  design. 


Tests  must  be  made  to  ensure  that  the  model  reproduces 
accurately  the  data  that  were  used  in  its  construction. 
This  process  of  verification,  as  it  is  sometimes  called, 
is  tautologous  in  principle  but  nevertheless  a  necessary 
check  that  the  mechanisms  of  the  model  are  in  fact  doing 
what  the  modeller  thinks  they  are  doing.  (Kitching 
1983:42) 


The  model  was  first  run  for  20  iterations,  covering  a  200- 
year  period  of  time,  without  fire.   The  movement  of  pixels  from 
class  to  class  each  iteration  was  carefully  tracked  to  verify 
that  the  model  rules  were  exhaustive  and  exclusive,  and  were 
functioning  as  planned. 

To  facilitate  the  analysis  of  the  model's  output,  the 
results  of  each  iteration  were  summarized  by  serai  stage. 
Pixels  in  the  28  successional  classes  in  each  output  base  map 
were  grouped  into  the  five  serai  stages  using  the  ERDAS  RECODE 
program.   All  statistics  presented  for  the  model  simulations  in 
the  remainder  of  this  chapter  are  calculated  from  base  maps 
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grouped  into  these  same  serai  stages.   For  each  iteration  in  the 
initial  model  run  without  fire,  the  percentage  of  the  base  map 
in  each  serai  stage  was  calculated  using  the  ERDAS  BSTATS 
program  and  plotted  (Figure  4-6). 
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Figure  4-6.   Composition  of  the  base  map  during  the  initial  run 
of  the  spatial  model  (percentage  of  map  in  each  serai  stage) 


As  succession  in  the  model  progressed  during  the 
simulation,  the  dominant  type  in  the  landscape  changed  from  wet 
prairie  to  sawgrass,  then  to  shrubs,  and  finally  to  trees.   By 
the  160-year  iteration,  trees  occupied  100%  of  the  base  map. 

The  progression  of  succession  in  the  initial  model  was 
unrealistically  rapid,  particularly  the  change  from  slough  and 
wet  prairie  to  sawgrass.   The  landscape  became  less 
heterogeneous  as  much  of  the  area  was  converted  to  a  single 
type,  which  progressed  through  the  successional  series.   The 
source  of  this  problem  was  identified  as  the  acceleration  of 
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succession  of  both  wet  prairie  pixels  and  slough/sparse  wet 
prairie  pixels  next  to  sawgrass.   The  acceleration  of  slough/wet 
prairie  succession  was  not  realistic,  since  sawgrass  is  not 
likely  to  spread  into  slough.   In  the  improved  model,  this  rule 
was  changed  to  apply  to  wet  prairie  pixels  only. 

The  initial  model  was  next  run  with  fires  simulated  at 
different  intervals,  to  test  the  burn  rules.   Five  model 
simulations  of  190  years  each  (19  iterations)  were  run  with  fire 
intervals  of  30  to  70  years.   A  burn  schedule  for  the 
simulations  is  shown  in  Table  4-3.   The  landscape  composition  of 
the  base  maps  at  the  end  of  each  simulation  is  shown  in  Figure 
4-7.   The  landscape  composition  was  measured  by  the  percentage 
of  the  base  map  in  each  serai  stage,  with  the  slough/sparse  wet 
prairie  and  wet  prairie  stages  grouped  for  plot  simplicity. 

Some  trends  in  the  dominance  of  the  serai  stages  at  the 
different  fire  intervals  had  been  expected.   Slough  and  wet 
prairie  were  expected  to  dominate  the  landscape  at  the  shorter 
fire  intervals.   Trees  were  expected  to  dominate  in  the  longer 
intervals . 

The  expected  trends  in  serai  stage  dominance  were  present 
in  the  simulation  results  in  Figure  4-7.   However,  the 
differences  between  the  landscape  compositions  of  the  five  fire- 
interval  simulations  were  greater  than  expected.   The  difference 
in  wet  prairie  coverage  between  the  30-year  and  40-year 
intervals,  for  example,  was  unreasonable.   This  was  partly 
attributable  to  the  accelerated  succession  of  slough/sparse  wet 
prairie  to  sawgrass,  as   mentioned  above.   Another  reason  was 
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Figure  4-7.   Initial  landscape  model  using  different  fire 
intervals:  percentage  of  base  map  in  each  serai  stage  after  190 
years 


that  vegetation  was  in  a  different  state  of  fire  recovery  at  the 
end  of  the  different  the  simulation  runs.   For  example,  when  the 
190-year  simulation  with  a  30-year  fire  interval  was  over,  the 
base  map  had  been  burned  during  the  previous  iteration.   After 
190  years  of  the  40-year  fire  interval  simulation,  a  fire  had 
not  burned  the  base  map  for  three  iterations.   In  the  first 
case,  vegetation  was  in  early  recovery,  and  in  the  second  case 
recovery  was  well  advanced. 

A  decision  was  made  to  address  this  problem  by  comparing 
the  latest  iteration  in  each  simulation  that  was  two  years  post- 
fire,  rather  than  comparing  the  last  iteration  in  each 
simulation.   These  selected  iterations  are  compared  in  the 
remainder  of  this  work,  and  are  marked  in  Table  4-3.   When  the 
selected  iterations  were  compared,  the  differences  between  the 
simulations  were  more  realistic  (see  Figure  4-8).   The 
distributions  of  the  serai  stages  across  the  range  of  fire 
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intervals  followed  more  gradual  curves.   However,  none  of  the 
simulations  produced  a  landscape  with  vegetation  as  evenly 
distributed  as  the  starting  base  map. 

The  landscape  spatial  heterogeneity  produced  by  the 
different  simulations  was  measured  by  the  number  of  patches  of 
vegetation  in  the  base  maps.   The  ERDAS  CLUMP  program  was  used 
to  identify  vegetation  patches  in  the  base  maps  after  they  had 
been  re-coded  from  the  28  successional  classes  to  the  five  serai 
stages.   The  patches  were  counted  using  the  ERDAS  BSTATS 
program,  and  are  presented  in  Figure  4-9.   The  sizes  of  the 
largest  patch  in  each  base  map  also  were  noted,  and  are 
presented  in  Figure  4-9  as  a  percentage  of  the  base  map  area. 
The  greater  the  size  of  the  largest  patch,  the  greater  the 
connectivity  of  the  dominant  serai  stage  in  the  base  map. 

Landscape  heterogeneity  was  lowest  when  fires  burned  most 
frequently,  increasing  the  amount  of  wet  prairie,  and  where 
fires  burned  least  frequently,  allowing  tree  and  brush  patches 
to  coalesce  (see  Figures  4-8  and  4-9) .   This  relationship 
between  disturbance  frequency  and  landscape  heterogeneity  has 
been  noted  in  landscape  ecology  literature  (Turner  1987)  . 
However,  in  all  the  simulations,  the  numbers  of  patches  were 
much  lower  and  the  sizes  of  the  largest  patches  were  much 
greater  than  in  the  starting  base  map.   Thus  the  model  was 
simplifying  the  landscape.   Some  simplification  was  inevitable, 
because  the  model  was  an  abstraction  of  the  ecosystem,  and 
lacked  a  mechanism  besides  succession  to  generate  of  new  patches 
(such  as  peat  batteries  or  alligators) .   However,  it  was  hoped 
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that    removing    slough-to-sawgrass    succession    and    improving    the 
model' s    fire    component   would   allow    the   model    base   maps    to    retain 
more    spatial    heterogeneity. 
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Figure  4-8.   Initial  model  results  using  different  fire 
intervals:  percentage  of  base  map  in  each  serai  stage  at 
selected  iterations 
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Improved  Model  Simulations 

Methods 

The  GIS  model  rules  were  revised  after  verification  of  the 
model,  and  are  shown  in  Table  4-4.  As  in  the  initial  model,  if 
there  was  not  a  fire  during  an  iteration,  pixels  in  each  class 
advanced  to  the  next  higher  successional  class.  Succession  was 
still  accelerated  in  wet  prairie  pixels  adjacent  to  sawgrass 
pixels,  but  slough/sparse  wet  prairie  pixels  were  excluded  from 
the  rule.  It  was  assumed  that  the  water  depth  in  slough  pixels 
was  too  low  for  sawgrass  colonization. 

To  improve  the  performance  of  the  fire  component  of  the 
model,  new  fire  rules  were  written  for  two  levels  of  fire 
severity.   Using  the  new  rules,  the  improved  model  could  be  run 
to  simulate  either  a  moderate  fire  without  tree  and  shrub 
mortality,  or  a  severe  fire  with  tree  and  shrub  mortality. 

The  improved  model  was  first  run  without  fire.   As  with 
the  initial  model,  dominance  of  the  landscape  progressed  from 
slough  and  wet  prairie  to  sawgrass,  then  to  brush,  then  trees 
(see  Figure  4-10) .   However,  the  landscape  composition  was  more 
even  throughout  the  simulation,  and  thus  more  realistic  than  the 
initial  model.   Trees  occupied  all  of  the  base  map  after  190 
years,  rather  than  160  years. 

Two  separate  sets  of  model  rules  were  written  for  fire  in 
the  improved  model.   The  first  set  simulated  a  moderate  fire 
without  shrub  and  tree  mortality.   If  there  was  a  fire  in  an 
iteration,  succession  occurred  only  in  core  island  pixels  over 
30  meters  from  an  island  edge  (i.e.  with  a  value  greater  than 
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Table  4-4 .   Improved  landscape  model  rules 

Succession  without  fire 

Condition 

Result 

fire  is  0  and  border  is  0 

veg  =  veg  +  1 

fire  is  0  and  border  is  1  and 
veg  <  10 

veg  =  veg  +  1 

fire  is  0  and  border  is  1  and 
veg  >=  10 

veg  =  veg  +  3 

fire  is  0  and  border  is  2 

veg  =  veg  +  1 

Succession  with  fire:  moderate 

Condition 

Result 

fire  is  1  and  distance  >=  3 

veg  =  veg  +  1 

fire  is  1  and  distance  <  3  and 
veg  <  4) 

veg  =  1 

fire  is  1  and  distance  <  3  and 
veg  >=  4  and  veg  <  2  0 

veg  =  veg  -  3 

fire  is  1  and  distance  <  3  and 
veg  >=  20  and  veg  <  24 

veg  =20 

fire  is  1  and  distance  <  3  and 
veg  >=  24) 

veg  =  2  2 

Succession  with  fire:  severe 

Condition 

Result 

fire  is  1  and  distance  >=  3 

veg  =  veg  +  1 

fire  is  1  and  distance  <  3  and 
veg  <  5) 

veg  =  1 

fire  is  1  and  distance  <  3  and 
veg  >=  5  and  veg  <  2  4 

veg  =  veg  -  4 

fire  is  1  and  distance  <  3  and 
veg  >=  24) 

veg  =  2  2 

141 


— e — 

-Slough/ sparse 
wet  prairie 

100  -, 
90  - 

-Wet   praire 

80  - 

70  ■ 

-Sawgrass 

60  • 
50  ■ 

— * — 

-Shrubs 

40  - 
30  - 



Trees 

20  - 

10  - 

0  ■ 

80  100  120  140  160  180 

Iteration 


Figure  4-10.   Composition  of  the  base  map  during  a  run  of  the 
improved  landscape  model  (percentage  of  map  in  each  serai  stage) 


three  on  the  distance  map) .   Slough/sparse  wet  prairie  pixels, 
wet  prairie  pixels,  and  sawgrass  pixels  (classes  1-19)  were 
decreased  three  classes  to  simulate  the  net  effect  of  succession 
during  the  iteration  followed  by  damage  from  fire.   Shrub  pixels 
were  set  back  to  the  lowest  shrub  class  (class  20)  to  simulate 
the  net  effect  of  succession  followed  by  fire  damage.   Tree 
pixels  were  set  back  to  the  middle  shrub  class  (class  22)  to 
simulate  the  effect  of  succession  followed  by  top-killing  of  the 
tree  by  fire.   It  was  assumed  that  trees  would  stump-sprout  and 
reach  tree  stature  after  20  years  (two  fire-free  iterations). 

The  second  set  of  fire  rules  simulated  a  severe  fire, 
allowing  for  shrub  and  tree  mortality,  but  without  peat  burning. 
If  there  was  a  fire  in  an  iteration,  succession  occurred  only  in 
core  island  pixels  over  30  meters  from  an  island  edge  (i.e.  with 
a  value  greater  than  three  on  the  distance  map) .   Slough/sparse 
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wet  prairie  pixels,  wet  prairie  pixels,  sawgrass  pixels,  and 
shrub  pixels  (classes  1-23)  were  decreased  four  classes  to 
simulate  the  net  effect  of  succession  during  the  iteration 
followed  by  severe  damage  from  fire.   Tree  pixels  were  set  back 
to  the  middle  shrub  class  (class  22) .   The  effect  of  this  last 
rule  was  that  trees  died  only  if  the  fire  interval  was  less  than 
thirty  years . 

As  with  the  initial  model,  the  improved  model  was  run  with 
fires  simulated  at  different  intervals.   A  20-year  interval,  not 
included  in  the  initial  model  verification,  was  added. 
Simulations  of  190  years  each  (19  iterations)  were  run  with  fire 
intervals  of  20,  30,  40,  50,  60  and  70  years.    One  set  of  fire 
interval  simulations  was  run  using  the  moderate  fire  rules,  and 
one  set  was  run  using  the  severe  fire  rules . 

The  results  of  the  simulations  were  evaluated  by  again 
comparing  the  latest  iteration  in  each  simulation  that  was  two 
years  post-fire,  with  the  exception  of  simulations  using  a  20- 
year  fire  interval  (see  Table  4-3) . 
Results 

The  compositions  of  the  base  maps  at  the  end  of  the  six 
simulations  using  the  moderate  fire  rules  are  shown  in  Figure  4- 
11.   The  landscapes  produced  by  the  improved  model  were  more 
evenly  composed  and  more  heterogeneous  than  those  produced  by 
the  initial  model  (see  Figures  4-11  and  4-13).   The  20-year  fire 
interval  simulation  most  retained  the  qualities  of  the  starting 
base  map  made  from  the  1987  satellite  image  of  Water 
Conservation  Area  1.   The  20-year  fire  interval  produced  a 
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Figure  4-11.   Improved  model  results  using  moderate  fire 
severity  and  different  fire  intervals:  percentage  of  base  map  in 
each  serai  stage  at  selected  iterations 


landscape  with  slightly  more  wet  prairie  and  shrubs,  and  fewer 
trees  and  sawgrass,  than  the  starting  base  map.  The  landscape 
contained  nearly  as  many  small  patches  as  the  starting  base  map. 

The  30-year  and  40-year  interval  simulations  resulted  in 
the  most  evenly  composed  landscapes  with  similar  patch  numbers. 
Sawgrass  was  dominant  in  both  simulations,  but  the  30-year 
interval  simulation  contained  more  wet  prairie  and  less  sawgrass 
than  the  40-year  simulation. 

In  the  fire  intervals  of  50  years  and  longer,  the 
percentage  of  the  base  map  occupied  by  trees  increased 
dramatically.   The  longer  fire-free  period  allowed  sawgrass  to 
succeed  to  brush,  yet  the  moderate  fire  rules  did  not  kill  brush 
and  convert  it  back  to  sawgrass.   The  longer  fire-free  period 
allowed  the  increasing  brush  area  to  succeed  to  trees.   While 
the  50-year  and  60-year  interval  base  maps  had  less  sawgrass  and 
more  trees  than  the  initial  model  produced  at  these  intervals, 
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they  also  retained  some  small  patches  of  slough  that  preserved 
more  landscape  heterogeneity. 

The  compositions  of  the  base  maps  in  each  serai  stage  at 
the  end  of  the  simulations  using  the  severe  fire  rules  are  shown 
in  Figure  4-12.   All  the  severe  fire  simulations  produced  less 
evenly  composed  landscapes  than  their  moderate  counterparts. 
Brush  mortality  specified  by  the  severe  fire  rules  reduced  or 
eliminated  brush  at  all  intervals.   Trees  were  not  eliminated, 
because  the  fire  rules  gave  them  more  resilience  than  brush,  and 
island  cores  were  protected  from  fire. 
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Figure  4-12.   Improved  model  results  using  high  fire  severity 
and  different  fire  intervals:  percentage  of  base  map  in  each 
serai  stage  at  selected  iterations 


Severe  fires  at  the  20-year  interval  converted  the 
landscape  almost  entirely  to  wet  prairie.   Only  124  small 
patches  of  trees  and  sawgrass  remained  after  170  years.   Unlike 
the  moderate  fire  simulation,  most  of  the  30-year  fire  interval 
was  occupied  by  wet  prairie.   The  40-year  interval  simulation 
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was  dominated  by  sawgrass,  as  in  the  moderate  simulation.   In 
the  severe  simulations  with  fire  intervals  greater  than  40, 
sawgrass  was  more  prevalent  than  in  the  moderate  simulations, 
owing  to  the  mortality  of  brush.   Areas  that  were  brush  and 
trees  in  the  moderate  50-year  and  60-year  simulations  were 
converted  to  sawgrass  in  the  severe  simulations. 
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Figure  4-13.  Improved  model  results  using  moderate  and  severe 
fires  and  different  fire  intervals:  landscape  heterogeneity  at 
selected  iterations 

Because  the  severe  fire  simulations  changed  the  balance 

between  the  woody  and  non-woody  serai  stages,  they  produced 

different  levels  of  heterogeneity  than  their  moderate 

counterparts.   Patch  density  in  the  moderate  simulations  was 

highest  in  the  20-year  interval,  because  this  simulation  changed 

the  starting  base  map  the  least.   The  frequent  fires  continually 

returned  pixels  that  had  raised  their  successional  class  back  to 

their  previous  states.   Yet  the  fires  did  not  damage  vegetation 

enough  to  eliminate  many  patches. 
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Some  heterogeneity  was  lost  in  the  moderate  simulations 
with  fire  intervals  longer  than  20  years,  because  much  of  the 
heterogeneity  of  the  starting  base  map  lay  in  the  interspersion 
of  slough/sparse  wet  prairie  and  wet  prairie  vegetation.   These 
types  decreased  with  longer  fire  intervals,  so  their 
contribution  to  heterogeneity  also  decreased.   At  the  same  time, 
longer  fire  intervals  allowed  succession  to  progress  further 
between  fires.   New  patches  of  brush  and  trees  appeared,  created 
by  the  succession  of  sawgrass  and  brush.   Heterogeneity 
increased  from  the  30-year  to  the  50-year  interval.   By  the  60- 
year  interval,  so  many  patches  of  brush  and  trees  had  appeared 
that  they  were  consolidating  into  larger  patches,  reducing 
heterogeneity. 

In  the  severe  fire  simulations,  these  patterns  of  change 
were  altered  by  the  higher  level  of  fire  damage.   The  greater 
damage  reduced  all  wet  prairie  pixels  to  slough/sparse  wet 
prairie  classes.   As  a  result,  no  wet  prairie  heterogeneity  was 
retained  in  any  of  the  simulations.   The  number  of  brush  patches 
was  reduced  at  all  fire  intervals,  further  eroding 
heterogeneity.   The  40-year  interval's  heterogeneity  was 
particularly  low,  because  sawgrass  was  so  prevalent  that 
sawgrass  patches  coalesced. 

The  severe  fire  rules  slowed  the  spread  of  trees.   The 
merging  of  tree  patches  seen  in  the  moderate  simulations  with 
intervals  of  50  years  and  more  was  delayed  in  the  severe 
simulations.   The  simulations  with  severe  fires  at  these 
intervals  were  more  heterogeneous  than  the  moderate  simulations. 
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Figure  4-14.   Tree  islands  in  landscapes  simulated  by  the 
spatial  model 
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Figure  4-14--continued 
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Tree  islands  in  the  spatial  model 

Because  tree  islands  are  the  focus  this  work,  tree  islands 
in  the  model  simulations  were  given  additional  study.   The  base 
maps  presented  in  the  previous  section  under  the  moderate  and 
severe  fire  rules  were  again  re-coded  using  the  ERDAS  RECODE 
program.   This  time  everything  but  tree  pixels  were  re-coded  to 
zero.   These  maps  are  shown  in  Figure  4-14,  with  tree  islands 
appearing  in  black. 

Figure  4-14  illustrates  the  changing  patterns  of  tree 
islands  in  the  landscape  under  different  fire  regimes.   As  fires 
became  less  frequent,  succession  advanced  until  tree  islands 
merged  together.   As  fires  became  more  frequent  or  severe,  small 
islands  were  reduced  in  size  and  disappeared.   The  largest 
islands  in  the  landscape  persisted  and  maintained  an  elongate 
shape  even  at  high  fire  frequency  and  severity. 

The  patches  of  trees  produced  by  the  different  simulations 
were  identified  and  counted  using  the  ERDAS  CLUMP  and  BSTATS 
programs.   These  patch  data  are  presented  in  Figure  4-15.   The 
loss  of  tree  patches  at  the  shortest  fire  interval  is  apparent 
in  Figure  4-15,  as  is  the  drop  in  the  number  of  patches  from 
island  consolidation  at  the  longer  intervals.   The  Figure 
suggests  that  a  fire  interval  between  20  and  30  years  would 
produce  a  density  of  tree  islands  close  to  the  starting  base 
map. 

While  most  vegetation,  including  small  islands,  was 
converted  to  slough/sparse  wet  prairie  by  severe  fires  in  the 
20-year  interval  simulation,  large  tree  islands  persisted.   To 
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verify  that  these  large  islands  remained  because  the  cores  of 
large  islands  were  protected  from  fire,  the  model  was  run  with  a 
20-year  fire  interval  and  severe  fire  without  the  core 
protection  rule.   At  the  end  of  the  simulation  without  the 
protection  rule,  only  slough/sparse  wet  prairie  was  present  (see 
Tables  4-5  and  4-6) .   The  protection  of  island  cores  was  the 
mechanism  keeping  large  islands  in  the  landscape  at  this  fire 
frequency. 
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Figure  4-15.  Improved  model  results  using  moderate  and  severe 
fires  and  different  fire  intervals:  number  of  tree  islands  at 
selected  iterations 


Table  4-5.   Landscape  composition  after  20-year  fire  interval 
simulations,  (percentage  of  area  by  serai  stage) 


Fire    severity 

Slough 

Wet    pr. 

Serai 
stage 

Sawgrass 

Shrubs 

Trees 

Moderate 

69.83 

3.  18 

13.22 

10.7 

3.07 

High 

97  .  93 

0.  49 

0.  62 

0.27 

0.68 

High   without 
core    protection 

100 

0 

0 

0 

0 
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Table  4-6.   Number  of  patches  after  20-year  fire  interval 
simulations 


Fire  severity 

Patches 

Moderate 

18368 

High 

125 

High  without 
core  protection 

1 

Comparison  of  Model  with  Independent  Data 


Model  validation 

The  validation  of  a  model  is  achieved  by  comparing  the 
output  of  the  model  with  independent  data  not  used  in  the 
development  of  the  model.   A  model  should  be  used  as  a 
predictive  tool  only  if  it  is  validated. 

Models  are  validated  by  running  them  on  a  different  site 
or  in  a  different  time  period  from  which  they  were  developed. 
Assume  for  example  that  a  forest  stand  model  was  developed  and 
verified  using  data  from  some  time  period  at  a  certain  site.   To 
validate  the  model,  it  could  be  run  over  the  same  time  period 
using  input  variables  from  a  different  site.   The  output  would 
then  be  compared  with  independent  stand  data  collected  at  the 
new  site.   Alternatively,  a  hydrologic  model  developed  and 
verified  using  data  from  some  time  period  at  a  certain  site 
could  be  validated  by  using  climatic  variables  from  a  different 
period  of  time  to  run  the  model,  and  comparing  the  output  with 
water  level  recordings  from  the  site. 

The  GIS  model  described  in  this  work  could  not  be 
validated  by  running  it  in  a  different  time  period.   The  rate  of 
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succession  and  the  model' s  time-step  are  so  long  compared  to  the 
period  of  time  during  which  landscape  data  have  been  recorded 
that  this  was  impossible. 

Although  not  an  optimal  strategy,  an  area  of  Water 
Conservation  Area  1  separate  from  the  5000  hectare  modeled  area 
could  be  used  to  run  the  model  for  validation.   Data  from  all  of 
Water  Conservation  Area  1  was  used  in  the  development  of  the 
conceptual  model,  but  only  the  5000-hectare  portion  of  the  area 
was  used  to  verify  the  model. 

Data  used  for  validation  should  be  independent  from  the 
data  used  in  model  development.   An  independently-produced  map 
with  information  comparable  to  the  1987  model  base  map,  or 
statistics  derived  from  such  a  map,  is  required  for  validation 
of  the  spatial  model.   Here  the  time  scale  of  the  model  again 
poses  a  problem.   To  validate  just  one  time-step  of  the  model 
beginning  with  the  1987  classified  satellite  image  would  require 
map  data  from  1997.   An  earlier  map  could  be  used  as  the  model 
base  map,  run  until  1987,  and  compared  with  the  1987  classified 
image,  but  this  would  require  a  map  from  approximately  1947, 
1957,  1967  or  1977. 

There  are  no  existing  maps  of  Water  Conservation  Area  1 
that  fulfill  these  validation  requirements.   Existing  vegetation 
maps  are  too  general  in  their  classes,  dividing  the  interior  of 
the  area  into  one  or  two  types.   Pre-1982  satellite  data  are  too 
coarse  in  their  resolution  to  be  compared  with  the  model. 
Aerial  photographs  from  as  far  back  as  1948  are  available, 
although  the  difficulties  in  interpreting  and  registering  old 
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photography  in  Water  Conservation  Area  1  (as  discussed  at  the 
end  of  Chapter  2)  are  great.   Because  the  interpretation  of 
aerial  photographs  is  subjective,  the  builder  of  the  model  could 
not  avoid  bias  in  interpreting  a  new  vegetation  map.   Therefore 
vegetation  would  have  to  be  specially  mapped  for  validation 
purposes  by  an  independent  source.   No  funds  were  available  for 
this  purpose. 
Comparison  of  model  and  satellite  data 

Because  validation  of  the  GIS  model  of  Water  Conservation 
Area  1  was  not  possible,  an  attempt  was  made  to  compare 
independent  data  with  the  burn  function  of  the  model.   The  area 
chosen  for  the  comparison  was  a  subset  of  Water  Conservation 
Area  1  outside  the  5000-hectare  area  used  for  model 
verification,  and  within  the  area  burned  by  the  1989  fire. 
Unclassified  satellite  data  were  chosen  for  comparison  with  the 
model,  to  avoid  the  bias  associated  with  data  classification. 
Unclassified  1987  satellite  data  from  the  subset  area  was  to  be 
compared  with  unclassified  1990  satellite  data  to  detect  burn 
damage  from  the  1989  fire.   The  model  was  to  be  run  in  the 
subset  area  for  a  single  iteration  with  fire,  using  the  1987 
base  map.   The  simulated  burn  damage  was  to  be  compared  with  the 
real  burn  damage  detected  by  the  raw  satellite  data. 

Unfortunately,  the  change  detection  procedure  run  on  the 
1987  and  1990  satellite  images  was  unsuccessful  in  showing  the 
effects  of  the  1989  fire,  and  could  not  be  compared  with  burn 
damage  simulated  by  the  model.   The  difference  in  water  levels 
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in  Water  Conservation  Area  1  between  the  two  raw  satellite 
images  caused  the  failure  of  the  burn  damage  detection. 

Soil  moisture  and  standing  water  strongly  influence  the 
spectral  reflectance  measured  by  the  SPOT  satellite.   Rainfall 
in  1987  was  slightly  above  normal,  while  1990  was  a  severe 
drought  year.   The  change  in  water  level  was  greatest  between 
the  two  dates  in  the  area  burned  by  the  fire,  since  fire  spread 
throughout  the  driest  areas  in  1989.   The  spectral  difference 
between  the  1987  image  and  1990  image  due  to  water  level  was 
much  greater  than  the  spectral  difference  between  the  two  images 
due  to  fire.   Although  the  extent  of  the  fire  was  clearly 
apparent  on  the  1990  image,  and  there  were  measurable 
differences  within  the  image  in  reflectance  between  damaged  and 
undamaged  areas,  these  differences  were  overwhelmed  by  the  much 
greater  variation  in  water-derived  reflectance  between  the  two 
images . 

The  first  change-detection  technique  tried  on  the  two 
satellite  images  was  simple.   With  this  method,  a  band  from  one 
image  is  displayed  on  the  computer  monitor  in  one  color  with  a 
band  from  the  other  image  in  another  color.   The  pixels  that 
differ  between  the  two  images  appear  in  different  color  from 
those  that  do  not  differ. 

Pairs  of  bands  from  each  image  were  displayed,  using  all 
four  bands  (bands  1,  2,  3,  and  the  NDVI  band) .   However,  none  of 
the  band  comparisons  showed  the  effects  of  fire. 

A  second  technique  for  change  detection  was  also  used  on 
the  images.   This  method  subtracts  the  pixel  values  of  one  image 


156 


from  the  pixel  values  of  the  other  image.   The  differences  are 
kept  as  positive  numbers  by  adding  a  constant  to  them. 
Differences  are  calculated  for  each  band  separately.   The 
greater  the  difference,  the  greater  the  change. 

When  an  environmental  condition  that  affects  the  image 
evenly,  such  as  haze,  varies  between  two  images,  the  difference 
it  introduces  can  be  corrected  for.   Rather  than  calculate  the 
difference  between  the  band  values  for  each  image,  the  values  in 
each  band  are  first  re-coded  as  the  difference  from  the  mean  for 
that  band.   The  differences  from  the  mean  are  then  used  in  the 
difference  calculation  (ERDAS  1991). 

Using  the  ERDAS  BSTATS,  RECODE,  and  ALGEBRA  programs 
(ERDAS  1991),  the  differences  between  the  two  images  were 
calculated  using  the  differences  from  the  mean  for  each  image. 
However,  the  fire  damage  was  not  highlighted  by  this  technique 
either.   The  correction  using  the  image  mean  was  not  effective, 
because  unlike  haze  or  sun-angle,  water  level  did  not  change 
evenly  throughout  the  image . 

Discussion 

The  design  of  this  spatial  model  of  Water  Conservation 
Area  1  was  innovative,  because  it  used  existing  GIS  software  to 
iteratively  model  landscape  change.   In  creating  the  model,  it 
was  necessary  to  simplify  or  ignore  some  processes  in  the 
ecosystem.   The  model  could  not  be  checked  against  independent 
data.   All  available  data  from  Water  Conservation  Area  1  had  to 
be  used  in  model  development,  because  data  were  so  scarce.   This 
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left  no  independent  data  for  validation.   For  all  of  the  above 
reasons,  the  spatial  model  presented  in  this  chapter  is  not 
meant  to  be,  nor  should  it  be,  used  as  a  predictive  management 
tool . 

In  spite  of  its  limitations,  the  model  served  to  clarify 
possible  interactions  between  fire,  succession  and  landscape 
pattern  in  Water  Conservation  Area  1 .   The  model  was  a 
quantitative  and  visual  expression  of  a  theory.   The  maps  it 
produced  made  visualization  of  the  theory  easier  than  imagining 
it  in  abstract.   The  process  of  creating  the  model  also  required 
a  thorough  review  of  all  the  factors  that  might  be  influencing 
the  system,  even  though  most  of  these  factors  were  not  included 
in  it.   The  model  is  best  suited  to  generate  new  ideas  and 
questions  to  be  investigated. 

The  simulated  landscapes  produced  by  the  GIS  model  showed 
that  both  fire  frequency  and  severity  strongly  influenced  the 
simulated  landscape  pattern.   The  simulations  using  a  moderate 
fire  severity  produced  landscapes  more  similar  to  the  landscape 
pattern  we  currently  observe  in  Water  Conservation  Area  1.   Of 
the  fire  intervals  simulated  with  moderate  fire,  the  optimum 
frequency  for  maintaining  a  vegetation  pattern  close  to  the 
current  pattern  of  Water  Conservation  Area  1  was  20  years.  A 
frequency  between  20  and  30  years  would  come  even  closer.   This 
suggests  that  an  optimal  frequency  range  of  moderate  fires  may 
exist  for  maintaining  the  historical  pattern  of  the  real 
ecosystem. 
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No  simulation  using  severe  fire  produced  a  vegetation 
pattern  close  to  the  present  pattern  in  Water  Conservation  Area 
1 .   The  entire  range  of  severe  fire  intervals  reduced  the 
evenness  of  the  landscape  composition,  because  severe  fire 
favored  sawgrass  over  brush  and  reduced  the  spatial  variation  in 
the  successional  state  of  wet  prairie. 

Fire  controlled  the  pattern  of  tree  islands  produced  by 
the  model  simulations.   Without  fire,  tree  islands  in  the  model 
grew  until  they  consolidated  into  a  single  patch.   Fires  at 
intervals  of  40  to  60  years  resulted  in  the  greatest  numbers  of 
tree  islands,  far  more  than  are  found  in  the  current  pattern  of 
Water  Conservation  Area  1.   Under  the  strongest  fire  regime 
(severe  fires  at  a  20-year  interval),  the  protection  of  island 
cores  built  into  the  model  prevented  the  largest  tree  islands  in 
the  landscape  from  disappearing. 


CHAPTER  5 
SUMMARY  AND  CONCLUSIONS 


As  the  twentieth  century  draws  to  a  close,  the  resources  of 
the  planet  are  becoming  ever  scarcer.   Awareness  of  human 
dependence  upon  natural  systems  is  growing,  as  the  benefits 
provided  by  ecosystems  disappear.   The  management  of  ecosystems 
only  for  species  or  functions  that  benefit  humans,  rather  than  for 
the  maintenance  of  the  natural  system,  has  led  to  the  collapse  of 
some  ecosystems.   In  response  to  a  growing  urgency  for 
conservation,  the  1990 's  began  with  a  new  resource  management 
paradigm  in  the  United  States:  ecosystem  management. 
Implementation  of  ecosystem  management,  however,  has  been  more 
difficult  than  the  decision  to  pursue  that  approach. 

Methods  of  evaluating  whole  ecosystems  or  landscapes  are 
still  being  developed  and  tested.   The  computer  and  remote  sensing 
technology  that  made  the  inventory  of  ecosystems  possible  are 
recent  developments.   This  technology  continues  to  improve  at  a 
rapid  rate.   Many  of  the  technological  limitations  to  the  work 
described  in  this  dissertation  (barriers  of  cost,  computer  disk 
space,  and  software)  have  been  resolved  during  the  completion  of 
this  document.   More  geographic  data  can  now  be  processed  faster, 
on  less  expensive  computers,  to  produce  a  wide  variety  of  spatial 
statistics.   However,  the  body  of  theoretical  and  methodological 
knowledge  to  correctly  apply  these  statistics  to  questions  of 
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and  methodological  knowledge  to  correctly  apply  these  statistics 
to  questions  of  ecosystem  management  is  growing  at  a  slower 
rate.   The  academic  field  of  landscape  ecology  is  still  in  its 
infancy. 

This  doctoral  research  was  an  effort  to  study  an  ecosystem 
using  many  approaches  and  tools,  including  remote  sensing  and 
computing  technology  available  between  1988  and  1993.   The 
ecosystem  under  study,  South  Florida's  Water  Conservation  Area 
1,  was  not  well-studied,  well-traveled,  nor  well-understood.   As 
part  of  the  Everglades,  it  was  the  object  of  both  human-induced 
environmental  change,  and  human  pressure  for  restoration  and 
conservation.   The  existing  pattern  of  its  landscape  was  a 
critical  aspect  of  wildlife  habitat,  particularly  for  endangered 
wading  birds.   The  ecosystem  exemplified  both  the  great  need  for 
and  the  great  difficulties  facing  an  ecosystem  or  landscape 
approach  to  management. 

The  work  described  in  this  dissertation  focused  on  the 
pattern  of  the  landscape  of  Water  Conservation  Area  1,  rather 
than  on  the  plant  and  animal  species  in  it.   The  pattern  of  the 
landscape  is  the  pattern  of  habitats  upon  which  both  plants  and 
animals  are  dependent.   This  research  took  a  holistic  approach 
to  gaining  knowledge  of  the  ecosystem  at  the  landscape  scale. 
The  investigation  of  how  this  ecosystem  produced  the  unique 
pattern  found  in  its  landscape  required  piecing  together  data 
from  many  parts  of  the  system,  at  different  scales,  gathered  by 
many  people. 
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Descriptions  of  vegetation  communities  and  vegetation 
studies  by  other  researchers  were  used  to  define  the  vegetation 
types  making  up  the  landscape.   Studies  of  the  distribution  of 
plant  species  with  respect  to  hydrology,  and  theories  of 
vegetation  succession  and  ecosystem  evolution,  led  to  the 
description  of  a  successional  sequence  in  Water  Conservation 
Area  1.   Studies  of  peat,  including  peat  formation,  oxidation, 
burning,  erosion  and  deposition,  led  to  the  concept  of  not  only 
vegetation  in  flux,  but  the  topography  itself,  changing  more 
rapidly  than  other  geologic  processes.   Observations  and 
theories  of  tree  island  formation  suggested  several  mechanisms 
responsible  for  the  development  of  raised  areas  of  peat, 
although  no  evidence  was  conclusive. 

All  of  these  pieces  of  information  together  led  to  the 
conclusion  that  islands  grow  in  area  over  time.   If  succession 
is  operating,  if  peat  is  being  formed,  if  floating  peat  and 
vegetation  accumulate,  then  islands  must  be  growing,  unless  some 
other  process  limits  growth. 

The  landscape  pattern  of  Water  Conservation  Area  1 
observed  on  the  ground  reinforced  the  ideas  of  vegetation 
succession  and  island  growth.   Core  forested  areas  of  tree 
islands  graded  into  brush  and  sawgrass  at  island  edges,  and 
island  edges  were  fringed  with  sawgrass.   Patches  of  vegetation 
appearing  to  be  in  various  stages  of  succession  were  abundant: 
sawgrass  patches,  sawgrass  patches  with  several  shrub 
individuals  at  their  centers,  brush  patches  fringed  with 
sawgrass,  brush  patches  fringed  with  sawgrass  with  one  or  two 
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trees  at  their  centers,  etc.   Yet  the  process  of  succession  was 
too  slow  to  perceive  in  real  time. 

The  comparison  of  photographs  taken  in  1952  with  satellite 
imagery  from  1990  showed  vegetation  over  a  38  year  time  span, 
long  enough  to  detect  succession.   The  vegetation  changes 
observed  over  this  time  span  reinforced  the  existing  evidence 
that  hydrology  drives  succession  in  the  Everglades.   Where 
hydroperiod  had  increased  over  the  38  year  period,  thousands  of 
hectares  of  sawgrass,  brush  and  trees  were  lost.   Where 
hydroperiod  had  decreased,  hundreds  of  hectares  of  brush  had 
appeared.   In  the  areas  where  hydroperiod  had  not  greatly 
changed,  vegetation  changes  could  not  be  measured,  because  the 
amount  of  error  in  matching  the  photographs  and  satellite  image 
was  as  large  as  the  changes  in  vegetation. 

Because  succession  in  the  non-degraded  parts  of  Water 
Conservation  Area  1  could  not  be  measured  directly,  a  different 
approach  to  the  subject  of  island  growth  was  taken.   The 
distribution  of  tree  islands  of  different  sizes,  as  recorded  in 
a  1987  classified  satellite  image  was  examined.   The  satellite 
image  showed  many  more  small  islands  than  large  ones.   The 
frequency  of  islands  of  each  size  class  decreased  exponentially 
as  their  sizes  increased.   When  the  distribution  of  the  size 
classes  of  the  islands  was  examined,  it  was  found  that  the  size 
classes  of  small  islands  were  separated  by  the  smallest  possible 
increment,  one  pixel.   However,  when  islands  grew  above  a 
certain  size  (1.59  hectares)  the  difference  between  one  island's 
size  and  the  next  larger  island's  size  increased  with  island 
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size,  following  an  exponential  curve.   When  island  sizes  reached 
another  threshold  (2.18  hectares)  the  increase  in  the  difference 
between  one  island  and  the  next  larger  island  followed  a  steeper 
curve,  creating  a  break  or  discontinuity  in  the  curve  of  the 
distribution. 

The  size  and  frequency  distributions  of  tree  islands 
recorded  by  the  satellite  image  raised  three  questions.   Could 
the  ongoing  formation  and  growth  of  islands  create  a  size 
distribution  like  that  found  in  the  landscape  of  Water 
Conservation  Area  1?   Why  were  more  small  islands  than  large 
islands  found  in  the  satellite  image  of  the  landscape?   What 
could  cause  the  discontinuity  in  the  distribution  of  island 
sizes? 

The  first  question  was  addressed  with  simple  models  built 
to  generate  island  size  and  frequency  distributions  for 
different  growth  functions.   The  modeling  exercise  demonstrated 
that  the  creation  and  growth  of  islands  over  time  could  create  a 
distribution  like  that  found  in  the  landscape  of  Water 
Conservation  Area  1.   A  growth  function  dependent  upon  island 
area  produced  the  distribution  most  similar  to  that  observed  in 
the  landscape. 

The  second  question  was  also  addressed  when  the  models 
demonstrated  that  the  large  numbers  of  small  islands  could  be 
the  result  of  a  growth  process  that  changed  with  island  size  or 
scale.   Because  the  metric  used  for  measuring  islands  has  a 
coarser  resolution  than  the  smallest  scale  of  the  process,  many 
small  islands  of  different  sizes  were  grouped  together  in  size 
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classes  separated  by  the  smallest  unit  of  the  metric.   The 
result  was  a  frequency  distribution  by  size  similar  to  that 
measured  in  the  satellite  image  of  the  landscape. 

The  third  question  was  not  easily  answered.   Theories  of 
the  relationship  between  ecosystem  processes  and  ecosystem 
architecture  suggested  that  a  change  of  some  process  affecting 
the  system  could  cause  a  discontinuity  in  the  island  size 
distribution. 

A  conceptual  model  proposing  a  mechanism  for  the 
discontinuity  in  the  island  size  distribution  was  constructed. 
The  hypothetical  system  was  devised  from  a  synthesis  of 
information  on  the  Everglades  ecosystem,  the  Okefenokee  swamp 
ecosystem,  and  fire  behavior.   The  conceptual  model  proposed 
that  in  the  absence  of  fire,  vegetation  follows  a  successional 
sequence  from  wet  prairie  to  trees,  and  that  fire  damage  moves 
vegetation  back  along  the  sequence.   By  the  mechanism  of 
succession,  accelerated  by  peat  accumulation  along  the  margins 
of  incipient  and  developed  islands,  tree  islands  grow  in  area 
over  time.   Fire  damage  sets  back  the  process  of  island  growth. 
The  cores  of  large  islands  develop  a  wetter  microclimate  than 
the  rest  of  the  landscape,  and  so  resist  fire.   Once  islands 
grow  large  enough  for  this  microclimate  effect  to  operate,  their 
growth  rate  is  boosted  by  a  reduction  in  fire  damage.   Island 
growth  and  resistance  to  fire  form  a  positive  feedback  loop, 
producing  a  steep  rise  in  growth  rate,  appearing  as  a 
discontinuity  in  the  island  size  distribution. 
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The  conceptual  model  was  converted  into  a  spatial  model 
capable  of  producing  landscape  maps.   The  spatial  model  was  used 
to  simulate  landscapes  under  different  fire  frequency  and 
severity  regimes.   Changing  the  strength  and  frequency  of  fire 
changed  the  patterns  of  the  modeled  landscapes,  as  measured  by 
the  vegetation  composition  and  the  number  of  vegetation  patches. 
In  the  abstract  ecosystem  of  the  model,  fire  controlled  the 
number  and  sizes  of  tree  islands. 

The  end  product  of  this  investigation  of  the  landscape  of 
Water  Conservation  Area  1--  the  synthesis  of  existing  research 
pertaining  to  the  ecosystem,  the  collection  and  analysis  of 
landscape  data,  the  simulation  of  different  island  growth 
functions,  the  creation  of  a  conceptual  landscape  model,  and  the 
spatial  simulation  of  landscape  patterns--  is  a  theory  of 
landscape  dynamics.   This  theory  is  depicted  in  the  conceptual 
model  in  chapter  4 .   The  theory  has  not  been  proven  or  disproven 
here.   It  may  never  be  possible  to  test  this  theory  the  way  a 
hypothesis  of  the  behavior  of  a  single  species  or  a  small  set  of 
variables  could  be  tested.   The  inability  to  test  hypotheses  is 
a  common  critique  of  ecosystem  and  landscape  ecology. 

However,  theories  addressing  large-scale  questions  can 
only  be  formulated  by  thinking  at  a  large  scale.   Community- 
level  or  species-level  hypotheses  are  more  easily  tested,  but 
the  answers  to  individual  small-scale  questions  cannot  answer 
large-scale  questions.   Many  small-scale  answers  must  be  fit 
together  to  create  a  larger  picture.   Once  a  landscape-scale 
theory  is  devised,  pieces  of  it  can  then  be  tested  using 
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accepted  scientific  methods,  but  it  cannot  be  devised  without 
holistic  thinking. 

This  investigation  has  raised  questions  that  need  to  be 
addressed  to  fill  in  missing  pieces  of  the  picture.   One  of 
these  missing  pieces  is  the  actual  rotation  time  of  fires  in 
Water  Conservation  Area  1 .   This  cannot  be  known  until  detailed 
fire  records  are  kept  for  Water  Conservation  Area  1 .   It  is 
recommended  that  the  time,  location,  and  severity  of  all  fires 
in  the  area  be  recorded,  so  that  future  researchers  can  evaluate 
fire  frequency. 

Another  unanswered  question  concerns  the  formation  of  tree 
islands,  particularly  the  role  of  peat  batteries  and  alligators 
in  the  creation  and  maintenance  of  these  landf orms .   A  long-term 
study  tracking  peat  batteries  is  needed.   Studies  of  alligator 
behavior  with  respect  to  the  site-selection  and  excavation  of 
alligator  holes,  and  the  peat-transporting  effects  of  alligator 
activities  are  needed. 

A  more  difficult  missing  piece  of  the  picture  is  the 
effect  of  sheet  flow  on  the  landscape  of  Water  Conservation  Area 
1.   While  we  can  no  longer  observe  sheet  flow  in  action,  perhaps 
new  technology  applied  to  the  observation  of  tree  island  shapes 
over  time  can  yield  clues  to  its  effect.   If  sheet  flow  was 
responsible  for  the  elongate  shape  of  large  tree  islands,  as  is 
frequently  assumed,  their  shapes  may  be  changing  in  the  absence 
of  sheet  flow.   And  finally,  the  succession  of  vegetation 
leading  to  the  expansion  of  tree  island  margins  (induced  and 
theorized  in  this  work)  remains  to  be  tested. 
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Since  science  has  become  increasingly  reductionist  over 
time,  holistic  approaches  have  been  rejected  rather  than 
improved.   Because  ecologists  have  avoided  rather  than  developed 
holistic  approaches  to  large-scale  questions,  resource  managers 
are  now  being  forced  to  devise  ecosystem  management  strategies 
without  the  benefit  of  an  established  scientific  methodology. 
Prejudice  against  holistic  thinking  among  ecologists  must  be 
overcome  if  the  needed  theories  and  methods  are  to  be  developed. 
This  work  is  a  contribution  toward  the  revival  of  holistic 
thinking,  and  it  is  hoped  that  this  depiction  of  the  ecosystem 
of  Water  Conservation  Area  1  will  be  useful  to  the  ecosystem's 
managers . 
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